Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,721 A→C 13.4% 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,7210AC13.4% 342.8 / 11.2 128intergenic (+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 (78/30);  new base C (2/18);  total (80/48)
Fisher's exact test for biased strand distribution p-value = 2.01e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGG  >  NC_000913/3608572‑3608869
                                                                                                                                                      |                                                                                                                                                    
tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                                      >  1:1805995/1‑150 (MQ=255)
tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                                      >  1:510540/1‑150 (MQ=255)
 ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                                                     >  1:2489659/1‑150 (MQ=255)
 ggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                                                     >  1:2539870/1‑150 (MQ=255)
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                    <  2:2473982/150‑1 (MQ=255)
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                    <  2:2598391/150‑1 (MQ=255)
  gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                                    >  1:1820445/1‑150 (MQ=255)
   ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCACCATTGAGCGCGAGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGACAAg                                                                                                                                                   >  2:2879306/1‑150 (MQ=255)
    tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                                                  >  2:3114728/1‑150 (MQ=255)
     ggcAGTGGTGGCAGATTCGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCATCGGGGGATAGGGTTAGGGTGAGGTTAAAAggc                                                                                                                                                 <  2:2252765/150‑1 (MQ=255)
       ctgTGCTGGCATTTAGGGGGGCGAGGGTGGTGGTTACTGCTATGTGTTTTTGTTTTTTGCGCTGGTGTTAAGGCAACCGTATCGAAACTTTTAGCGCGTTGGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGTAAAAggcgg                                                                                                                                               <  2:91333/148‑1 (MQ=255)
       cAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcgg                                                                                                                                               >  2:1062862/1‑150 (MQ=255)
        aGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcggc                                                                                                                                              >  1:1584551/1‑150 (MQ=255)
         gTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCa                                                                                                                                             >  2:1988702/1‑150 (MQ=255)
           gctggctgATACGGGGGCGACGGTGCTGGTGACCGCGATTGCGTTAATATTGTTGCGCAGGAGATATGGCAAACCGATCGCAACATTGAGCGCGTTGGGTCCCCTCGCCCCTCTGGGGATTGGGTTAGGGTGAGGGGAAAAGGCGGCATc                                                                                                                                           <  2:763931/150‑1 (MQ=38)
            cTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCg                                                                                                                                          >  1:2303365/1‑150 (MQ=255)
             tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                                         >  1:1416995/1‑150 (MQ=255)
             tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                                         >  1:1416932/1‑150 (MQ=255)
                ctgttcggggggCGAGGGTGGTGGTGACAGCGAATGCGTTAAGTTTGTTGCGCAGGAGATAAGGCAAACCGATCGAAACATTGAGCGCGATCGCTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAGc                                                                                                                                      <  2:2045113/143‑1 (MQ=38)
                      cGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                                                >  1:2709044/1‑150 (MQ=255)
                      cGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                                                >  1:1886866/1‑150 (MQ=255)
                      cGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCa                                                                                                                                >  1:2709018/1‑150 (MQ=255)
                              cgggTGCTGGTTACAGCTATTGCTTTATGTTTGTTGCGCCGTATATAAGGCAACCCGATCTAAACTTTGAGGGCGATGGGTCCCCTCGCCCCTCTGGGGTGAGGGTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTAt                                                                                                                        <  2:1694954/148‑1 (MQ=25)
                                  tGCTGGTGACAGCCAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                                                    >  2:1875487/1‑150 (MQ=255)
                                   gggggTGCCAGCGCAGGGGTTAATATTGTGGCGCCGTAGATATGGCAACCCGTGCGCAACTTTGAGCGCGTTCGGCCCCCTCGCCCCTGTGGGGTGTGGGTTGGGGGGGGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                                                   <  2:953912/147‑1 (MQ=255)
                                       gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                                               >  1:2296119/1‑150 (MQ=255)
                                       gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                                               >  2:632704/1‑150 (MQ=255)
                                        tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCtt                                                                                                              >  1:2867088/1‑150 (MQ=255)
                                          aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGCTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                                            >  2:1235727/1‑150 (MQ=255)
                                          aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGTCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                                            >  1:1940613/1‑150 (MQ=255)
                                           cAGCGAATGCGTTAAGATGGTTGCGCCGGAGATACGGCAACCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGGGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                                           <  2:2714638/150‑1 (MQ=38)
                                            aGCGAATGCGTTAAGTTTGTTGCGCCGGTGATATGCCAAACCGATCGAAAATTTGAGCGCGTTCGTTCCCCTCGCCCCTCTGGGGTGAGTGTTTGGTTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACg                                                                                                          <  2:1088392/150‑1 (MQ=255)
                                            aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTGTACg                                                                                                          >  2:2769992/1‑150 (MQ=255)
                                            aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGAAAAGGCGGCATCGAAGCCACTCAGCCCCTATCAACCGCCGTTACg                                                                                                          >  2:116576/1‑150 (MQ=255)
                                               gattGCGTTATTATTGTTGCGCTGGAGATAAGCCAAACCCATCGCAACTTTGAGCGCGAGCGGTCCCCCCGCCCCTCTGGGTTGCGGTTTTGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                                                       <  1:814753/147‑1 (MQ=255)
                                                   gCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaaa                                                                                                   >  1:1187434/1‑150 (MQ=255)
                                                     gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                                                 >  2:432882/1‑150 (MQ=255)
                                                       ttAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                                                >  1:2551289/1‑150 (MQ=255)
                                                             ttgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                                          >  1:405989/1‑150 (MQ=255)
                                                             ttgttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATaa                                                                                          >  1:1117447/1‑150 (MQ=255)
                                                               gttgCGCAGTAGATAAGGCAACCCTATCGCAACTTTGAGCGCGTTGGGTCCCCTCGCCCCTCTGGGTAGTGGGTTCGGGTGAGGGGACAAGGCGGCATCGAAGCCAACCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                                                        <  2:1820445/150‑1 (MQ=255)
                                                                   cgcAGGATATAAGGAAACCCGATCGCAACTTTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTCGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                                                    <  2:760774/150‑1 (MQ=38)
                                                                   cgcAGGAGATATGGCAAACCTATCGCAACATTGTGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                                                    <  2:3194668/150‑1 (MQ=39)
                                                                      aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                                                 >  2:125939/1‑150 (MQ=255)
                                                                      aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGGTAGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAGCAAATAACGAGAAGGCCGTCCa                                                                                 >  2:125941/1‑150 (MQ=39)
                                                                           aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                                            >  2:2797353/1‑150 (MQ=255)
                                                                            aaaGCCAACCCGATCGCAACATTGAGCGCGATCGGTCCCCGCGCCCCTCTGGGTAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtc                                                                           <  1:2412911/149‑1 (MQ=39)
                                                                                 caccccGAGCGAAACATTGCGCGCGATGGGTCCCCTCCCCCCTCTGGGGAGGGGGTGTGGGTGCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                                                      <  1:3103732/146‑1 (MQ=255)
                                                                                  acaCCGATCGCAACTTTGCGCGCGATGGGCCCCCTCGCCCCTGTGGGGTTGGGGTTGGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     <  1:233691/148‑1 (MQ=255)
                                                                                  aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  2:3011412/1‑150 (MQ=255)
                                                                                  aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  1:1011761/1‑150 (MQ=255)
                                                                                  aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                                                     >  2:3011480/1‑150 (MQ=255)
                                                                                   aaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                                                    >  1:1110445/1‑150 (MQ=255)
                                                                                   aaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                                                    >  1:1426939/1‑150 (MQ=255)
                                                                                      cGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                                                 >  2:1225312/1‑150 (MQ=255)
                                                                                      cGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                                                 >  2:572371/1‑150 (MQ=255)
                                                                                      cGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                                                 >  2:572426/1‑150 (MQ=255)
                                                                                       taTCGAAACTTTGCGCGCGATCGGCCCCCGCGCCCCTGTGGGGAGTGGGTTTGGGGGTGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                <  1:23323/149‑1 (MQ=255)
                                                                                       gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                >  2:2905659/1‑150 (MQ=255)
                                                                                       gATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                                                >  1:221176/1‑150 (MQ=255)
                                                                                         tCGCAACTTTGAGCGCGATGGGGCCCCTCGCCCCTCTGGGTAGTGGTTTCGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                                              <  2:3066780/150‑1 (MQ=255)
                                                                                          cGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTGTACGAATCAAATAACGATAAGGCCGGCCATCCGTCTCTTTAGCACCCAg                                                             >  2:2371775/1‑150 (MQ=39)
                                                                                            cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                           >  1:730335/1‑150 (MQ=255)
                                                                                            cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                           >  1:2361917/1‑150 (MQ=255)
                                                                                            cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                           >  1:730327/1‑150 (MQ=255)
                                                                                            cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                                           >  1:607285/1‑150 (MQ=255)
                                                                                             aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                                          >  1:3191785/1‑150 (MQ=255)
                                                                                             aaCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                                                          >  1:1594465/1‑150 (MQ=255)
                                                                                                    aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                   >  2:1919847/1‑150 (MQ=255)
                                                                                                    aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                                                   >  2:1919955/1‑150 (MQ=255)
                                                                                                        cgctCGGTCCCCTCGCCCCTGTGGGGTGCGGTTTTGGGTGTGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                               <  1:2480607/147‑1 (MQ=255)
                                                                                                        cgATCGGGCCCCTCGCCCCTCTGGGTAGAGGGTTCGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                               <  2:775047/150‑1 (MQ=255)
                                                                                                        cgATCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGACTCAAATCACGATAAGGCAGTCCATCCTTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                               <  2:3191160/150‑1 (MQ=39)
                                                                                                        cgATCGGCCCCCCCGCCCCTGGGGGGTGCGGGTTTGGGGGTGGGGCCAAGGCGGCATCGAAGCCAATCAGCCCCTATCACCCGCCTTTACGAATCACATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaa                                               <  2:2494945/150‑1 (MQ=255)
                                                                                                         gaccGGTCCCCTCGCCCGTCGGGGGATTGGGTTAGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                              <  1:2104483/147‑1 (MQ=255)
                                                                                                         gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCCGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                              >  1:2650227/1‑150 (MQ=255)
                                                                                                         gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                              >  1:869176/1‑150 (MQ=255)
                                                                                                         gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                                              >  2:1143994/1‑150 (MQ=255)
                                                                                                         gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGGAGGGGGAGGGGAAAGGGCGGCCTCGAAGCCGATCCGCCCCGATCAACCGCCTGTACGAATAAAATAACGATAAGGCAGGCCTTCCGTCTCTGTAGCAGCCAGTTCGGGGTcccgaac                                              >  2:1143987/1‑145 (MQ=25)
                                                                                                          aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                                             >  2:2904538/1‑150 (MQ=255)
                                                                                                              gcccccGCGCCCCTGTGGGGTGTGGGTTCGGGGGAGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTAc                                         <  1:837523/148‑1 (MQ=255)
                                                                                                                 cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                      >  1:2125519/1‑150 (MQ=255)
                                                                                                                 cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaa                                      >  2:3154979/1‑150 (MQ=255)
                                                                                                                 cccTCGCCCCTCGGGGGGGAGGGGTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCCTCCGTCTCTTTAGCAACCAGTTCGTGTTCCAGAAAGGTACaaa                                      >  2:2951784/1‑150 (MQ=255)
                                                                                                                  ccTCCCCCCTCTGGGTAGAGGTTTAGGGTGAGGCCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa                                     <  2:1029787/150‑1 (MQ=255)
                                                                                                                   cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc                                    >  2:2063304/1‑150 (MQ=255)
                                                                                                                           tctGGGGAGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                            <  1:613386/150‑1 (MQ=255)
                                                                                                                           tctGGGGAGAGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                            <  2:1639866/150‑1 (MQ=255)
                                                                                                                           tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                            >  2:1615852/1‑150 (MQ=255)
                                                                                                                           tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAat                            >  2:679818/1‑150 (MQ=255)
                                                                                                                             tGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAatat                          >  2:2963372/1‑150 (MQ=255)
                                                                                                                                ggAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATcgc                       >  2:2113514/1‑150 (MQ=255)
                                                                                                                                    agGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGt                   >  1:2447156/1‑150 (MQ=255)
                                                                                                                                    agGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGt                   >  2:1770394/1‑150 (MQ=255)
                                                                                                                                       gTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGt                >  1:1135846/1‑150 (MQ=255)
                                                                                                                                        tttgggTGCGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCCGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTg               <  1:1928327/147‑1 (MQ=255)
                                                                                                                                           gggggtggggAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGcc            <  2:2706269/144‑1 (MQ=255)
                                                                                                                                           gggTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGcc            <  2:2394397/150‑1 (MQ=255)
                                                                                                                                                gaggggAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtc       <  2:2759509/148‑1 (MQ=255)
                                                                                                                                                   ggACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGAtcgtc    <  1:2324142/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATATCGCGGGTAGTGGCCGGATCGTCGG  >  NC_000913/3608572‑3608869

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: