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→GGT) mnmAtRNA‑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

GCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  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)
                                                                                                                            aCCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGc                                                                                                                >  2:338530/1‑251 (MQ=255)
                                                                                                                            aCCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGc                                                                                                                >  1:484542/1‑251 (MQ=255)
                                                                                                                              cGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCtt                                                                                                              <  1:676477/251‑1 (MQ=255)
                                                                                                                               gctACCCAGACCTTTACGCTGACCGTGATGGTGATTCAGCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGTAATTTACGGTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCGTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGTGAACA‑GGCttt                                                                                                             <  2:59891/249‑1 (MQ=255)
                                                                                                                                 tACCCAGACCTTTACGCGGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCAACTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCATTgc                                                                                                           >  2:428045/1‑251 (MQ=255)
                                                                                                                                  aCCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTgcc                                                                                                          >  2:160992/1‑250 (MQ=255)
                                                                                                                                               cGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGCCGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCCATCTGCTCATg                                                                                             >  2:353773/1‑251 (MQ=255)
                                                                                                                                                gCTGACCGATAGTGTGATACATCAGCCACTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATgg                                                                                            <  1:548206/251‑1 (MQ=255)
                                                                                                                                                gCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATgg                                                                                            <  1:39065/251‑1 (MQ=255)
                                                                                                                                                  tGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCt                                                                                          >  1:514810/1‑251 (MQ=255)
                                                                                                                                                   gACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTGGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCt                                                                                          <  1:415716/250‑1 (MQ=255)
                                                                                                                                                       gagaGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCACTCTGCTCATGGATGAGCg                                                                                     >  1:643450/1‑251 (MQ=255)
                                                                                                                                                               gATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAg                                                                              <  1:438372/250‑1 (MQ=255)
                                                                                                                                                                    aTCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAACAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCCGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGCGCGTATTACGGAAGt                                                                         >  2:6839/1‑250 (MQ=255)
                                                                                                                                                                         ccccTGGAGCTCACCAATTTCATCGACATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTAAGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTAGCGGTGACCAGCCCAAGAACTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCTCAATCTGCTCAGGGCTGCGCGTTTAAAGGAAGAAGctcc                                                                   >  2:471256/1‑250 (MQ=255)
                                                                                                                                                                                   tCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTAtt                                                         >  1:119686/1‑251 (MQ=255)
                                                                                                                                                                                    cACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTc                                                               >  1:440451/1‑244 (MQ=255)
                                                                                                                                                                                    cACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTg                                                        <  1:353773/251‑1 (MQ=255)
                                                                                                                                                                                    cACCAAGTTCATCGCCATCGACGTTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTGCGCCGACCGTGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCTTATAAAGGAAGTAGCTCTGGTc                                                               <  2:440451/244‑1 (MQ=255)
                                                                                                                                                                                     aCCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGc                                                       >  1:579369/1‑251 (MQ=255)
                                                                                                                                                                                        aaTTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGt                                                    >  2:164057/1‑251 (MQ=255)
                                                                                                                                                                                             cATCGCCATCGACTGGAATGAGTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCGTCGCGGTGACCAGACCAAGATCTTCAGCCATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCt                                                                                                    <  2:191745/198‑1 (MQ=255)
                                                                                                                                                                                             cATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCt                                                                                                    >  1:191745/1‑198 (MQ=255)
                                                                                                                                                                                             aaTCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCGCTGGTCTTTAGTGCTGTCCAGa                                               <  2:579369/250‑1 (MQ=255)
                                                                                                                                                                                                    aTCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGCAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgc                                         >  1:238908/1‑250 (MQ=255)
                                                                                                                                                                                                    aTCGACGGTAATGATTTTGCCAGGTTGCGCAGGGAGCTAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGATGCAAATGCAGGTAGCGTCTTTTTTCTTCTCGGGGACCAGACCAAGATATTCAGCAATCTTCCGCACCTGCGGTTTTTCCAGTTCGCCGACCGG‑AACAAGGCTTTGCGCAATCTGCTCATGGATGCGCATCTAAAGACAGTATCTATGGTCTTTCTTGCTGTGCAGACCACga                                         >  2:530356/1‑249 (MQ=255)
                                                                                                                                                                                                     acGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcag                                       <  1:317742/250‑1 (MQ=255)
                                                                                                                                                                                                      cGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcag                                       <  1:405003/250‑1 (MQ=255)
                                                                                                                                                                                                            gaaTGATTTTGCCCGGGGGCGCCGGGAGACATCGGCTAAGGAACGCGCGGAATTTTCGGTCGCCATTGAAGCAAATGCCGGTAGAGTTTTTTTTCCTCGCGGTGGCCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTGCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCGGAGCGTAGAAAGGAAGTAGCTCTGGTATTTATTGCTTTCCAGCCCACGCAGCAggcgg                                <  2:174732/250‑1 (MQ=255)
                                                                                                                                                                                                             aaTGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAg                                                >  1:326232/1‑234 (MQ=255)
                                                                                                                                                                                                             aaTGATTTGGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTGCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTGTTTCCAGTGCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCGCATGGCTGAGCGTATAAAGGAAGTAGCGCTGGGCTTTATTGCTGTCCAg                                                <  2:326232/234‑1 (MQ=255)
                                                                                                                                                                                                                    ttGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAAGCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGc                         <  1:645570/250‑1 (MQ=255)
                                                                                                                                                                                                                         cgggtGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGa                   <  2:371178/247‑1 (MQ=255)
                                                                                                                                                                                                                          ggTTGCGCCGGGCGATAACGACCCAGTAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGGAGCTCTGGTCTTTATTGCTGTCCAACCCACGCAGCAGGCGGCTCTTGccctcacc                  >  1:596795/1‑245 (MQ=255)
                                                                                                                                                                                                                          ggTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATc                                                                                                     >  1:283013/1‑168 (MQ=255)
                                                                                                                                                                                                                          ggTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATc                                                                                                     <  2:283013/168‑1 (MQ=255)
                                                                                                                                                                                                                               cgcCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCGTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCgg             <  1:566429/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aaCGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAACGGAAGTCGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCCTCGACATCGGCCCgctggacg  >  2:312816/1‑244 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
GCTGGGCAATCAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACA‑GGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGCAGGCGGCTCTTGCCATCGACATCGGCCCGACGTACG  >  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: ‑