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

CGGTAAAGAGAGCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTA  >  NC_000913/2374453‑2374883
                                                                                                                                                                                                                           |                                                                                                                                                                                                                     
ccgTAAAGAGAGCAACCAATCATCGTGACCCACAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTccc                                                                                                                                                                                                               <  2:252576/225‑1 (MQ=255)
          agCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCCTGTGCGCCTTCAACCACCTGAAAGTGCCGCAGCGTTAATGCCGGGAACGGGTTTTCCCCATCCGGccc                                                                                                                                                                                                     >  2:640513/1‑225 (MQ=255)
          agCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTTAACCAGCTGAATGTGCAGGGACGTTAATTCCGGAAACGGGCTTTCCCCATCCGGCAc                                                                                                                                                                                                     >  2:41632/1‑227 (MQ=255)
          agCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAAACAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGGGTTTTCCCCAACCGGccc                                                                                                                                                                                                     >  2:591611/1‑225 (MQ=255)
          agCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGGGCGCCTTCAACCAGCGGCAAGTGCAGCGACGTTAATGCCGGGAACCGGGTTTCCCCATCCCGCAcc                                                                                                                                                                                                    >  2:68687/1‑228 (MQ=255)
           acAACCAATCCCTGTGATCGCCAAACGGAATTCGCGATAACCAACCTTTCGGACTGGCAAGATGGGCCTGGTAAGGATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCa                                                                                                                                                                                                   <  1:185123/227‑1 (MQ=255)
                  atcatcGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCATGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGTATGTGCCGAGACCGTTAAGCCGGAAACCTGTTTTTCCCATCTCGCACCAGCGCaaa                                                                                                                                                                                            >  2:257960/1‑228 (MQ=255)
                                              gATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATTGGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGACCCCAAGCAAATTGCGGCGTCGGATTGGGCAGCAAttct                                                                                                                                                               >  2:71820/1‑228 (MQ=255)
                                                 aaCACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGGTAATGCCGGGAACGTGTTTTCTCCATCCCGCCCCAAGGCAAAATGCACGGCCAGATTGGGCATGCAttcttc                                                                                                                                                             >  2:639569/1‑227 (MQ=255)
                                                       ccTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGCCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATATGCAGCGCCGGTAAAGCCGGAAACTGGTTTTCCCCACCCGGCACCACACCAACAAGCAACGCCAGAAAGGGCAAAAAAACTTCaacaaaa                                                                                                                                                      >  1:285983/1‑223 (MQ=255)
                                                          tGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGGAACGCGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCCATTCTTCAAGCAACCGt                                                                                                                                                   >  2:499396/1‑228 (MQ=255)
                                                              ggcagagttgtataagagcaaaaGTAGGGACAGCGGCTTTCGGCAAACCGGCAGAACCCGAGGTAAACTGCCTTGAGCACAGACGGGTCGGCTGCCACGTAGCGGCATGGGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTCCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCg                                                                                                                                               <  1:569973/207‑1 (MQ=255)
                                                                 tGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGATGAATGTGCAGCGCCGCTAATGCCGGAAACGGGTTTTTCCCACCCCGCCCCAGAGAACATTGCACGCGCAGATTGGGCAGCAATTTTTCAAGCAAccggcgcgcca                                                                                                                                            >  1:492906/1‑218 (MQ=255)
                                                                  ggCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGGAACGGGTTTTCCCCATCCGGCCCCAGAGCAAATTGCAGCGCAAGATTGGGCCGCAATTCTTCAAGACACGGTTGCGGCAg                                                                                                                                           >  2:300254/1‑228 (MQ=255)
                                                                     aaGGTGGGCCTGAAAAGTATGGACAGCGGCTTTCGGCAAACCGGTAAAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTg                                                                                                                                        <  1:640513/228‑1 (MQ=255)
                                                                             ccTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGGAACGTGGTTTCCCCATCCGGCACCAGAGCCAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCacgtcgggggtc                                                                                                                                 >  2:656640/1‑216 (MQ=255)
                                                                             ccTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGCGCGCCTTACACCCGCCGAATGTGCAGCGACGTTACTGCCGGAAACGGGTTTTCCCCACCCGGCACCACAACACATTCCAGCGTGCGATTGTGCCCCAAATCTTCCAACAACGCTTCGGGCAGATGAGGGGTc                                                                                                                                 >  1:437852/1‑227 (MQ=255)
                                                                              cggatagttATGGACAGCGGGTTTCGGCAAACCGGTAGGACCCGAGGTCAACGTCCTTGAGCCCAGGCGCGTCGGCTGCCACGGTGCGGCATGTGCGCCTTCAACCAGCTGAATGTTCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTTCGGCAGCTGAGGGTTCa                                                                                                                                <  1:105151/219‑1 (MQ=255)
                                                                                           aCAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGGGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCGAATTGGGCAGCAATTCTTCAAGCAACGGATGCGGCAGCGGAGGGGTCCCGGGCAAAAACcgg                                                                                                                  >  1:176828/1‑227 (MQ=255)
                                                                                            actgcgttttCGGCAAACCGGTTAAACCCCGTGGCAACCGCCTTGAGCACCGGCGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCcgcg                                                                                                                 <  1:300254/221‑1 (MQ=255)
                                                                                             aGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGCGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCccgc                                                                                                                >  2:157422/1‑228 (MQ=255)
                                                                                              cgcgCTTTCGGCACACCGGTAGACCCCGGGGGCAACGTCATTGAGCACAGACGCGTCGGCTGCCCGGTAGCGGCATGTGCGCCGTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGcc                                                                                                               <  1:248841/225‑1 (MQ=255)
                                                                                               cGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGAAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGCgggcggtcccttttacacaccccc                                                                                                               >  2:569973/1‑206 (MQ=255)
                                                                                               cGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGGTAATGCCGGAAACGGGTTTTCCCCCTCCGGCCCCCGAGCAAATTGCAGGGTCCGATTGGGCAGCAATTCTTCCAGCAACGGTTGCGGCCGCCGAGGGTTTACGGGGAaccacccccccc                                                                                                              >  1:398694/1‑218 (MQ=255)
                                                                                               cGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGGGCGCCTTCAACCAGCTGAATGGGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGAGTCCGCTTGGGCAGCAATTCTTCAAACAACGGGTTCGGCAGCTGAGGGTTTACGGGCAACAccccacccc                                                                                                              >  2:433149/1‑222 (MQ=255)
                                                                                                ggCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGGATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCCGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGGTCACCGGCAACACCCGCGcccc                                                                                                             >  2:385982/1‑228 (MQ=255)
                                                                                                              cGGTACAACCCGGGGTCACCGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAGCGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACg                                                                                               <  1:378867/228‑1 (MQ=255)
                                                                                                                               aaCGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGGGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGGGTTTTCCCCCTCCCGCACCCAACCAAATTGCCACGTCAGATTGGGCAGCAATTCTTCAACCAACGGGTGCGGCAGCCTCGGGGTCCAGGGCAAACCCCCCACCCCCCCATGCAGTaaaccccccccagccgaaaga                                                                              >  1:59821/1‑209 (MQ=255)
                                                                                                                                    ccttggaccacgccgcgTCGGCTGCCACGGAGGGGCATTGGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTg                                                                         <  1:656640/215‑1 (MQ=255)
                                                                                                                                         ggCACGGACGCGTTGCCTCCCACGTAGCCGCATGTGCGCCGTCAACCAGCTGAATGGGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACCCCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTgcggcg                                                                    <  1:385982/227‑1 (MQ=255)
                                                                                                                                         agcccaggcggctCGGGTTCCCCGTGGCGGCATGGGCGCCTTCAACCAGCGGAATGTGCGGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCGGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTgcggcg                                                                    <  1:71820/216‑1 (MQ=255)
                                                                                                                                           cacaGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGCTTTCGGCAGCCTAGGGTTCACGGGCAACACCCGACCCCCCCCTTTCAGTAAACCCAGCCCGCCGCGCAACGttttgggggt                                                                   >  1:260534/1‑220 (MQ=255)
                                                                                                                                               gACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTcc                                                              <  2:390200/228‑1 (MQ=255)
                                                                                                                                                 cggctCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCAc                                                            <  2:533016/224‑1 (MQ=255)
                                                                                                                                                        gCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGTCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCCCGGGCAACACCCCCGCCCCCCCTTGCAGTAACGCCAGCCAGGCGAGCACCGTTTTGGGCGCATTCCCCGcccccc                                                     >  1:570861/1‑225 (MQ=255)
                                                                                                                                                               cgttgcggcctggtcggtttCAACCCGCCTAATGTGCAGCGTCTTTAATTGCGGAAACGGTTTTTCCCCATCCGGCACCAGGGCAAATTGCAGCGTTAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCGAGCCAGGCGAGCCGCGTTTGCGCCGTATTCCACGCCCGAAACATCa                                               <  1:499396/209‑1 (MQ=255)
                                                                                                                                                                         tggcccccTGCAACCAGCTGAAGGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCGATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGcc                                     <  2:378108/221‑1 (MQ=255)
                                                                                                                                                                                cttctaccatcggaattgtgaGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCCCCAGAGCAAATTGCAGCGTCGGATTGGGGAGCAATTCGTCAAGCAACGGTTGCGGCAGCTGAGGGTTCCCGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAacc                              <  1:536456/207‑1 (MQ=255)
                                                                                                                                                                                    ccacaaggtgaaggtagagcaccGTTAAGTCTGGAAACGGTCTTTCCCCATCCGGCCCCAGTGCAACGTGCAGCGTCAGTTTGGGCAGCAATTCTACAAGCCAGGGTTGCGGCATGCGGGGGTTCACGGGCAACCCCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCTTATTCCACGCCCGAAACATCACGCCGCTGCCCTCAaccacc                           <  1:506054/205‑1 (MQ=255)
                                                                                                                                                                                          cTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTACGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTCACGCCAGCCAGCCGAACAGCGTTTGCGGCGTATTCCACGCCCGCCACATCACGCCGCTGCACCCAACCACCCCCTGCAc                   >  2:506914/1‑228 (MQ=255)
                                                                                                                                                                                              aTGGGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCCCCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCGTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCa               <  2:398694/228‑1 (MQ=255)
                                                                                                                                                                                                    agagACGTTTATGGCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCg         <  1:288153/225‑1 (MQ=255)
                                                                                                                                                                                                     gCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGGGCAGCGTTTGCCGCGTATTCCCCGCCCGCAACACCAACCCGCGACCCCCACCCCCCCCCTTCCCCGGCAATccca        >  2:308884/1‑226 (MQ=255)
                                                                                                                                                                                                     agggcgtttAATGCCGGAAACGTGTTTCCCCCAGCCGGCCCCGGAGCAAGTTGCAGCGTACGGTTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCCGCTGAGGGTTCCCGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCg         <  1:68687/221‑1 (MQ=255)
                                                                                                                                                                                                            ttAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCCCTGCCCTCAACCACCCCCTGCCCCGCAAATCCGGGGGcca  >  2:588868/1‑226 (MQ=255)
                                                                                                                                                                                                                           |                                                                                                                                                                                                                     
CGGTAAAGAGAGCAACCAATCATCGTGATCGCCAAACGGAATCAGCGATAACACACCTTGCGCACTGGCAAGATGGGCCTGATAAGTATGGACAGCGGCTTTCGGCAAACCGGTAGAACCCGAGGTCAACGTCATTGAGCACAGACGCGTCGGCTGCCACGTAGCGGCATGTGCGCCTTCAACCAGCTGAATGTGCAGCGACGTTAATGCCGGAAACGTGTTTTCCCCATCCGGCACCAGAGCAAATTGCAGCGTCAGATTGGGCAGCAATTCTTCAAGCAACGGTTGCGGCAGCTGAGGGTTCACGGGCAACACCCGCGCCCCGCATTGCAGTAACGCCAGCCAGGCGAGCAGCGTTTGCGGCGTATTCCACGCCCGCAACATCACGCCGCTGCCCTCAACCACCCCCTGCACCGCAAATCCGGAGGCTA  >  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: