Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,138,152 C→T T261I (ACC→ATC)  flgJ → putative peptidoglycan hydrolase FlgJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,138,1520CT98.6% 209.0 / NA 70T261I (ACC→ATC) flgJputative peptidoglycan hydrolase FlgJ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (32/37);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGACCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAAT  >  NC_000913/1137906‑1138395
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
cacgggAATCTGGTTGGGGCAAACGGCAACTCCGCCGCGAAAACGGTGAGCCGAGCGATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATTAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTACCCCTCGC‑TACGCCGCCGTGACGAtcgcc                                                                                                                                                                                                                                                  <  2:270064/247‑1 (MQ=255)
cACTGGAATCTGGTGGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGAtcgcc                                                                                                                                                                                                                                                  <  1:128576/251‑1 (MQ=255)
        tCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGc                                                                                                                                                                                                                                          >  1:233591/1‑251 (MQ=255)
           gttggggggCACCGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑a                                                                                                                                                                                                                                      <  1:163626/247‑1 (MQ=255)
                 ggacAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACATGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCTCGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAggg                                                                                                                                                                                                                                <  1:518137/248‑1 (MQ=255)
                  ggCAAGGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAgggg                                                                                                                                                                                                                               <  1:164002/251‑1 (MQ=255)
                     aaCGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGGTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAACATTGTCGGATTACGTTGGGCTGTTAACGAGTAAACCGCGC‑TACGCCGCCGTGACGTTCGACGCGAGTGGGG‑AACAGGGGgcg                                                                                                                                                                                                                            >  2:633648/1‑251 (MQ=255)
                       cGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCa                                                                                                                                                                                                                          <  1:395251/251‑1 (MQ=255)
                            aaTCCGCCGCGAAAACGGCGAGCCGCGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACAGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGATCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCATAACCCGCGC‑TACGCCACCGTGACGATCGCCGCGAGTGGGG‑AACGGTGGGCGAAGGccc                                                                                                                                                                                                                     >  2:571744/1‑251 (MQ=255)
                                 gccgcGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGGAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGGGCGGCAAC‑GGGGGCGCCGGCCCtcacg                                                                                                                                                                                                                >  2:193247/1‑247 (MQ=255)
                                     cgAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGAcg                                                                                                                                                                                                            <  1:546767/251‑1 (MQ=255)
                                                         ataACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATcc                                                                                                                                                                                         <  1:257439/250‑1 (MQ=255)
                                                           aaCCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTc                                                                                                                                                                                       >  1:526468/1‑250 (MQ=255)
                                                           aaCCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGCGCTATGCCACCGATCCTcc                                                                                                                                                                                      >  2:230581/1‑250 (MQ=255)
                                                            aCCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCAc                                                                                                                                                                                     >  1:671268/1‑251 (MQ=255)
                                                             ccTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTccct                                                                                                                                                                                    >  1:250871/1‑248 (MQ=255)
                                                               tGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTccctct                                                                                                                                                                                  >  2:333474/1‑246 (MQ=255)
                                                                gttcggtttCAAAGCTTCGGGAAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTAAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCGACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTa                                                                                                                                                                                   <  2:553603/241‑1 (MQ=255)
                                                                gTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTa                                                                                                                                                                                   <  1:59398/249‑1 (MQ=255)
                                                                gTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATg                                                                                                                                                                                 >  2:189999/1‑251 (MQ=255)
                                                                          aaaGCCTCTGGCAACTGGACAGGGCCAGTTACTGAAATCACCCCGACTGAATATGAAAACGGCGAAGCGAAGAAAGAAAAAGCGAAGTTTCGCGTCTAAAGCTAGTATCTGGAAGCCTTGTCGGATTACGTTGGGATGTTAACGCGTCACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGAGCGG‑AAAGGGGGGCGCCGGCCCTACAGGACGCGGGCAATGCCACCGATCCTCACTATACCCCCAAAc                                                                                                                                                                        >  2:403926/1‑250 (MQ=255)
                                                                          aaaGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACt                                                                                                                                                                       >  1:209083/1‑251 (MQ=255)
                                                                                         tGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCa                                                                                                                                                                                      <  1:526597/221‑1 (MQ=255)
                                                                                                  ccAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGCATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGCTTCAGCAGATGa                                                                                                                                                >  1:84325/1‑250 (MQ=255)
                                                                                                        aCTGAAATCACCACGACTGAATATGAAAACGGCGCAGCGAAGACAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGAATACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTCTGCCCGCAAACTCACCAACATGATTCAGCAGATGAACTCGAt                                                                                                                                         >  2:82851/1‑251 (MQ=255)
                                                                                                        aCTGAAACCACCACGGTCGAATATGAAAACGGCGAATCGAAGAAAGTAAAAGCGAAGTTTCGCGTTTACATCTCGTATCTGGAACCCTTGTCGGATTCCGTTTGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGCCGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGa                                                                                                                                          <  2:246766/250‑1 (MQ=255)
                                                                                                         cggaaacccccccggttGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCTTCTACATCTCGTATCTGGAAGCCTTGGCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATa                                                                                                                                        <  2:547810/235‑1 (MQ=255)
                                                                                                          tGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATa                                                                                                                                        >  1:489790/1‑250 (MQ=255)
                                                                                                           gAAATCACCACGACTGAAAATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGTGTTTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGGTGGGCTGTTACCGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATTAAATCGATa                                                                                                                                        <  2:680632/249‑1 (MQ=255)
                                                                                                           gAAAACACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTGCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGGTGGGCTGTTAACGCGTAACCCGGGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATaa                                                                                                                                       <  1:396992/250‑1 (MQ=255)
                                                                                                            aaaTCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGc                                                                                                                                     >  1:185552/1‑251 (MQ=255)
                                                                                                                caccacGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTAGCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGAc                                                                                                                                  <  1:406705/250‑1 (MQ=255)
                                                                                                                      gACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTGGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAg                                                                                                                                                                                                                <  1:570248/166‑1 (MQ=255)
                                                                                                                      gACTGAATATGAAAAAGGCGCAGCGAAGAAAGTAACAGCGAAGTTTCGAGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGAGGGCGCAGGACCTTCAg                                                                                                                                                                                                                >  2:570248/1‑166 (MQ=255)
                                                                                                                             cacgaaaactggcaatcgaagAAAGACAAAGTGAAGTGTGGCGTCTACAGCCGGGAGCTGGAACCCTTGACGTATTACGTGGGCCTGTTAGCGCGTAACC‑GTGCCTACTCCGCCGCGACGACCGCCGCGTGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTAGGCCACCGATCCTCACTATGCCCGCAAACTCGCAAACATGATGCAGCAGAGGAAAGCGCTAAGCCACAAGGTGAGcggaa                                                                                                                     <  2:301341/234‑5 (MQ=255)
                                                                                                                                gAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCt                                                                                                                  >  1:481879/1‑250 (MQ=255)
                                                                                                                                 aaaaCGGCGAAGCGAACAAAGTAAAAGCGAAGTTGCGCGTTGACAGCTCGTTACTGGAAGCCTAGTCGTATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACATGGTGAGCAAATCCTAc                                                                                                                <  2:150061/251‑1 (MQ=255)
                                                                                                                                            gcgaagAAATTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACAt                                                                                                     <  2:489790/251‑1 (MQ=255)
                                                                                                                                               aagaaAGTAAAAGCGAAGGTTCGCGTCTACAGCTCGTATCGGGAAGCCTTGACGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCGACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGa                                                                                                  <  1:315208/251‑1 (MQ=255)
                                                                                                                                                agaaAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTACCCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGAt                                                                                                 <  1:63409/251‑1 (MQ=255)
                                                                                                                                                agaaAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGACCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGAGACAAGGTGAGCAAAACCTACAGTATGAACATTGAt                                                                                                 >  2:478712/1‑251 (MQ=255)
                                                                                                                                                   aaGTAAAAGCTAAGTTGCGAGTCTGCAGCCCGAATCAGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑GACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCGCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAt                                                                                              <  2:255959/251‑1 (MQ=255)
                                                                                                                                                   aaGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAt                                                                                              <  1:193247/251‑1 (MQ=255)
                                                                                                                                                      aaaaatggatgTTTCGCGTCTCCAGCTCTTATCGGGAAGCCTTGTCGGTTTACGTTGGGCTGTTAACTCGTAACCCGCGC‑TCCGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCt                                                                                            <  2:569945/240‑1 (MQ=255)
                                                                                                                                                       aaaaGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGt                                                                                          >  1:133728/1‑251 (MQ=255)
                                                                                                                                                        aaaGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGGGCAAAACCTACAGTATGAACCTTTATAATCTGtt                                                                                         >  2:561995/1‑251 (MQ=255)
                                                                                                                                                        aaaGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGtt                                                                                         <  1:614261/251‑1 (MQ=255)
                                                                                                                                                         aaGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTAAGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTc                                                                                        >  1:528039/1‑251 (MQ=255)
                                                                                                                                                              aaGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAAt                                                                                   <  1:559829/251‑1 (MQ=255)
                                                                                                                                                                gTTTCGCGTTTACAGCCCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTACCCCGCGC‑TCCGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATaa                                                                                 <  2:526468/251‑1 (MQ=255)
                                                                                                                                                                    cgcgTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCa                                                                             <  1:566448/251‑1 (MQ=255)
                                                                                                                                                                        tCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTc                                                                         >  2:519613/1‑251 (MQ=255)
                                                                                                                                                                           aCAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTcc                                                                        <  1:230581/249‑1 (MQ=255)
                                                                                                                                                                             aGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTCCAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCcggcg                                                                    >  1:453967/1‑251 (MQ=255)
                                                                                                                                                                                 cGTATCTGGAAGCCTGGTCGGCTTACGGTGGGCTGTTAACGAGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGGCGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTc                                                                <  1:251021/251‑1 (MQ=255)
                                                                                                                                                                                    aTCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGc                                                              <  1:554875/250‑1 (MQ=255)
                                                                                                                                                                                        ggAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTg                                                            <  1:403926/248‑1 (MQ=255)
                                                                                                                                                                                        ggAAGCCTTGTCGGATTACGTTGGGCTATTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGcc                                                          <  1:226963/250‑1 (MQ=255)
                                                                                                                                                                                         gAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCg                                                         >  1:313207/1‑250 (MQ=255)
                                                                                                                                                                                           aGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGATATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGCCCCTGCCGat                                                       >  1:295903/1‑250 (MQ=255)
                                                                                                                                                                                                      aTTACGTTGGGCTGTTAACGCGTAACCCGCTA‑TCCGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACATGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCACACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACaat                                                                                                                               <  2:406533/167‑2 (MQ=255)
                                                                                                                                                                                                      aTTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAt                                                                                              <  1:433488/200‑1 (MQ=255)
                                                                                                                                                                                                      aTTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAt                                                                                              >  2:433488/1‑200 (MQ=255)
                                                                                                                                                                                                      aTTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGACGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAg                                                                                                                               >  1:406533/1‑167 (MQ=255)
                                                                                                                                                                                                                gCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTAc                                                                                                                >  2:129543/1‑172 (MQ=255)
                                                                                                                                                                                                                gCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTAc                                                                                                                <  1:129543/172‑1 (MQ=255)
                                                                                                                                                                                                                   gTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCCGATGAAATCGATAAGCGCCAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGGTAATACTCTGTAATTGAAGGCTtat                               >  2:628151/1‑250 (MQ=255)
                                                                                                                                                                                                                             aaCCCGCGC‑TACGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTcc                    >  1:169544/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aCGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCt          <  1:551426/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aCGCCGCCGTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAAGACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCt          <  1:75406/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cgccgTGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGAt      <  1:109683/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 tGACGATCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATtaat  >  1:13508/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
CACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGC‑TACGCCGCCGTGACGACCGCCGCGAGTGCGG‑AACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAAT  >  NC_000913/1137906‑1138395

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: ‑