Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,721 A→C 6.1% intergenic (+72/+30) 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,7210AC 6.1% 1671.1 / 19.0 553intergenic (+72/+30)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base A (342/164);  new base C (6/41);  total (348/205)
Fisher's exact test for biased strand distribution p-value = 1.50e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 574 total aligned reads displayed.

TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGC  >  NC_000913/3608572‑3608870
                                                                                                                                                     |                                                                                                                                                     
tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGaa                                                                                                                                                       >  2:2534770/1‑150 (MQ=255)
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                     <  1:5437298/150‑1 (MQ=255)
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                     >  1:7919049/1‑150 (MQ=255)
    tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                   >  1:2020494/1‑150 (MQ=255)
      gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                                                 >  2:931400/1‑150 (MQ=255)
         gtggtgGCATATAGGGGGGCGACGTTGCTGGTGACAGCTATTGGGTTAATTTTTTTGCGCAGGATATATGGCAAACCTATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGTGAGTGTGTTTGGGTGAGGGGAAAAGGCGGCa                                                                                                                                              <  2:5915917/150‑1 (MQ=255)
          tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                                             >  2:8587474/1‑150 (MQ=255)
             tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGAGAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGACAAGGCGGCATCGa                                                                                                                                          >  2:8067487/1‑150 (MQ=255)
                 aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGcc                                                                                                                                      >  1:8184841/1‑150 (MQ=255)
                   aTACGGGGGCGAAGGTGCTGGTGACAGCGAATGCGTTAAGAATGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCGTCTGGGGAGAGGGTTAGCGAGAGGGGAAAAGGCGGCATCGAAGCCaa                                                                                                                                    >  2:7377642/1‑150 (MQ=255)
                     aCGGGGGCGACGGTGCTGGTGACAGCGAATGCTTTAATATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGAGCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATc                                                                                                                                  <  2:1529205/150‑1 (MQ=255)
                      cGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                                                 >  1:7349673/1‑150 (MQ=255)
                        ggggCGTCGGTGGTGGTGACAGCGAATGCGTTATGATGGTTGCGCCGGAGATACGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                                               <  2:6243998/150‑1 (MQ=255)
                         gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                                              >  1:207874/1‑150 (MQ=255)
                         gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGGTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                                              >  2:2729258/1‑150 (MQ=255)
                          ggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                                             >  1:3896300/1‑150 (MQ=255)
                           gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                                                                                            >  2:8448983/1‑150 (MQ=255)
                            cGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCt                                                                                                                           >  2:8359695/1‑150 (MQ=255)
                              ccGGTGCTGGGGCCCGCGATTGCGTTATTTTTGTTGCGATGGATTTATGGAAAACCGATCGCAACTTTGCGCGCGTTCGGTCCCCTCGCCCCTGTGGGTAGTGGGTTAGGGGGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAt                                                                                                                         <  2:22292/149‑1 (MQ=255)
                                ggTGCTGGTGACAGCGACGGCGTAAAGATTGTTGCCCTGGAGATACGCAAAACCGATCGCAACTTGGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGCGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCa                                                                                                                       <  2:6725095/150‑1 (MQ=255)
                                  ggggggTGACAGCGATTGCGTTATGTTTGTGGCGCTGGTGATAGGGAAACCCGATCGCAACATTGAGCGGGAGGGGGCCCCTCGCCCGTGTGGGGTGTGGGTTTGGTTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                     <  2:4519640/146‑1 (MQ=255)
                                    cTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCg                                                                                                                   >  1:4139373/1‑150 (MQ=255)
                                     tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                  >  1:9672688/1‑150 (MQ=255)
                                     tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                  >  2:4171400/1‑150 (MQ=255)
                                     tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                  >  1:2044104/1‑150 (MQ=255)
                                     tGGTGACAGCGAATGCGTGAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                                                  >  1:4837806/1‑150 (MQ=255)
                                        tGACAGCTATTGGGTTAAGATTGTTGCGCAGGAGATACGGAAAACCTTTCGCAACTTTGAGCGCGTTGGGTCCCCTGGCCCCTGTGGGTTGTGGTTTCGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                                               <  2:5140345/150‑1 (MQ=255)
                                        tGACAGCGAATGCGTTAAGATGGTTGCGCCGGTGTTATGGCAACCCGATCGCAACATTGAGCGCGATCGGTCCCCCCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                                               <  2:6485836/150‑1 (MQ=38)
                                         gACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCttt                                                                                                              >  1:6226354/1‑150 (MQ=255)
                                          aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtttc                                                                                                             >  2:2957159/1‑149 (MQ=255)
                                          aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                                             >  1:6765978/1‑150 (MQ=255)
                                           cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATACGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                            <  2:7250183/150‑1 (MQ=255)
                                            aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                           >  2:3822990/1‑150 (MQ=255)
                                            aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                           >  1:64472/1‑150 (MQ=255)
                                             gCGAATGGGTTAAGTTTGTTGCGCAGGAGATAGGGCAACCGGAGCGCAACATTGAGCGGGATCGGTCCCCTCGCCCCTCGGGGGTGTGGGTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGa                                                                                                          <  2:638536/150‑1 (MQ=255)
                                              cGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGaa                                                                                                         >  1:2881045/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                  >  1:7193647/1‑150 (MQ=255)
                                                      ttatttttgttgCGCTGTATATATGGCAAACCGTGGGCAACTTTGCGGGCGTTGGGCCCCCTCGCCCCTGTGGGGTGTGGTTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                                 <  2:7756958/144‑1 (MQ=255)
                                                      ttATGATTGTTGCGAAGGAGATACGGCAACCCGAGCGCAACATTGAGCGCGAGCGGTCCCCGCGCCCCTCTGGGGAGCGGGTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                                 <  2:7577140/150‑1 (MQ=255)
                                                      ttAAGATTGTTGCGCAGGAGATATGGCAAACCGATCGCAACTTTGAGCGCGATGGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                                 <  2:4906723/150‑1 (MQ=39)
                                                       tAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                                                >  1:5382677/1‑150 (MQ=255)
                                                           aTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                                            >  1:3243072/1‑150 (MQ=255)
                                                                  cgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                                                     >  1:8560059/1‑150 (MQ=255)
                                                                     aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                                                  >  1:6190548/1‑150 (MQ=255)
                                                                       gagaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                                                >  1:594647/1‑150 (MQ=255)
                                                                       gagaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                                                >  2:5344770/1‑150 (MQ=255)
                                                                       gagaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                                                >  1:1232437/1‑150 (MQ=255)
                                                                          aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGGGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                                             >  2:6143546/1‑150 (MQ=255)
                                                                          aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTATGGGGAGAGGGTTCGGTTGAGCGGAAAAGCCGGCATCGAAGCCAAAATGCCCCTATCAACCGCCTTAACGAATCAACTAACGATAAGGCAGTCCATCCGt                                                                             >  2:5527151/1‑150 (MQ=255)
                                                                           tAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                                            >  1:8316574/1‑150 (MQ=255)
                                                                               gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAAGGCGGCAGCGAAGCCCATCAGCCCCTAGCAACCGCCTGTACGCATCAAATAACGATACGGCCGGCCCTCCGtctcgt                                                                        >  2:9415348/1‑148 (MQ=37)
                                                                               gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                        >  1:9477556/1‑150 (MQ=255)
                                                                               gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                                                        >  2:6848130/1‑150 (MQ=255)
                                                                                 acaCCGATCGCAACATTGTGCGCGATGGGTCCCCTCGCCCCTGTGGGGTGTGGTTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                      <  2:1117987/148‑1 (MQ=255)
                                                                                 aaaCCTATCGCAACTTTGAGCGCGATCGGCCCCCTCGCCCCTCGGGGTAGTGGGTTTGGGTGTGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                      <  2:605519/150‑1 (MQ=255)
                                                                                     cGATCGCAACATTGAGCGCGATCGGTCCCCTCCCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                                                  >  1:3508573/1‑150 (MQ=255)
                                                                                      gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                 >  1:2600758/1‑150 (MQ=255)
                                                                                         cGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAg                                                              >  2:4402175/1‑150 (MQ=255)
                                                                                            aaCATTGAGCGCGAGCGGTCCCCTCGCCCCTCGGGGGAGAGGGTGAGGGTGAGGGGAAAAGGCGGCATCGAAGCCCATCAGCCCCTCTCCACCGCCTTTACGAATCAAATAACGAGAAGGCCGTCCATCCGTCTCGTGAGCACCCAGtgc                                                           >  2:1228328/1‑148 (MQ=255)
                                                                                                 tgcgcgcGTTGGGCCCCCCCCCCCCTCGGGGGTGTGGGTTCGGTTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                                                      <  1:1720581/147‑1 (MQ=255)
                                                                                                 ggAGCGCGATCGGTCCCCCCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGtt                                                      <  1:1030635/149‑1 (MQ=255)
                                                                                                  gcgcgcgcTCGGGCCCCGCGCCCCGCGGGGGAGGGGTTTCGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                                                     <  2:2900268/148‑1 (MQ=255)
                                                                                                     cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGACCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                                                  >  1:6959960/1‑150 (MQ=255)
                                                                                                        gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                               >  2:1918665/1‑150 (MQ=255)
                                                                                                        gATCGGTCCCCTCGCCCCTATGGTGAGTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                               <  1:9580911/150‑1 (MQ=255)
                                                                                                            ggTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa                                           <  1:7822300/150‑1 (MQ=255)
                                                                                                            ggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCAGTCTCTTTAGCAACCAGTTCGAGTTCCATAAAGGTa                                           >  1:4495540/1‑150 (MQ=255)
                                                                                                              tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTCCGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                                         >  1:5751046/1‑150 (MQ=255)
                                                                                                              tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa                                         >  1:151310/1‑150 (MQ=255)
                                                                                                               ccccTCGCCCCTGTGGGGTGTGGGTTTGGTTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa                                        <  1:3775753/150‑1 (MQ=255)
                                                                                                                cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                       >  2:8638538/1‑150 (MQ=255)
                                                                                                                cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                       >  2:1747936/1‑150 (MQ=255)
                                                                                                                cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                       >  1:9401017/1‑150 (MQ=255)
                                                                                                                  cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                     >  2:6582703/1‑150 (MQ=255)
                                                                                                                  cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                     >  2:3245676/1‑150 (MQ=255)
                                                                                                                  cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                     >  2:9250597/1‑150 (MQ=255)
                                                                                                                     gCCCCTCTGGTGTGTGGTGTTGGGGGAGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa                                  <  1:8944320/150‑1 (MQ=255)
                                                                                                                     gCCCCTCTGGGTTGTGGTTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa                                  <  2:7963836/150‑1 (MQ=255)
                                                                                                                      ccccTGTGGGGAGTGGTTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg                                 <  1:5176631/150‑1 (MQ=255)
                                                                                                                       cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg                                >  2:336591/1‑150 (MQ=255)
                                                                                                                       cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAgg                                >  1:1650677/1‑150 (MQ=255)
                                                                                                                          tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                             >  1:2283/1‑150 (MQ=255)
                                                                                                                             ggggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATc                          >  2:4033307/1‑150 (MQ=255)
                                                                                                                               ggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgc                        <  1:9188386/150‑1 (MQ=255)
                                                                                                                                 tgtgggTTTGGGGGGGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCgg                      <  2:1374945/147‑1 (MQ=255)
                                                                                                                                 agagGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCgg                      >  2:4696257/1‑150 (MQ=255)
                                                                                                                                      gTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGt                 <  1:5569479/150‑1 (MQ=255)
                                                                                                                                       ttAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg                >  1:4988177/1‑150 (MQ=255)
                                                                                                                                       ttAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg                >  2:3601626/1‑150 (MQ=255)
                                                                                                                                           ggTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCg            >  2:9067477/1‑150 (MQ=255)
                                                                                                                                             tGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGa          >  2:2352302/1‑150 (MQ=255)
                                                                                                                                              gtgggtAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAt         <  1:7512330/148‑1 (MQ=255)
                                                                                                                                              gcggggACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATCACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAt         <  1:119374/148‑1 (MQ=255)
                                                                                                                                              gcgggcAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAt         <  1:5303784/148‑1 (MQ=255)
                                                                                                                                               aGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc        >  1:3035819/1‑150 (MQ=255)
                                                                                                                                                 gggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCTCGGGTAGTGGCCGGAtcgt      >  2:8923506/1‑150 (MQ=255)
                                                                                                                                                 gggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgt      <  1:5219471/150‑1 (MQ=255)
                                                                                                                                                 gggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgt      >  1:6580725/1‑150 (MQ=255)
                                                                                                                                                  ggACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGCCACAACTCACGATACGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGGTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGGGGCCGGAtcgtc     <  1:4648603/150‑1 (MQ=255)
                                                                                                                                                  ggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtc     >  1:1610549/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGGC  >  NC_000913/3608572‑3608870

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: