Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,374,671 T→G 24.4% T111P (ACG→CCG)  menE ← O‑succinylbenzoate‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,374,6710TG24.4% 64.2 / 2.2 45T111P (ACG→CCG) menEO‑succinylbenzoate‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (13/18);  new base G (11/2);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 1.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTAAAGAGAGCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTA  >  NC_000913/2374453‑2374883
                                                                                                                                                                                                                           |                                                                                                                                                                                                                     
ccgTAAAGAGAGCAACCAATCATCGTGACCCACAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTccc                                                                                                                                                                                                               <  2:252576/225‑1 (MQ=255)
          agCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCCTGTGCGCCTTCAACCACCTGAAAGTGCCGCAG‑CGTTAATGCCGGGAACGGGTTTTCCCCATCCGGccc                                                                                                                                                                                                     >  2:640513/1‑225 (MQ=255)
          agCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTTAACCAGCTGAATGTGCAGGGA‑CGTTAATTCCGGAAACGGGCTTTCCCCATCCGGCAc                                                                                                                                                                                                     >  2:41632/1‑227 (MQ=255)
          agCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAAACAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGGGTTTTCCCCAACCGGccc                                                                                                                                                                                                     >  2:591611/1‑225 (MQ=255)
          agCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGGGCGCCTTCAACCAGCGGCAAGTGCAGCGA‑CGTTAATGCCGGGAACCGGGTTTCCCCATCCCGCAcc                                                                                                                                                                                                    >  2:68687/1‑228 (MQ=255)
           acAACCAATCCCTGTGATCGCCAAACGGAATTCGCGATAACCAACCTTTCGGACTGGCAAGATGGGCCTGGTAAGGATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCa                                                                                                                                                                                                   <  1:185123/227‑1 (MQ=255)
                  atcatcGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCATGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGTATGTGCCGAGACCGTTAA‑GCCGGAAACCTGTTTTTCCCATCTCGCACCAGCGCaaa                                                                                                                                                                                            >  2:257960/1‑228 (MQ=255)
                                              gATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATTGGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGACCCCAAGCAAATTGCG‑GCGTCGGATTGGGCAGCAAttct                                                                                                                                                               >  2:71820/1‑228 (MQ=255)
                                                 aaCACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGGTAATGCCGGGAACGTGTTTTCTCCATCCCGCCCCAAGGCAAAATGCA‑CGGCCAGATTGGGCATGCAttcttc                                                                                                                                                             >  2:639569/1‑227 (MQ=255)
                                                       ccTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGCCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATATGCAGCGC‑CGGTAAAGCCGGAAACTGGTTTTCCCCACCCGGCACCACACCAACAAGCA‑ACGCCAGAAAGGGCAAAAAAACTTCaacaaaa                                                                                                                                                      >  1:285983/1‑223 (MQ=255)
                                                          tGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGGAACGCGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCCATTCTTCAAGCAACCGt                                                                                                                                                   >  2:499396/1‑228 (MQ=255)
                                                              ggcagagttgtataagagcaaaaGTAGGGACAGCGGCTTTCGGCAAACCGGCAGAACCCGAGGTAAACTGCCTTGAGCACAGACGGGTCGGCTGCCACGTAGCGGCATGGGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTCCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCg                                                                                                                                               <  1:569973/207‑1 (MQ=255)
                                                                 tGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGATGAATGTGCAGCGC‑CGCTAATGCCGGAAACGGGTTTTTCCCACCCCGCCCCAGAGAACATTGCACGCG‑CAGATTGGGCAGCAATTTTTCAAGCAAccggcgcgcca                                                                                                                                            >  1:492906/1‑218 (MQ=255)
                                                                  ggCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGGAACGGGTTTTCCCCATCCGGCCCCAGAGCAAATTGCA‑GCGCAAGATTGGGCCGCAATTCTTCAAGACACGGTTGCGGCAg                                                                                                                                           >  2:300254/1‑228 (MQ=255)
                                                                     aaGGTGGGCCTGAAAAGTATGGACAGCGGCTTTCGGCAAACCGGTAAAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTg                                                                                                                                        <  1:640513/228‑1 (MQ=255)
                                                                             ccTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGGAACGTGGTTTCCCCATCCGGCACCAGAGCCAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCacgtcgggggtc                                                                                                                                 >  2:656640/1‑216 (MQ=255)
                                                                             ccTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGCGCGCCTTACACCCGCCGAATGTGCAGCGA‑CGTTACTGCCGGAAACGGGTTTTCCCCACCCGGCACCACAACACATTCCA‑GCGTGCGATTGTGCCCCAAATCTTCCAACAACGCTTCGGGCAGATGAGGGGTc                                                                                                                                 >  1:437852/1‑227 (MQ=255)
                                                                              cggatagttATGGACAGCGGGTTTCGGCAAACCGGTAGGACCCGAGGTCAACGTCCTTGAGCCCAGGCGCGTCGGCTGCCACGGTGCGGCATGTGCGCCTTCAACCAGCTGAATGTTCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTTCGGCAGCTGAGGGTTCa                                                                                                                                <  1:105151/219‑1 (MQ=255)
                                                                                           aCAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGGGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCGAATTGGGCAGCAATTCTTCAAGCAACGGATGCGGCAGCGGAGGGGTCCCGGGCAAAAACcgg                                                                                                                  >  1:176828/1‑227 (MQ=255)
                                                                                            actgcgttttCGGCAAACCGGTTAAACCCCGTGGCAACCGCCTTGAGCACCGGCGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCcgcg                                                                                                                 <  1:300254/221‑1 (MQ=255)
                                                                                             aGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGCGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCccgc                                                                                                                >  2:157422/1‑228 (MQ=255)
                                                                                              cgcgCTTTCGGCACACCGGTAGACCCCGGGGGCAACGTCATTGAGCACAGACGCGTCGGCTGCCCGGTAGCGGCATGTGCGCCGTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGcc                                                                                                               <  1:248841/225‑1 (MQ=255)
                                                                                               cGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGAAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGCgggcggtcccttttacacaccccc                                                                                                               >  2:569973/1‑206 (MQ=255)
                                                                                               cGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGGTAATGCCGGAAACGGGTTTTCCCCCTCCGGCCCCCGAGCAAATTGCA‑GGGTCCGATTGGGCAGCAATTCTTCCAGCAACGGTTGCGGCCGCCGAGGGTTTACGGGGAaccacccccccc                                                                                                              >  1:398694/1‑218 (MQ=255)
                                                                                               cGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGGGCGCCTTCAACCAGCTGAATGGGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GAGTCCGCTTGGGCAGCAATTCTTCAAACAACGGGTTCGGCAGCTGAGGGTTTACGGGCAACAccccacccc                                                                                                              >  2:433149/1‑222 (MQ=255)
                                                                                                ggCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGGATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCCGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGGTCACCGGCAACACCCGCGcccc                                                                                                             >  2:385982/1‑228 (MQ=255)
                                                                                                              cGGTACAACCCGGGGTCACCGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAGCGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACg                                                                                               <  1:378867/228‑1 (MQ=255)
                                                                                                                               aaCGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGGGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGGGTTTTCCCCCTCCCGCACCCAACCAAATTGCC‑ACGTCAGATTGGGCAGCAATTCTTCAACCAACGGGTGCGGCAGCCTCGGGGTCCAGGGCAAACCCCCCACCCCCCCATGCAGTaaaccccccccagccgaaaga                                                                              >  1:59821/1‑209 (MQ=255)
                                                                                                                                    ccttggaccacgccgcgTCGGCTGCCACGGAGGGGCATTGGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTg                                                                         <  1:656640/215‑1 (MQ=255)
                                                                                                                                         ggCACGGACGCGTTGCCTCCCACGTAGCCGCATGTGCGCCGTCAACCAGCTGAATGGGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACCCCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTgcggcg                                                                    <  1:385982/227‑1 (MQ=255)
                                                                                                                                         agcccaggcggctCGGGTTCCCCGTGGCGGCATGGGCGCCTTCAACCAGCGGAATGTGCGGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCG‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTgcggcg                                                                    <  1:71820/216‑1 (MQ=255)
                                                                                                                                           cacaGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGCTTTCGGCAGCCTAGGGTTCACGGGCAACACCCGACCCCCCCCTTTCAGTAAACCCAGCCCGCCGCGCAACGttttgggggt                                                                   >  1:260534/1‑220 (MQ=255)
                                                                                                                                               gACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTcc                                                              <  2:390200/228‑1 (MQ=255)
                                                                                                                                                 cggctCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCAc                                                            <  2:533016/224‑1 (MQ=255)
                                                                                                                                                        gCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGTCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCCCGGGCAACACCCCCGCCCCCCCTTGCAGTAACGCCAGCCAGGCGAGCACCGTTTTGGGCGCATTCCCCGcccccc                                                     >  1:570861/1‑225 (MQ=255)
                                                                                                                                                               cgttgcggcctggtcggtttCAACCCGCCTAATGTGCAGCGT‑CTTTAATTGCGGAAACGGTTTTTCCCCATCCGGCACCAGGGCAAATTGCA‑GCGTTAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCGAGCCAGGCGAGCCGCGTTTGCGCCGTATTCCACGCCCGAAACATCa                                               <  1:499396/209‑1 (MQ=255)
                                                                                                                                                                         tggcccccTGCAACCAGCTGAAGGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCGATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGcc                                     <  2:378108/221‑1 (MQ=255)
                                                                                                                                                                                cttctaccatcggaattgtgaGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCCCCAGAGCAAATTGCA‑GCGTCGGATTGGGGAGCAATTCGTCAAGCAACGGTTGCGGCAGCTGAGGGTTCCCGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAacc                              <  1:536456/207‑1 (MQ=255)
                                                                                                                                                                                    ccacaaggtgaaggtagagcaccGTTAAGTCTGGAAACGGTCTTTCCCCATCCGGCCCCAGTGCAACGTGCA‑GCGTCAGTTTGGGCAGCAATTCTACAAGCCAGGGTTGCGGCATGCGGGGGTTCACGGGCAACCCCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCTTATTCCACGCCCGAAACATCACGCCGCTGCCCTCAaccacc                           <  1:506054/205‑1 (MQ=255)
                                                                                                                                                                                          cTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTACGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTCACGCCAGCCAGCCGAACAGCGTTTGCGGCGTATTCCACGCCCGCCACATCACGCCGCTGCACCCAACCACCCCCTGCAc                   >  2:506914/1‑228 (MQ=255)
                                                                                                                                                                                              aTGGGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCCCCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCGTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCa               <  2:398694/228‑1 (MQ=255)
                                                                                                                                                                                                    agagA‑CGTTTATGGCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCg         <  1:288153/225‑1 (MQ=255)
                                                                                                                                                                                                     gCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGGGCAGCGTTTGCCGCGTATTCCCCGCCCGCAACACCAACCCGCGACCCCCACCCCCCCCCTTCCCCGGCAATccca        >  2:308884/1‑226 (MQ=255)
                                                                                                                                                                                                     agggcgtttAATGCCGGAAACGTGTTTCCCCCAGCCGGCCCCGGAGCAAGTTGCA‑GCGTACGGTTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCCGCTGAGGGTTCCCGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCg         <  1:68687/221‑1 (MQ=255)
                                                                                                                                                                                                            ttAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCCCTGCCCTCAACCACCCCCTGCCCCGCAAATCCGGGGGcca  >  2:588868/1‑226 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                     
CGGTAAAGAGAGCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGA‑CGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCA‑GCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTA  >  NC_000913/2374453‑2374883

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑