Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,398,027 A→T 1.5% intergenic (+331/+1567) RCAP_rcc01307 → / ← RCAP_rcc01309 hypothetical protein/IS66 family transposase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,398,0270A→T 1.5% 3.0 0.0 ‑0.5264intergenic (+331/+1567)RCAP_rcc01307/RCAP_rcc01309hypothetical protein/IS66 family transposase
Reads supporting (aligned to +/- strand):  new base (2/2):  ref base (132/128):  total (134/130)
Fisher's exact test strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.34e-01

Only 100 of 291 total aligned reads displayed.

ATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGCAC  >  NC_014034/1397927‑1398127
                                                                                                    |                                                                                                       
aTCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAACGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCaagaa                                                                                                         <  1:2807219/101‑1 (MQ=255)
 tCGAAACTCTATGGGGGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCTCCGCGCAAAAAGCCCGCCaagaag                                                                                                        >  1:4878720/1‑101 (MQ=255)
  cGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGc                                                                                                       >  1:2115356/1‑101 (MQ=255)
   gAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCg                                                                                                      >  2:2105395/1‑101 (MQ=255)
    aaaCTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCgg                                                                                                     >  2:3364819/1‑101 (MQ=255)
     aaCTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..a                                                                                                  >  2:4015382/1‑101 (MQ=255)
       ctctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGc                                                                                                <  2:5360576/101‑1 (MQ=255)
        tctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa                                                                                               >  1:4329449/1‑101 (MQ=255)
        tctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa                                                                                               >  2:2542963/1‑101 (MQ=255)
        tctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa                                                                                               >  2:4736212/1‑101 (MQ=255)
        tctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa                                                                                               >  2:4156431/1‑101 (MQ=255)
        tctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa                                                                                               >  2:4422851/1‑101 (MQ=255)
         ctATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCa.a                                                                                             >  2:1977353/1‑101 (MQ=255)
          tATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.Ag                                                                                            <  2:1197715/101‑1 (MQ=255)
          tATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.Ag                                                                                            >  1:3859367/1‑101 (MQ=255)
           aTGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGt                                                                                           <  1:1799828/101‑1 (MQ=255)
           aTGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGt                                                                                           >  1:854999/1‑101 (MQ=255)
            tGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTc                                                                                          >  2:5587591/1‑101 (MQ=255)
            tGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTc                                                                                          <  1:3147309/101‑1 (MQ=255)
                 ggctagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCAAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcc                                                                                     <  1:5352391/101‑1 (MQ=255)
                  gctagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTccc                                                                                    >  2:157567/1‑101 (MQ=255)
                  gctagggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTccc                                                                                    <  2:1939131/101‑1 (MQ=255)
                     agggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTcccccc                                                                                 <  2:584119/101‑1 (MQ=255)
                       ggaggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCcgc                                                                               >  2:5687864/1‑101 (MQ=255)
                         aggctaCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCgg                                                                             >  2:4390677/1‑101 (MQ=255)
                             taCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGc                                                                         >  1:751215/1‑101 (MQ=255)
                             taCGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGc                                                                         >  1:1100677/1‑101 (MQ=255)
                               cGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGt                                                                       >  1:152853/1‑101 (MQ=255)
                                 cACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTgg                                                                     >  1:4719248/1‑101 (MQ=255)
                                   ccTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   >  2:312/1‑101 (MQ=255)
                                   ccTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGAt                                                                   <  2:5667172/101‑1 (MQ=255)
                                    cTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  >  1:4065089/1‑101 (MQ=255)
                                    cTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  >  2:3743021/1‑101 (MQ=255)
                                     ttCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCg                                                                 >  2:2746258/1‑101 (MQ=255)
                                      tgcggcgTCGCTCCCCCCGCGGGGGCGTGGATCGAAA.CGCTGCTGACACCGCCAACAAATATGTGGTCAACACTGTCGCTCCCCCCGCGGGGGCGTGGATc                                                                  <  1:2512478/95‑1 (MQ=2)
                                       cAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAACCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTa                                                               <  2:5221852/101‑1 (MQ=255)
                                       cAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGGGGATCGTa                                                               >  1:2599025/1‑101 (MQ=255)
                                         aaaGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAg                                                             >  1:785133/1‑101 (MQ=255)
                                          aaGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGc                                                            <  1:4656067/101‑1 (MQ=255)
                                            gTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCgg                                                          >  2:3298320/1‑101 (MQ=255)
                                               gCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGc                                                       >  2:2145595/1‑101 (MQ=255)
                                                cTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCa                                                      >  2:5646770/1‑101 (MQ=255)
                                                cTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCa                                                      >  1:3250909/1‑101 (MQ=255)
                                                cTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCa                                                      >  2:2920953/1‑101 (MQ=255)
                                                 tCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCAt                                                     >  1:3967480/1‑101 (MQ=255)
                                                  ccccccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATc                                                    >  2:2456415/1‑101 (MQ=255)
                                                  ccccccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATc                                                    >  1:1319014/1‑101 (MQ=255)
                                                  ccccccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATc                                                    <  1:1949584/101‑1 (MQ=255)
                                                   cccccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCg                                                   <  2:4224928/101‑1 (MQ=255)
                                                   cccccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCg                                                   >  2:1117873/1‑101 (MQ=255)
                                                     cccGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGc                                                 >  2:5083289/1‑101 (MQ=255)
                                                       cgcgGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGccc                                               >  1:504424/1‑101 (MQ=255)
                                                        gcgGGGGCGTGGATCGAACCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCg                                              <  1:1024253/101‑1 (MQ=255)
                                                        gcgGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCg                                              >  1:1745807/1‑101 (MQ=255)
                                                          gggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGAt                                            >  2:3511591/1‑101 (MQ=255)
                                                          gggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGAt                                            >  1:2422232/1‑101 (MQ=255)
                                                           gtgtcGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATa                                           <  2:5531245/97‑1 (MQ=255)
                                                           ggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATa                                           >  1:5456600/1‑101 (MQ=255)
                                                           ggggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATa                                           <  1:5519522/101‑1 (MQ=255)
                                                            gggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATaa                                          <  1:4423123/101‑1 (MQ=255)
                                                            gggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATaa                                          >  2:482868/1‑101 (MQ=255)
                                                             ggCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATaaa                                         <  1:2296165/101‑1 (MQ=255)
                                                              gCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAAc                                        <  1:5647062/101‑1 (MQ=255)
                                                                 tGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAAcccc                                     <  1:5597524/101‑1 (MQ=255)
                                                                  ggATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCg                                    >  1:67425/1‑101 (MQ=255)
                                                                  ggATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCg                                    <  2:3450872/101‑1 (MQ=255)
                                                                  ggATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCg                                    <  2:2972194/101‑1 (MQ=255)
                                                                  ggATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGGGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCg                                    >  2:1845163/1‑101 (MQ=255)
                                                                       gAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGAtttt                               <  1:2088185/101‑1 (MQ=255)
                                                                        aaaCCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTc                              <  1:2168401/101‑1 (MQ=255)
                                                                            cgccgcGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGcc                          >  1:1253022/1‑101 (MQ=255)
                                                                            cgccgcGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGcc                          <  2:4691024/101‑1 (MQ=255)
                                                                             gccgcGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGccc                         <  2:2180743/101‑1 (MQ=255)
                                                                              ccgcGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGcccc                        >  2:1424847/1‑101 (MQ=255)
                                                                               cgcgCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCt                       <  1:4759586/101‑1 (MQ=255)
                                                                               cgcgCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCt                       >  2:357529/1‑101 (MQ=255)
                                                                                gcgcAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTg                      <  2:854930/101‑1 (MQ=255)
                                                                                gcgcAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTg                      <  2:3967236/101‑1 (MQ=255)
                                                                                gcgaAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTg                      <  2:5456310/101‑1 (MQ=255)
                                                                                 cgcAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGa                     >  1:3738749/1‑101 (MQ=255)
                                                                                  gcAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTgcg                    >  1:3910777/1‑99 (MQ=255)
                                                                                   cAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgg                   <  2:679336/101‑1 (MQ=255)
                                                                                   cAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAgg                   >  2:5186830/1‑101 (MQ=255)
                                                                                    aaaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGCggg                  >  2:2585579/1‑101 (MQ=255)
                                                                                    aaaaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAggg                  <  1:1046921/101‑1 (MQ=255)
                                                                                      aaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGa                >  1:2930165/1‑101 (MQ=255)
                                                                                      aaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGa                <  1:2785032/101‑1 (MQ=255)
                                                                                      aaaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGa                >  1:1199097/1‑101 (MQ=255)
                                                                                       aaGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGaa               <  1:2453068/101‑1 (MQ=255)
                                                                                        aGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGaaa              <  2:1523487/101‑1 (MQ=255)
                                                                                         gCCCGCCAAGAAGCGG..AGCG.AGTCGCCCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAg             <  1:3478150/101‑1 (MQ=255)
                                                                                          cccGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGa            <  1:5224210/101‑1 (MQ=255)
                                                                                           ccGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGAt           <  2:1318926/101‑1 (MQ=255)
                                                                                            cGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGAtt          <  1:5140694/101‑1 (MQ=255)
                                                                                             gCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTc         <  1:936842/101‑1 (MQ=255)
                                                                                               cAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCtt       >  1:4221495/1‑101 (MQ=255)
                                                                                               cAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCtt       <  1:1129853/101‑1 (MQ=255)
                                                                                                  gaagCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGc    >  2:4586871/1‑101 (MQ=255)
                                                                                                    agCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGcac  >  1:3097347/1‑101 (MQ=255)
                                                                                                    agCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGcac  >  2:3940353/1‑101 (MQ=255)
                                                                                                    |                                                                                                       
ATCGAAACTCTATGGGAGGCTAGGGAGGCTACGCACCTTCAAAAGTCGCTCCCCCCGCGGGGGCGTGGATCGAAACCGCCGCGCAAAAAGCCCGCCAAGAAGCGG..AGCA.AGTCGCTCCCCCCGCGGGGGCGTGGATCGTAAGCGGTGCATCGGCCCGATAAACCCCGATTTTCGGCCCCTGAGGGGAAAGATTCTTGGCAC  >  NC_014034/1397927‑1398127

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