Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 10.9% 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,6950CG10.9% 1375.0 / 45.1 437intergenic (+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 (243/138);  new base G (3/49);  total (246/191)
Fisher's exact test for biased strand distribution p-value = 1.36e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 450 total aligned reads displayed.

CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  NC_000913/3608546‑3608844
                                                                                                                                                     |                                                                                                                                                     
 cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                      >  1:1015876/1‑150 (MQ=255)
        ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               >  1:13488/1‑150 (MQ=255)
         tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  1:5444823/150‑1 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                             <  1:3096891/150‑1 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             <  2:6688428/150‑1 (MQ=255)
          aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                             <  1:2563931/150‑1 (MQ=255)
           gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                            <  1:2742701/150‑1 (MQ=255)
             gATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCTCCCCAGAGGGGCGAgggga                                                                                                                                          <  2:5159034/150‑3 (MQ=255)
              aTGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTa                                                                                                                                         >  2:1103456/1‑150 (MQ=255)
               tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                        <  2:830974/150‑1 (MQ=255)
               tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                        <  1:6310621/150‑1 (MQ=255)
                   cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                    >  1:6814261/1‑150 (MQ=255)
                    gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   <  2:139514/150‑1 (MQ=255)
                      ggtgtgGCTGGCAGTGCTGGCAGAGACGGGGGCGCCGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 <  2:4598305/148‑1 (MQ=255)
                      gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 <  2:5529207/150‑1 (MQ=255)
                             ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                          >  1:6033147/1‑150 (MQ=255)
                              tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGATAAgg                                                                                                                         <  2:3752415/150‑1 (MQ=255)
                              tggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAgg                                                                                                                         >  2:6511611/1‑150 (MQ=255)
                                gcAGTGGTGGCATATCGGGGGGCGAGGGTGCTGGTGACAGGGTTTGCGTTAATATTGTTGCGCAGTATATATGGCAAACCGATCGAAACTTTGAGCGCGATGGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGACAAggcg                                                                                                                       <  2:1882300/150‑1 (MQ=25)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                >  2:5816620/1‑150 (MQ=255)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                >  1:379616/1‑150 (MQ=255)
                                       tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                                >  2:6120444/1‑150 (MQ=255)
                                        ggCAGATACGGGGGCGCGGGTGCTGGTTACAGCTATGGGGTTAAGATTGTGGCGCAGTATAAAAGGAAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGTGTGGGTTAGGGTGAGGGTAAAAGGCGGCATCGaa                                                                                                               <  2:7006839/150‑1 (MQ=255)
                                          cAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGc                                                                                                             >  1:1637300/1‑150 (MQ=255)
                                             aTACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCaa                                                                                                          >  2:2847859/1‑150 (MQ=255)
                                               aCGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATc                                                                                                        >  1:1765589/1‑150 (MQ=255)
                                                   gggCGAGGGTGCTGGTGACTGCTATTGGTATATGATTTGTGCGCTTTAGATACGCCAAACCGATGGAAAATTTGCGCGCGATGGGGCCCCGCGCCCCTGGGGGGAGTGGTTTGGGGGGCGGGCAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                    <  2:4306953/150‑1 (MQ=255)
                                                   gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTCGGGTGAGGGGAAAAGGCGGCCTCGAAGCCAATCAGcc                                                                                                    >  2:6442986/1‑150 (MQ=255)
                                                   gggCGACGGTGCAGGTGACAGCGAATGCGTTAAGATTGTAGCGCAGGAGATAAGGCAAACCGATCCCAACAAAGAGCGCGATCCGTCCCCACCCCCCTCAGGGGACAGGGATAGCGAGACGCCAAAAGCCCCCATCGAACCCAAACAGcc                                                                                                    >  2:3904759/1‑150 (MQ=255)
                                                     gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                                                                  >  1:3406673/1‑150 (MQ=255)
                                                           gTGCTGGTGACAGCTATTGGTTTATTTTTGTGGCGCTGTTTATACGGCAAACCTATCGAAACATTGCGGGCTAGGGGTCCCCCCGCCCCTCTGGGGTGTGGTTTCGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCaa                                                                                            <  2:2113894/150‑1 (MQ=255)
                                                            tGCTGGTTACAGCGAATGCTTTATGATTGTTGCGCAGGTGATATGGCAAACCTAGCGCAACATTGCGCGCTATCGGTCCCACCGCCCCTCTGGGGAGTGGTTTAGGTTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                           <  2:5837525/150‑1 (MQ=255)
                                                            tGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAc                                                                                           >  1:3945871/1‑150 (MQ=255)
                                                             gCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                          >  2:333530/1‑150 (MQ=255)
                                                             gCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAAcc                                                                                          >  1:6117781/1‑150 (MQ=255)
                                                               tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                        >  2:4175253/1‑150 (MQ=255)
                                                               tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                        >  1:867178/1‑150 (MQ=255)
                                                               tGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGc                                                                                        >  1:1976443/1‑150 (MQ=255)
                                                                 gTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                      >  1:5062189/1‑150 (MQ=255)
                                                                    aCAGCTAATGCTTTAATTTTGTTGCGCTGGAGATAAGGAAAACCTATCGAAACATTGAGCGCTATGGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                   <  2:2992064/150‑1 (MQ=37)
                                                                    aCAGCTAATGCTTTAAGATTGTTGCGCTGTAGATACGGCAAACCTATCGCAACATTGTGCGCTATCGGTCCCGTCGCCCCTCTGGGGTGAGGGTTTGGTTGAGGGGAAAAGTCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                   <  2:7136194/150‑1 (MQ=255)
                                                                     cAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                  >  1:4481000/1‑150 (MQ=255)
                                                                       gCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGACAAGGCGGCATCGAAGCCACTCAGCCCCTATCAACCGCCTTTACGa                                                                                >  1:4855177/1‑150 (MQ=255)
                                                                         gattGCTTTAAGATTGTGGCGCAGGAGATATGCCAAACCGAGGGCAACATTGAGCGCTCTGGGCCCCCGCGCCCCTGTGGGTTGTGGGTTCGGGTGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                              <  1:4309293/147‑1 (MQ=255)
                                                                         gAATGCGTTAAGATTGTTGCTCAGGAGATCAGGCAAACCGATCGCAACATTGAGCGCGTTCGGTCCCCTCGCCCCTCTGGGTTGAGGGTTATGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAAt                                                                              <  2:5383573/150‑1 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  1:6179391/1‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  2:727841/1‑150 (MQ=255)
                                                                               gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATa                                                                        >  2:2395161/1‑150 (MQ=255)
                                                                                ttatgtttgttgCGCAGTATATAAGCCAAACCGATCTCAACTTTGAGCGCGAGCGGTCCCCTCGCCCCTGTGGGTTGTGGGTTTGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATaa                                                                       <  1:5774263/144‑1 (MQ=255)
                                                                                   aGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                    >  1:1727538/1‑150 (MQ=255)
                                                                                     aTTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                                  >  1:74767/1‑150 (MQ=255)
                                                                                           gcgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAg                                                            >  1:119054/1‑150 (MQ=255)
                                                                                           gcgcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAg                                                            >  1:5704439/1‑150 (MQ=255)
                                                                                             gctgtagaTAAGGAAAACCGTTCGCAACTTTGTGCGCGTTCGGTCCCCTCGCCCCTGTGGGGTGTGGTTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTc                                                          <  2:168863/145‑1 (MQ=255)
                                                                                               aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                        >  1:3539737/1‑150 (MQ=255)
                                                                                               aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                        >  1:5569629/1‑150 (MQ=255)
                                                                                                 gagaTAAGGCAAACCGATCGCACCTTTGCGCGCGATCGGTCCCCCCGCCCCTGTGGGGTGTGGTTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                      <  1:7205780/150‑1 (MQ=255)
                                                                                                   gaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTGTACGAATCAAATAACGATACGGCAGTCCATCCg                                                    >  2:1293490/1‑150 (MQ=255)
                                                                                                      accccAAACCGATCCAAAAATGGACCGCGAGCGGACCCCCCCCCCCACTGGGCAGAGGTTTAGGCTGAGGCCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGACTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                                 <  2:2573707/146‑1 (MQ=255)
                                                                                                         gcaccccGCTCGCACCATGGCGCGCGACCGGCCCCCTCGCCCCTTTGGGGTGCGGGTTAGGGTGCGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              <  2:4970385/145‑1 (MQ=255)
                                                                                                         gaaaaCCTACCGCAACTTTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGTTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              <  2:592993/148‑1 (MQ=255)
                                                                                                         gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  1:4036402/1‑150 (MQ=255)
                                                                                                         gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              >  2:1160756/1‑150 (MQ=255)
                                                                                                         gCAAACCGATCGCAAAATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGTTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                              <  2:1784775/150‑1 (MQ=255)
                                                                                                          cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                             >  1:4747881/1‑150 (MQ=255)
                                                                                                           caaCCGATCGCAACATTGCGCGCGATCGGGCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGGGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                            <  2:493820/149‑1 (MQ=255)
                                                                                                            aaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGATGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAg                                           <  2:3612386/150‑1 (MQ=255)
                                                                                                             aCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGc                                          >  2:3446256/1‑150 (MQ=255)
                                                                                                               cGATCTAAACTTTGAGCGGCATCGCGCCCCGCGCCCGTGTGGGGTGTGGGTTTGTGCGCGGTCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCaa                                        <  2:5736969/150‑1 (MQ=255)
                                                                                                                gATCGCAACTTTGCGCGGGATGGGTCCCCTCGCCCCTGGGGGTTTTGGTTTCGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                       <  2:2656194/150‑1 (MQ=255)
                                                                                                                 agcgaaaCTTTGCGCGCGCTCGGTCCCCTCCCCCCTGGGGGTTGGGGGTTCGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                      <  1:3326893/145‑1 (MQ=255)
                                                                                                                 aTCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAcc                                      >  2:1762528/1‑150 (MQ=255)
                                                                                                                  tCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCa                                     >  1:2172902/1‑150 (MQ=255)
                                                                                                                    gCAACATTGCGCGCGATCGGTCCCCGCGCCCCTCTGGGGAGGGGGTTAGGGTGAGGGGAAAAGGCGGCAGCGAAGCCAATCAGCCCCTGTCCACCGCCTTTACGACTCAAATAACGAGAAGGCAGTCCATCCGTCTCTTTCGCTGCCAGt                                   >  2:4471921/1‑150 (MQ=37)
                                                                                                                     cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  >  1:4790881/1‑150 (MQ=255)
                                                                                                                     cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                                  >  1:5365711/1‑150 (MQ=255)
                                                                                                                       acttggcgcgcgcgcGGGCCCCTCGCCCCTGTGGGTTGTGGGTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCg                                <  2:4871746/143‑1 (MQ=255)
                                                                                                                             cgcgcGAGCGGCCCCCGCGCCCCTGTGGGTTGTGGGTGCGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                          <  2:6249683/149‑1 (MQ=255)
                                                                                                                              gcgcGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa                         >  1:3691312/1‑150 (MQ=255)
                                                                                                                              gcgcGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCa                         >  2:6295805/1‑150 (MQ=255)
                                                                                                                               cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                        >  2:4708463/1‑150 (MQ=255)
                                                                                                                               cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                        >  1:1795919/1‑150 (MQ=255)
                                                                                                                                  gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATaaa                     >  2:2281982/1‑150 (MQ=255)
                                                                                                                                   aTCGGTCCCCTCGCCCCTCTGGGTAGAGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    <  1:4136877/150‑1 (MQ=255)
                                                                                                                                   aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    <  1:6006478/150‑1 (MQ=255)
                                                                                                                                   aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                    >  1:1454473/1‑150 (MQ=255)
                                                                                                                                    tCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                   >  1:701007/1‑150 (MQ=255)
                                                                                                                                     cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt                  >  1:2635215/6‑150 (MQ=255)
                                                                                                                                     cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTGCGTGTTCCATAAAGGt                  >  2:3989467/6‑150 (MQ=255)
                                                                                                                                        tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa               >  1:5151636/1‑150 (MQ=255)
                                                                                                                                        tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACa               >  1:4923595/1‑150 (MQ=255)
                                                                                                                                        tCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaa              >  1:245817/1‑150 (MQ=255)
                                                                                                                                           ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCACCCAGTTCGTGTTCCATAAAGGTACaaaa            >  1:6558208/1‑150 (MQ=255)
                                                                                                                                            cccGCCCCTGTGGGGTGTGGGTTCGGGTGGGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCGTTTACGAATCAATTAACGATAAGGCAGTCCTTCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAc           <  1:6707028/148‑1 (MQ=255)
                                                                                                                                             tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc          >  1:4392068/1‑150 (MQ=255)
                                                                                                                                               gCCCCTCTGGGGTGTGGGTTAGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  2:5659333/150‑1 (MQ=255)
                                                                                                                                               gCCCCTCTGGGGAGAGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa        <  2:5999958/150‑1 (MQ=255)
                                                                                                                                                ccccTGGGGGGAGTGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg       <  1:524197/150‑1 (MQ=255)
                                                                                                                                                ccccTCTGGGTTTAGGTTTTTGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg       <  2:531070/150‑1 (MQ=255)
                                                                                                                                                  ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa     >  1:6891810/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAATA  >  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: