Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,434,387 C→A 16.5% intergenic (‑166/‑792) RCAP_rcc01334 ← / → RCAP_rcc01335 hypothetical protein/Fis family GAF modulated sigma54 specific transcriptional regulator

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,434,3870C→A16.5% 149.0 ‑1.3 0.0376intergenic (‑166/‑792)RCAP_rcc01334/RCAP_rcc01335hypothetical protein/Fis family GAF modulated sigma54 specific transcriptional regulator
Reads supporting (aligned to +/- strand):  new base (37/26):  ref base (158/152):  total (195/181)
Fisher's exact test strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 1.00e+00

Only 100 of 396 total aligned reads displayed.

TATATGCACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTC  >  NC_014034/1434287‑1434487
                                                                                                      |                                                                                                     
tataTGCACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTcc                                                                                                       <  2:2545748/101‑1 (MQ=255)
tataTGCACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTcc                                                                                                       <  1:4117648/101‑1 (MQ=255)
 atatGCACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCg                                                                                                      >  1:2799916/1‑101 (MQ=255)
 atatGCACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCg                                                                                                      >  2:3275198/1‑101 (MQ=255)
 atatGCACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCg                                                                                                      <  1:5229816/101‑1 (MQ=255)
  tatGCATCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGc                                                                                                     <  1:941067/101‑1 (MQ=255)
  tatGCACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGc                                                                                                     >  1:5552557/1‑101 (MQ=255)
      cACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTg                                                                                                 <  1:4637778/101‑1 (MQ=255)
      cACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTg                                                                                                 <  1:3735551/101‑1 (MQ=255)
        ccGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCt                                                                                               <  1:1170238/101‑1 (MQ=255)
         ctcTTACCGTGGATGGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTTGGTCCCTCCCCCCGCCGGGGCGTTGATCGATAC..GTCTCATCTTTGCtt                                                                                              >  2:4399931/3‑101 (MQ=255)
            ttACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCAc                                                                                           >  2:1340524/1‑101 (MQ=255)
            ttACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCAc                                                                                           <  1:3995735/101‑1 (MQ=255)
            ttACCGTGGACCGAAACAGGCTGCTTTGGGCGGAACCCCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCAc                                                                                           <  2:709956/101‑1 (MQ=255)
              aCCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACg.g                                                                                        >  1:2921210/1‑101 (MQ=255)
                cGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCt                                                                                      >  2:3385104/1‑101 (MQ=255)
                cGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCt                                                                                      <  1:779832/101‑1 (MQ=255)
                 gTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTg                                                                                     >  2:2831617/1‑101 (MQ=255)
                 gTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTg                                                                                     <  2:861734/101‑1 (MQ=255)
                  tGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGc                                                                                    <  1:1769473/101‑1 (MQ=255)
                   ggATCGAAGCAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCg                                                                                   <  2:5280351/101‑1 (MQ=255)
                   ggATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCg                                                                                   <  2:1717013/101‑1 (MQ=255)
                     aTCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTc                                                                                 >  1:2494070/1‑101 (MQ=255)
                      tCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCa                                                                                <  2:5628684/101‑1 (MQ=255)
                      tCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCGCTGCTTCACG.GCTGCGTCa                                                                                <  1:2760109/101‑1 (MQ=255)
                       cGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCaa                                                                               <  1:4029698/101‑1 (MQ=255)
                       cGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCaa                                                                               >  2:380158/1‑101 (MQ=255)
                        gAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAAc                                                                              >  1:1046784/1‑101 (MQ=255)
                         aaaCAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAAcc                                                                             >  2:5398693/1‑101 (MQ=255)
                         aaaCAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAAcc                                                                             >  2:5735901/1‑101 (MQ=255)
                            cAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTg                                                                          <  2:1971847/101‑1 (MQ=255)
                            cAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTg                                                                          <  1:354953/101‑1 (MQ=255)
                            cAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGCCGTGGATCGAAAC..ATCTCCGCTCTCCTTCACA.GCTCCGCCAACCTTg                                                                          >  1:228588/1‑101 (MQ=37)
                              ggCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTc                                                                        <  2:3837137/101‑1 (MQ=255)
                               gctgctTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCg                                                                       >  1:4428335/1‑101 (MQ=255)
                               gctgctTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCg                                                                       <  2:4345645/101‑1 (MQ=255)
                               gctgctTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCg                                                                       <  1:1382177/101‑1 (MQ=255)
                               gctgctTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCg                                                                       <  1:2726376/101‑1 (MQ=255)
                                ctgctTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGc                                                                      <  2:2716271/101‑1 (MQ=255)
                                    tttGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTccc                                                                  >  2:2065281/1‑101 (MQ=255)
                                    tttGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTccc                                                                  <  2:1854095/101‑1 (MQ=255)
                                     ttGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTcccc                                                                 <  1:3298909/101‑1 (MQ=255)
                                      tGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTccccc                                                                >  1:19375/1‑101 (MQ=255)
                                      tGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTccccc                                                                >  1:2231096/1‑101 (MQ=255)
                                      tGGGCGGAACACCTGATCCTGCTGGGTAGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTccccc                                                                >  2:4136451/1‑101 (MQ=255)
                                       gggCGGAACACCTGCTCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTcccccc                                                               <  2:598040/101‑1 (MQ=255)
                                       gggCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTcccccc                                                               <  1:2049738/101‑1 (MQ=255)
                                       gggCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTcccccc                                                               <  2:1618609/101‑1 (MQ=255)
                                           ggAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTcccccccgcg                                                           >  2:4877709/1‑97 (MQ=255)
                                           ggAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCgg                                                           >  2:672775/1‑101 (MQ=255)
                                            gAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCggg                                                          <  2:808406/101‑1 (MQ=255)
                                             aaCACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCgggg                                                         <  2:286935/101‑1 (MQ=255)
                                              acacCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCggggg                                                        >  2:220893/1‑101 (MQ=255)
                                              acacCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCggggg                                                        >  1:2603282/1‑101 (MQ=255)
                                                   tGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTgg                                                   >  2:2128525/1‑101 (MQ=255)
                                                     aTCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGAt                                                 <  1:1602863/101‑1 (MQ=255)
                                                      tCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                <  2:1199827/101‑1 (MQ=255)
                                                       gggctgtcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGGGGa                                                  >  2:4224076/9‑101 (MQ=0)
                                                       gggctatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGa                                                  <  2:1283625/93‑1 (MQ=2)
                                                           tatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                              >  2:3965041/5‑101 (MQ=2)
                                                           tatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                              <  2:975871/97‑1 (MQ=2)
                                                           tatcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGa                                              >  1:4953106/5‑101 (MQ=2)
                                                            atcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGGGGATCGaa                                             >  1:5706284/4‑101 (MQ=2)
                                                             gggTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacac                                         >  1:3080208/1‑101 (MQ=255)
                                                              cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                           >  1:4583036/2‑101 (MQ=2)
                                                              cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                           <  2:5055095/100‑1 (MQ=2)
                                                              cgTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                           >  2:481564/2‑101 (MQ=2)
                                                              cgTCGCCCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAac                                           <  1:5608186/100‑1 (MQ=2)
                                                               gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacg                                          <  1:1650775/101‑2 (MQ=2)
                                                                tCGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacgc                                         >  1:640131/1‑99 (MQ=2)
                                                                 cGCTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacgct                                        >  1:1493945/1‑98 (MQ=2)
                                                                 cGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGGGGATCGAAacaccggt                                     >  1:4285308/1‑97 (MQ=255)
                                                                  gCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGAtt                                    >  1:1831991/1‑101 (MQ=255)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgc                                      >  2:2985961/1‑96 (MQ=2)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgc                                      <  1:3212315/101‑6 (MQ=2)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgc                                      >  2:4391724/1‑96 (MQ=2)
                                                                   cTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTg                                   <  2:2476035/101‑1 (MQ=255)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgct                                     >  2:468022/1‑95 (MQ=2)
                                                                    tCCCCCCGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgct                                     >  2:4479162/1‑95 (MQ=2)
                                                                      cccccGCGGGGGCGTGGATCGAAACGGGTCGAAGACCTGCTTGACGCGCTGGATCTGCGGCGTCGCTCCCCCCGCGGGGGCGTGGATCGAAacgctgctga                                   <  1:4502407/101‑9 (MQ=2)
                                                                       ccccGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATg                               >  2:1933485/1‑101 (MQ=255)
                                                                          cgcgGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGgcg                            >  2:2355870/1‑101 (MQ=255)
                                                                           gcgGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGgcgc                           <  2:5206216/101‑1 (MQ=255)
                                                                            cgGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGcc                          >  2:1010741/1‑101 (MQ=255)
                                                                             gggggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCg                         <  1:3917733/101‑1 (MQ=255)
                                                                             gggggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCg                         >  2:4570945/1‑101 (MQ=255)
                                                                             gggggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCg                         <  1:4379360/101‑1 (MQ=255)
                                                                              ggggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGa                        <  1:3860377/101‑1 (MQ=255)
                                                                                ggCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAgg                      <  2:1067088/101‑1 (MQ=255)
                                                                                   gTGGATCGAACC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGAc                   <  1:2992390/101‑1 (MQ=255)
                                                                                   gTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGAc                   <  1:5342668/101‑1 (MQ=255)
                                                                                     ggATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAg                 >  2:1457480/1‑101 (MQ=255)
                                                                                      gATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAgg                >  1:4817199/1‑101 (MQ=255)
                                                                                       aTCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGc               >  1:5663459/1‑101 (MQ=255)
                                                                                         cGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCtt             >  2:153626/1‑101 (MQ=255)
                                                                                             aC..GTCTCCGCTCCGATTCACG.GCTGCATCAGCCTTGTCGCTCCCACCGCGGGGGCGTGGATCGAAACTCGGATTGCATGGCGCCGAGGGACAGGCTTTGCa         <  2:2259302/101‑1 (MQ=17)
                                                                                              c..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAg        >  1:3259286/1‑101 (MQ=255)
                                                                                                  tctcCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTc      <  1:3385309/101‑1 (MQ=255)
                                                                                                   ctcCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCg     <  2:2783174/101‑1 (MQ=255)
                                                                                                   ctcCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCg     >  1:501010/1‑101 (MQ=255)
                                                                                                      |                                                                                                     
TATATGCACCGCTTACCGTGGATCGAAACAGGCTGCTTTGGGCGGAACACCTGATCCTGCTGGGTCGCTCCCCCCGCGGGGGCGTGGATCGAAAC..GTCTCCGCTCTGCTTCACG.GCTGCGTCAACCTTGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACACGGATTGCATGGCGCCGAGGGACAGGCTTTGCAGTCGCTC  >  NC_014034/1434287‑1434487

Base quality scores:  ATCG  < 3 ≤  ATCG  < 19 ≤  ATCG  < 27 ≤  ATCG  < 33 ≤  ATCG  < 40 ≤  ATCG