Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 12.5% 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,6950CG12.5% 607.7 / 25.3 206intergenic (+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 (94/83);  new base G (4/25);  total (98/108)
Fisher's exact test for biased strand distribution p-value = 8.25e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGA  >  NC_000913/3608547‑3608841
                                                                                                                                                    |                                                                                                                                                  
cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATCAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                   >  1:3176823/1‑150 (MQ=255)
cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCACCATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                   >  2:3080683/1‑150 (MQ=255)
cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                   >  1:1840381/1‑150 (MQ=255)
cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                   >  1:215732/1‑150 (MQ=255)
cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                   >  2:377870/1‑150 (MQ=255)
 acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                  <  2:3557148/150‑1 (MQ=255)
 acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACAGTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                  >  2:3078726/1‑150 (MQ=255)
   gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                <  1:3633719/150‑1 (MQ=255)
   gCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAAGGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTggg                                                                                                                                                <  1:3507452/150‑1 (MQ=255)
        tAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGTTCGCTACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGtgag                                                                                                                                           <  1:2462326/150‑1 (MQ=255)
         aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                          >  1:2595003/1‑150 (MQ=255)
         aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                          >  1:3751252/1‑150 (MQ=255)
          gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                         >  2:3937602/1‑150 (MQ=255)
          gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                         >  2:3937320/1‑150 (MQ=255)
          gggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                         >  1:1930763/1‑150 (MQ=255)
           ggATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                        >  2:3030976/1‑150 (MQ=255)
                cccGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                   <  2:339362/149‑1 (MQ=255)
                aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                   >  1:1781776/1‑150 (MQ=255)
                  cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTATGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCTACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                 >  1:3945818/1‑150 (MQ=255)
                  cGGGTTGTGGCTGGCAGTGCTGGCAAATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                 >  1:222667/1‑150 (MQ=255)
                      ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                             <  2:2796657/150‑1 (MQ=255)
                      ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                             <  1:3780571/150‑1 (MQ=255)
                        gtgGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                           >  1:2122594/1‑150 (MQ=255)
                         tggctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                          >  1:3379982/1‑150 (MQ=255)
                            ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGCGCGCAACATTGCGCGCGATCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAAg                                                                                                                       <  2:2401119/150‑1 (MQ=39)
                            ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCATCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                       <  2:1962421/150‑1 (MQ=255)
                            ctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                       >  1:3941146/1‑150 (MQ=255)
                              ggcAGTGGTGGCAGATACGGGGGCGACGGTGCTGGTGACCGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGTGCGCGATCGGTCCCCTCGCCCCTGTGGGTTGAGGGTTAGGGTGAGGGGAAAAggc                                                                                                                     <  2:50447/150‑1 (MQ=39)
                                cAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcgg                                                                                                                   >  1:1190559/1‑150 (MQ=255)
                                   tGCTGGCAGATTCGGGGGCGCGGGTGGTGGTGACAGCTATTGCGTTAATTTTGTTGCGCTGTAGATAAGGCAAACCTATCGAAACATTGAGCGCTATGGGTCCCCTCGCCCCTCTGGGGTGTGGTTTTGGGTGAGGGGAAAAGGCGGCAt                                                                                                                <  2:1277562/150‑1 (MQ=255)
                                   tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCACCATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                >  2:3709288/1‑150 (MQ=255)
                                   tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                >  1:3012408/1‑150 (MQ=255)
                                   tGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGACAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCAt                                                                                                                >  2:2284023/1‑150 (MQ=255)
                                    gctggctgAACCGGGGGCGACGGTGCTGGTGACAGCGAAGGCTTTAAGTTTGTTGCGCAGTATATAAGGCAAACCGTTCGCAACTTTGAGCGAGATCGGTCCCCTCGCCCCTCGGGGGCGCGGGTTCGGGTTAGGGTAAAAGGCGGCATc                                                                                                               <  2:255615/150‑1 (MQ=255)
                                      tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                             >  2:2722765/1‑150 (MQ=255)
                                      tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGa                                                                                                             >  1:1702397/1‑150 (MQ=255)
                                        gCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAg                                                                                                           >  2:2218866/1‑150 (MQ=255)
                                          aGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGcc                                                                                                         >  1:2465751/1‑150 (MQ=255)
                                              aCGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATc                                                                                                     >  1:1636035/1‑150 (MQ=255)
                                                 ggggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCTGGAGATAAGGCAAACCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                  <  2:646052/150‑1 (MQ=255)
                                                  gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATCAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                 >  2:1162585/1‑150 (MQ=255)
                                                  gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcc                                                                                                 >  2:3221542/1‑150 (MQ=255)
                                                      gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTa                                                                                             >  1:2797515/1‑150 (MQ=255)
                                                                    cAGCTAATGCGTTAAGATTGTTGCGCTGGATATAAGGCAAACCGATCGCAACATTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                               <  2:59195/150‑1 (MQ=255)
                                                                    cAGCGAATGCGTTAAGATGGGTGCGCGGGTGATACGGCAAACCGATGGCAACCTTGAGCGCGAGCGGTCCCCTCGCCCCTGTGGTGTGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                               <  2:3740706/150‑1 (MQ=37)
                                                                           tGCGTTATTTTTGTTGCGCAGTTTATAAGGCAACCCGATCTAAACTTTGAGCGCGAGCGGTCCCCTCGCCCCTCTGGGGTGTGGTTTTGGGTGTGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaa                                                                        <  1:1382145/150‑1 (MQ=255)
                                                                                  aGATTGTTGCGCAGGAGATAAGGCAAACCGATCTCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                 >  1:859778/1‑150 (MQ=255)
                                                                                  aGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                 >  1:694047/1‑150 (MQ=255)
                                                                                    ttttttgCGCTGGAGATAAGCAAAACCTATCGCAAATTTGAGCGCTATCGGCCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                               <  2:3868652/146‑1 (MQ=255)
                                                                                    tttgttgCGCTGGAGTTAAGGAAAACCGATCGCAACTTTGCGCGCTAGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGTGCGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                               <  2:237986/149‑1 (MQ=38)
                                                                                    aTTGTTGCGCTGGTGATAAGGCAAACCGATCGCAACTTTGCGCGCGTTGGGGCCCCTCGCCCCTCGGGGGAGAGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATa                                                               <  2:1714382/150‑1 (MQ=255)
                                                                                       gttgCGCAGGATATAAGGCAAACCTATCGCAAATTTGAGCGCTATCGGTCCCCTCGCCCCTCTGGGTATTGGTTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                            <  2:2465751/150‑1 (MQ=255)
                                                                                       gttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGACCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                            >  1:3260380/1‑150 (MQ=255)
                                                                                              aGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCa                                                     >  1:3870156/1‑150 (MQ=255)
                                                                                                gagaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                   >  1:2090300/1‑150 (MQ=255)
                                                                                                gagaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATc                                                   >  2:148836/1‑150 (MQ=255)
                                                                                                gagaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTGTACGAATCAAATAACGATAAGGCAGTCCATc                                                   >  1:3859238/1‑150 (MQ=255)
                                                                                                   atatGGCAACCCGATGGAAACATTGTGCGCGATCGTTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                <  2:2412617/150‑1 (MQ=255)
                                                                                                   aTAAGGCAAACCGATCGAAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGt                                                <  2:3733764/150‑1 (MQ=255)
                                                                                                       ggCAAACCGATCGCAACATTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGGGGGTTCGGGTGCGGTCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                            <  2:1872887/150‑1 (MQ=255)
                                                                                                       ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                            >  1:2787132/1‑150 (MQ=255)
                                                                                                       ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                            >  2:2849608/1‑150 (MQ=255)
                                                                                                       ggCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTctct                                            >  2:2597285/1‑150 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                           >  2:1189379/1‑150 (MQ=255)
                                                                                                        gCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                           >  1:1394868/1‑150 (MQ=255)
                                                                                                         cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTCACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                          >  1:2925674/1‑150 (MQ=255)
                                                                                                         aaacccGCTCGCAAATTGGAGCGCGACCGGTCCCCTCGCCCCTCGGGGGTGTGGTTTTGGTTGTGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCttt                                          <  2:159516/146‑1 (MQ=255)
                                                                                                          aaaCCTATCGCAACTTTGCGCGCGAGCGGACCCCGCGCCCCTCTGGGGTGAGGGTTGGGGTGAGGGTACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                         <  1:2039694/150‑1 (MQ=255)
                                                                                                          aaaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                         >  1:2720043/1‑150 (MQ=255)
                                                                                                          aaaCCGACCTCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                         <  1:2986459/150‑1 (MQ=255)
                                                                                                               gctCGAAAAATTGAGCGCGATCGGTCCCCCCGCCCGGCTGGGGAGTGGGTTGGGGGGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                    <  2:3988863/148‑1 (MQ=255)
                                                                                                               gATCGCAACATTGCGCGCGATCGGTCCCCTCGCCCCTCTGGGGTTTGGTTTTGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                    <  2:908234/150‑1 (MQ=255)
                                                                                                                   tcAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGTGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGt                                <  2:762120/149‑1 (MQ=255)
                                                                                                                    cAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGtt                               >  1:4067122/1‑150 (MQ=255)
                                                                                                                     aaCTTTGAGCGCGAGGGGGCCCCCCGCCCCTGTGGGGTGTGGTTTTGGGGGCGGGCACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                              <  1:1917717/150‑1 (MQ=255)
                                                                                                                     aaCATTGCGCGCTATCGGTCCCCGCGCCCCTGTGGGGAGTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTc                              <  1:434657/150‑1 (MQ=255)
                                                                                                                       ctttGCGCGCGATCGGACCCCTCGCCCCTCTGGGGAGAGGTTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                            <  1:3446828/148‑1 (MQ=39)
                                                                                                                       cATTGAGCGCGATGGGTCCCCGCGCCCCTGGGGTGAGTGGTTTTGGGTGTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgt                            <  1:125606/150‑1 (MQ=255)
                                                                                                                       cATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCTCCTTGACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTGTAGCAACCAGTTCgt                            >  2:773071/1‑150 (MQ=255)
                                                                                                                           gcgcgcGATCGGCCCCCGCGCCCCTCTGGGGAGTGGTTTTGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                        <  1:514274/148‑1 (MQ=255)
                                                                                                                           gAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTc                        >  2:1456520/1‑150 (MQ=255)
                                                                                                                            cgcgcTATGGGCCCCCCCGCCCCGGGGGGTTGTGGGTTCGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                       <  2:1809390/149‑1 (MQ=255)
                                                                                                                            aGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTcc                       >  1:2019527/1‑150 (MQ=255)
                                                                                                                              cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                     >  1:1073272/1‑150 (MQ=255)
                                                                                                                              cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                     >  1:728281/1‑150 (MQ=255)
                                                                                                                              cgcgATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCAt                     >  2:2164771/1‑150 (MQ=255)
                                                                                                                                  ttgggcccccTCGCCCGTGTGGGTTTTGGTTTTGGGTGTGGGCAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                 <  2:1926092/144‑1 (MQ=255)
                                                                                                                                  aTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                 <  1:3058108/150‑1 (MQ=255)
                                                                                                                                  aTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAg                 >  2:3339294/1‑150 (MQ=255)
                                                                                                                                    ggttCCCGTCGCCCCTCGGGGTTTCGGGTTTGGTTTTGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGACTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt               <  2:462491/147‑1 (MQ=255)
                                                                                                                                    cGGCCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGACAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGt               <  2:1240499/150‑1 (MQ=255)
                                                                                                                                     ctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTa              >  1:2151423/5‑150 (MQ=255)
                                                                                                                                          ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCCTCAAGGTACaaca         >  1:3537709/1‑148 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCTACCAGTTCGTGTTCCATAAAGGTACAAAAc        >  1:2544566/1‑150 (MQ=255)
                                                                                                                                            tCGCCCCTGTGGGGAGTGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc       <  2:208855/150‑1 (MQ=255)
                                                                                                                                            ccGCCCCTCGGGGGAGTGGGTTTGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAcc       <  2:1633489/149‑1 (MQ=255)
                                                                                                                                            acGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATACGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAACAcc       <  1:643479/149‑1 (MQ=255)
                                                                                                                                             cGCCCCTGTGGGTAGAGGTTTTGGTTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAAccc      <  2:1298064/150‑1 (MQ=255)
                                                                                                                                              gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa     >  1:571907/1‑150 (MQ=255)
                                                                                                                                                 ccTGTGGGGTGTGGGTTTGGGTGAGGGTAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGa  <  1:2218866/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGA  >  NC_000913/3608547‑3608841

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: