Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,193,141 G→A G211G (GGC→GGT mnmA ← tRNA‑specific 2‑thiouridylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,193,1410GA98.5% 185.3 / NA 65G211G (GGC→GGTmnmAtRNA‑specific 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (37/27);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

GCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  NC_000913/1192903‑1193387
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
gctgGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTc                                                                                                                                                                                                                                              <  1:629568/250‑1 (MQ=255)
        aTCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAAttt                                                                                                                                                                                                                                     <  1:35559/251‑1 (MQ=255)
           aaCCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGGGCTCACCAATTTCATCGCCATCGACGGGAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTAGCGGAATTTACg                                                                                                                                                                                                                                  >  2:232681/1‑251 (MQ=255)
            ccccGACAGACCTCAGCCGCGGGGGTTCATGGCCCTGAGCGAAACCCAGAATGTGGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGTTACCCTGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATGTTGCCCGGGTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGt                                                                                                                                                                                                                                 <  2:70766/250‑1 (MQ=255)
                     acaTAAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGGTGTCGACGTCTGTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTGGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATg                                                                                                                                                                                                                        <  2:253390/251‑1 (MQ=255)
                               cgGGTGTGCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATAACACTGTTCTTCGGTACCTGCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACAGCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTGGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAAt                                                                                                                                                                                                                <  1:286836/249‑1 (MQ=255)
                                 ggTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCg                                                                                                                                                                                                            <  1:382171/251‑1 (MQ=255)
                                   tgctcaaggCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGCGATAACGACCCAGGAACTCGCGGAATTTACGTTTGCCACTGAAGCAAATGCCGGt                                                                                                                                                                                                          >  2:415716/8‑251 (MQ=255)
                                         tggcccagcgcgccacacatactgttgtTTTCAACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCAACGATACCCAGACCTTTACGCTGACCGCGAGTGTGCTACATAAGCCCCTGGGGATCACCAATTTCATCGCCCTCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCGGGAACTCGCGGAATTTACGTTTGCCAATGAAGCAAATGCAGGGAGAGTc                                                                                                                                                                                                    >  2:467987/23‑251 (MQ=255)
                                          ggCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTc                                                                                                                                                                                                    <  1:17674/250‑1 (MQ=255)
                                                   cGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGCTTTTGCCCGGTTGCGCCGGGCGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTttcttc                                                                                                                                                                                          >  2:301142/1‑251 (MQ=255)
                                                   cGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGAAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAATCAAATGCCGGTAGAGTCTTTTttcttt                                                                                                                                                                                          >  2:98802/1‑250 (MQ=255)
                                                      cAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTcgcg                                                                                                                                                                                       <  1:301142/251‑1 (MQ=255)
                                                      cAACCAGAATGTTGTTTTCGACGTCTTTGCCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCCCCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCACATGCCGGTAGAGTCTTTTTTCTTcgcg                                                                                                                                                                                       >  1:582533/1‑251 (MQ=255)
                                                                gttgttTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTgaccaga                                                                                                                                                                              <  1:165780/250‑1 (MQ=255)
                                                                         gACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTTCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCtt                                                                                                                                                                    >  1:639891/1‑251 (MQ=255)
                                                                          aCGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTc                                                                                                                                                                   >  1:453964/1‑251 (MQ=255)
                                                                           cGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGCTACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGagctac                                                                                                                                                                   >  2:386389/1‑245 (MQ=255)
                                                                           cGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTc                                                                                                                                                                   >  1:59891/1‑250 (MQ=255)
                                                                                 tGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATc                                                                                                                                                            >  1:371178/1‑251 (MQ=255)
                                                                                             tccccGGTTCTTCGGTACCTCCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGAGACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTg                                                                                                                                                 <  2:321162/246‑1 (MQ=255)
                                                                                                 cGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGtt                                                                                                                                            >  2:405003/1‑251 (MQ=255)
                                                                                                  ggTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTACGGttt                                                                                                                                           >  2:253799/1‑251 (MQ=255)
                                                                                                  ggTTCTTCGGTACCTTCTTTGATGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCAACGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGGAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGttt                                                                                                                                           >  1:189951/1‑251 (MQ=255)
                                                                                                          ggTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGtt                                                                                                                                   >  2:566957/1‑251 (MQ=255)
                                                                                                               cttcttTGGAGCCACCGATACCCAGACCTTTACGCTGCCCGAGAGTGTGATACATCAACCCCTGGCGCTCACCAATTTCATCGCCATCGCCGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAAAGACCCAGGAACTCGCGGAATTTACGTTAGCCAATGAAGCAACTGCCGGACGAGTCTTTTTTCTTCGCTGTGACCCGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCg                                                                                                                              >  2:37885/1‑251 (MQ=255)
                                                                                                                       gTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGCGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGa                                                                                                                       >  1:65643/1‑250 (MQ=255)
                                                                                                                            aCCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGc                                                                                                                >  2:338530/1‑251 (MQ=255)
                                                                                                                            aCCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGc                                                                                                                >  1:484542/1‑251 (MQ=255)
                                                                                                                              cGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCtt                                                                                                              <  1:676477/251‑1 (MQ=255)
                                                                                                                               gctACCCAGACCTTTACGCTGACCGTGATGGTGATTCAGCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGTAATTTACGGTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCGTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGTGAACAGGCttt                                                                                                             <  2:59891/249‑1 (MQ=255)
                                                                                                                                 tACCCAGACCTTTACGCGGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCAACTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCATTgc                                                                                                           >  2:428045/1‑251 (MQ=255)
                                                                                                                                  aCCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTgcc                                                                                                          >  2:160992/1‑250 (MQ=255)
                                                                                                                                               cGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGCCGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCCATCTGCTCATg                                                                                             >  2:353773/1‑251 (MQ=255)
                                                                                                                                                gCTGACCGATAGTGTGATACATCAGCCACTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATgg                                                                                            <  1:548206/251‑1 (MQ=255)
                                                                                                                                                gCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATgg                                                                                            <  1:39065/251‑1 (MQ=255)
                                                                                                                                                  tGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCt                                                                                          >  1:514810/1‑251 (MQ=255)
                                                                                                                                                   gACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTGGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCt                                                                                          <  1:415716/250‑1 (MQ=255)
                                                                                                                                                       gagaGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCACTCTGCTCATGGATGAGCg                                                                                     >  1:643450/1‑251 (MQ=255)
                                                                                                                                                               gATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAg                                                                              <  1:438372/250‑1 (MQ=255)
                                                                                                                                                                    aTCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAACAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCCGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGCGCGTATTACGGAAGt                                                                         >  2:6839/1‑250 (MQ=255)
                                                                                                                                                                         ccccTGGAGCTCACCAATTTCATCGACATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTAAGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTAGCGGTGACCAGCCCAAGAACTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCTCAATCTGCTCAGGGCTGCGCGTTTAAAGGAAGAAGctcc                                                                   >  2:471256/1‑250 (MQ=255)
                                                                                                                                                                                   tCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTAtt                                                         >  1:119686/1‑251 (MQ=255)
                                                                                                                                                                                    cACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTc                                                               >  1:440451/1‑244 (MQ=255)
                                                                                                                                                                                    cACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTg                                                        <  1:353773/251‑1 (MQ=255)
                                                                                                                                                                                    cACCAAGTTCATCGCCATCGACGTTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTGCGCCGACCGTGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCTTATAAAGGAAGTAGCTCTGGTc                                                               <  2:440451/244‑1 (MQ=255)
                                                                                                                                                                                     aCCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGc                                                       >  1:579369/1‑251 (MQ=255)
                                                                                                                                                                                        aaTTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGt                                                    >  2:164057/1‑251 (MQ=255)
                                                                                                                                                                                             cATCGCCATCGACTGGAATGAGTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCGTCGCGGTGACCAGACCAAGATCTTCAGCCATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCt                                                                                                    <  2:191745/198‑1 (MQ=255)
                                                                                                                                                                                             cATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCt                                                                                                    >  1:191745/1‑198 (MQ=255)
                                                                                                                                                                                             aaTCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCGCTGGTCTTTAGTGCTGTCCAGa                                               <  2:579369/250‑1 (MQ=255)
                                                                                                                                                                                                    aTCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGCAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgc                                         >  1:238908/1‑250 (MQ=255)
                                                                                                                                                                                                    aTCGACGGTAATGATTTTGCCAGGTTGCGCAGGGAGCTAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGATGCAAATGCAGGTAGCGTCTTTTTTCTTCTCGGGGACCAGACCAAGATATTCAGCAATCTTCCGCACCTGCGGTTTTTCCAGTTCGCCGACCGGAACAAGGCTTTGCGCAATCTGCTCATGGATGCGCATCTAAAGACAGTATCTATGGTCTTTCTTGCTGTGCAGACCACga                                         >  2:530356/1‑249 (MQ=255)
                                                                                                                                                                                                     acGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcag                                       <  1:317742/250‑1 (MQ=255)
                                                                                                                                                                                                      cGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcag                                       <  1:405003/250‑1 (MQ=255)
                                                                                                                                                                                                            gaaTGATTTTGCCCGGGGGCGCCGGGAGACATCGGCTAAGGAACGCGCGGAATTTTCGGTCGCCATTGAAGCAAATGCCGGTAGAGTTTTTTTTCCTCGCGGTGGCCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTGCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCGGAGCGTAGAAAGGAAGTAGCTCTGGTATTTATTGCTTTCCAGCCCACGCAGCAggcgg                                <  2:174732/250‑1 (MQ=255)
                                                                                                                                                                                                             aaTGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAg                                                >  1:326232/1‑234 (MQ=255)
                                                                                                                                                                                                             aaTGATTTGGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTGCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTGTTTCCAGTGCGCCGACCGGGAACAGGCTTTGCGCAATCTGCGCATGGCTGAGCGTATAAAGGAAGTAGCGCTGGGCTTTATTGCTGTCCAg                                                <  2:326232/234‑1 (MQ=255)
                                                                                                                                                                                                                    ttGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAAGCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGc                         <  1:645570/250‑1 (MQ=255)
                                                                                                                                                                                                                         cgggtGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGa                   <  2:371178/247‑1 (MQ=255)
                                                                                                                                                                                                                          ggTTGCGCCGGGCGATAACGACCCAGTAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGGAGCTCTGGTCTTTATTGCTGTCCAACCCACGCAGCAGGCGGCTCTTGccctcacc                  >  1:596795/1‑245 (MQ=255)
                                                                                                                                                                                                                          ggTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATc                                                                                                     >  1:283013/1‑168 (MQ=255)
                                                                                                                                                                                                                          ggTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATc                                                                                                     <  2:283013/168‑1 (MQ=255)
                                                                                                                                                                                                                               cgcCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCGTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCgg             <  1:566429/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aaCGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAACGGAAGTCGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCCTCGACATCGGCCCgctggacg  >  2:312816/1‑244 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
GCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  NC_000913/1192903‑1193387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: