Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 11.1% intergenic (+46/+56) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,6950CG11.1% 1938.9 / 68.7 626intergenic (+46/+56)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (325/221);  new base G (9/69);  total (334/292)
Fisher's exact test for biased strand distribution p-value = 1.73e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 652 total aligned reads displayed.

CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608546‑3608844
                                                                                                                                                      |                                                                                                                                                     
ccacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                       <  2:7073235/150‑1 (MQ=255)
ccacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                       <  1:3288495/150‑1 (MQ=255)
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                      >  2:7232288/1‑150 (MQ=255)
  acGCTTTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGAGAATGCGTTAAGATTGTTGCTCAGGATATAAGGCAAACCTATCGAAAAA‑TTTAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                     <  1:8502329/150‑1 (MQ=255)
   cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  1:7391377/150‑1 (MQ=255)
   cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGACCCCTCGCCCCTCTgg                                                                                                                                                    <  2:9316709/150‑1 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   >  2:8491327/1‑150 (MQ=255)
    gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  1:10740784/150‑1 (MQ=255)
     cTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                  <  1:2979529/150‑1 (MQ=255)
      tGTTAGGGATGCCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                 <  1:10875828/150‑1 (MQ=255)
      tGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                 <  2:3974197/150‑1 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  1:1167865/150‑1 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               >  1:5784610/1‑150 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               >  1:9604895/1‑150 (MQ=255)
         tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  2:3414524/150‑1 (MQ=255)
         tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCGGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAACCCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  2:10466453/150‑1 (MQ=255)
           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                            <  2:123081/150‑1 (MQ=255)
                gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                       <  2:5580371/150‑1 (MQ=255)
                gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                       <  1:6331883/150‑1 (MQ=255)
                gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                       <  2:8392003/150‑1 (MQ=255)
                  ccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                     >  1:6902636/1‑150 (MQ=255)
                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                    >  2:8149778/1‑150 (MQ=255)
                    gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   >  2:7352552/1‑150 (MQ=255)
                      gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 >  2:8737409/1‑150 (MQ=255)
                           ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATATGGCTAACCGTTCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGaaa                                                                                                                            <  1:8420110/150‑1 (MQ=255)
                            gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                           <  2:9564776/150‑1 (MQ=255)
                             ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                          >  1:3751951/1‑150 (MQ=255)
                                     gCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCTTTAAGATTGTTGCGCTGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAAGGCGGCATc                                                                                                                  <  2:7399888/150‑1 (MQ=255)
                                           aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGAGTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGcc                                                                                                            >  2:6070415/1‑150 (MQ=255)
                                                   gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                    >  1:8865140/1‑150 (MQ=255)
                                                   gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                    >  2:4987709/1‑150 (MQ=255)
                                                   gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGAGGAGAGGGTTAGGGTTAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                    >  2:10178713/1‑150 (MQ=255)
                                                    ggCGCCGGTGCTGGTGACCGCGATTGGGTTAAGTTTGTTGCGCTGGTGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGTTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                   <  2:11040620/150‑1 (MQ=38)
                                                    ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGACATCGAAGCCAATCAGccc                                                                                                   >  2:8026303/1‑150 (MQ=255)
                                                      cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                 >  2:3081039/1‑150 (MQ=255)
                                                      cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGGTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                 >  1:4173254/1‑150 (MQ=255)
                                                       gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGAGAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCACGCAGCCCCTa                                                                                                >  2:2458849/1‑150 (MQ=39)
                                                         cGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATc                                                                                              >  1:9999646/1‑150 (MQ=255)
                                                              cTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTCTCAACCg                                                                                         >  1:7999349/1‑150 (MQ=255)
                                                              cTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCg                                                                                         >  2:1881680/1‑150 (MQ=255)
                                                               tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                        >  2:10328850/1‑150 (MQ=255)
                                                                 gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                      >  1:10305767/1‑150 (MQ=255)
                                                                  tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                     >  1:6768764/1‑150 (MQ=255)
                                                                      aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                 >  1:1100347/1‑150 (MQ=255)
                                                                        cgtaTGGGTTAAGTTTGTTGCGCAGGATATCCGGCACCCCTATCGCAACT‑GTGAGCGCTAGCGGGCCCCTCGCCCCTCTGGGTCGGGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                               <  2:979255/147‑1 (MQ=255)
                                                                        cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                               >  1:7310921/1‑150 (MQ=255)
                                                                         gAATGCGTTATGATTGTTGCGCAGGAGATAAGGCAACCCGATCGCAACT‑TTGCGCGCGATCGGACCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                              <  2:1796362/150‑1 (MQ=38)
                                                                              cGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TGGAGCGCGAGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGACGCCCATCAGCCCCTATCAACCGCCTTTACGAATCAAAt                                                                         >  1:8115914/1‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  1:5602124/1‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  2:3766123/1‑150 (MQ=255)
                                                                                ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                       >  1:10126848/1‑150 (MQ=255)
                                                                                ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                       >  1:4577586/1‑150 (MQ=255)
                                                                                 tctgcttgttgCGCCGGAGATAAGGCAACCCGATCG‑AACATTTGCGCGCGATCGGCCCCCCCGCCCCTGTGGGGTGTGTTTTTGGGGGTGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                      <  1:3507449/145‑1 (MQ=255)
                                                                                   attttgttgCGCTGGATTTACGGCAACCCTATGGAAACT‑TTGCGCGCGATGGGCCCCCTCCCCCCGGTGGGGTGTGGTTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                    <  2:4059036/147‑1 (MQ=255)
                                                                                    gATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACTAt                                                                   >  1:10000428/1‑150 (MQ=255)
                                                                                      ttgttcCCCAGCACATAAGCCAAACCCAACCCAACA‑TTCACCGCGAACGGACCCCACCCCCCTCTGGGCACAGGCTTAGGCTCAGGGCAAAACGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                 <  2:3518668/150‑1 (MQ=255)
                                                                                      tggttgCGATGGATATACTGCAAACCTCTCGCACCT‑TTGCGCGGGCTCGGCCCCCGGCCCCGTGTGGGGTGGGGTTTGGGGGTTGGTCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                 <  1:4761861/148‑1 (MQ=255)
                                                                                          tgCGCAGGCTATACGGAAACCCGATCTCAACT‑TTGAGCGGGATCGGTCCCCTCGCCCCTCTGGTTTTTGGTTTTGGTTTTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCa                                                             <  1:7137017/150‑1 (MQ=255)
                                                                                              ctggAGATAAGGCAACCCGTTCGCAACT‑TGGCGCGCGAGCGGTCCCCTCGCCCCTCTGGGGAGTGGTTTGGGTTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTcc                                                         <  1:5315061/148‑1 (MQ=255)
                                                                                                gtttatatGGCAACCCTATCGCAACT‑TTGCGCGCGTTGGGTCCCCGCGCCCCTCTGGGGTGTGGTTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAt                                                       <  1:5977747/146‑1 (MQ=255)
                                                                                                ggtgaTAAGGCAAACCGATCGCACCT‑TTGTGCGGTATCGGCCGCCCCGCCCCTCTGGTGAGTGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAt                                                       <  2:5805388/147‑1 (MQ=255)
                                                                                                ggAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCAt                                                       >  2:3232609/1‑150 (MQ=255)
                                                                                                   gaTAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCg                                                    >  1:980926/1‑150 (MQ=255)
                                                                                                      acggCAACCCTCTCGCAAAT‑TTGAGCGCGCTCGGCCCCCTCCCCCCTGGGGGGATTGGTTTTGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 <  1:7888982/148‑1 (MQ=255)
                                                                                                      acgccAAACCGATCCCACCT‑GTGCGCGCGTCGGGCCCCCGCCCCCCTGGGGGGTGTGGGTTCGGGGGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 <  1:3885782/146‑1 (MQ=255)
                                                                                                       aGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtctc                                                >  1:10422284/1‑150 (MQ=255)
                                                                                                        ggCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                               >  2:9079610/1‑150 (MQ=255)
                                                                                                         gCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGACTCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  2:10617876/1‑150 (MQ=255)
                                                                                                          cataCCGATCGCAACT‑TTGAGCGCGATCGGTCCCCTCGCCCGTCTGGGTAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             <  2:6445145/147‑1 (MQ=255)
                                                                                                          cAACCCGCTCGCAACA‑TTGAGCGCGCTCGGTCCCCCCGCCCCTCGGGGGAGAGGTTTAGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             <  2:6608851/150‑1 (MQ=255)
                                                                                                          cAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             >  1:2102911/1‑150 (MQ=255)
                                                                                                           aacccGAGCGAAACA‑TGGCGCGGGAGCGGTCCCCCCGCCCCGCTGGGGAGTGGTTTGGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            <  2:7879741/147‑1 (MQ=255)
                                                                                                           aaaCCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            >  2:4280430/1‑150 (MQ=255)
                                                                                                              ccGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGTGTGAGGGGAAAAGGAGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCa                                         >  2:2315989/1‑150 (MQ=255)
                                                                                                                 aTCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                      >  1:2511311/1‑150 (MQ=255)
                                                                                                                 aTCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                      >  1:8737396/1‑150 (MQ=255)
                                                                                                                  tCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                     <  2:2708442/150‑1 (MQ=255)
                                                                                                                     cAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTGTAGCAACCAGtt                                  >  2:1137738/1‑150 (MQ=255)
                                                                                                                      aaCA‑TTGCGCGCTACCGGACCCCCCCCCCCTCTGGGCAGAGTTTTAGGGCCAGGCCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                 <  1:1228339/150‑1 (MQ=255)
                                                                                                                        cA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                               <  2:6356472/150‑1 (MQ=255)
                                                                                                                               gcgcGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa                         >  1:2645376/1‑150 (MQ=255)
                                                                                                                                cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                        >  1:5046901/1‑150 (MQ=255)
                                                                                                                                cgcgATCGGTCCCCGCGCCCCTGGGGGTTGAGGGTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                        <  1:4349900/150‑1 (MQ=255)
                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                     >  2:6892483/1‑150 (MQ=255)
                                                                                                                                   gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                     >  2:2084489/1‑150 (MQ=255)
                                                                                                                                    cgcGGGCCCCGCGCCCCTGGGGGGTGGGGGTGCGGGTGCGGGGACCAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    <  2:8929533/148‑1 (MQ=255)
                                                                                                                                        ggccccGCGCCCCGGTGGGGTGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                <  2:1582940/148‑1 (MQ=255)
                                                                                                                                        gTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                >  2:2888224/1‑150 (MQ=255)
                                                                                                                                         tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa               >  1:1637256/1‑150 (MQ=255)
                                                                                                                                          acccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa              >  2:5375162/2‑150 (MQ=255)
                                                                                                                                             cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc           >  1:9269516/1‑150 (MQ=255)
                                                                                                                                             cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc           >  2:7720264/1‑150 (MQ=255)
                                                                                                                                               cGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc         <  2:10130893/150‑1 (MQ=255)
                                                                                                                                                gCCCCTCTGGGGTGTGGTTTAGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  1:1318943/150‑1 (MQ=255)
                                                                                                                                                gCCCCTCTGGGGAGAGGGTTCGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  2:1319834/150‑1 (MQ=255)
                                                                                                                                                gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        >  2:5943587/1‑150 (MQ=255)
                                                                                                                                                 ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg       <  1:8207688/150‑1 (MQ=255)
                                                                                                                                                  cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg      >  2:4676125/1‑150 (MQ=255)
                                                                                                                                                     tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat   >  1:2838595/1‑150 (MQ=255)
                                                                                                                                                     tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat   >  2:9182124/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACA‑TTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608546‑3608844

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑