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

CACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAAT  >  NC_000913/1137906‑1138395
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
cacgggAATCTGGTTGGGGCAAACGGCAACTCCGCCGCGAAAACGGTGAGCCGAGCGATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATTAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTACCCCTCGCTACGCCGCCGTGACGAtcgcc                                                                                                                                                                                                                                                  <  2:270064/247‑1 (MQ=255)
cACTGGAATCTGGTGGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGAtcgcc                                                                                                                                                                                                                                                  <  1:128576/251‑1 (MQ=255)
        tCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGc                                                                                                                                                                                                                                          >  1:233591/1‑251 (MQ=255)
           gttggggggCACCGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGa                                                                                                                                                                                                                                      <  1:163626/247‑1 (MQ=255)
                 ggacAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACATGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCTCGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAggg                                                                                                                                                                                                                                <  1:518137/248‑1 (MQ=255)
                  ggCAAGGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAgggg                                                                                                                                                                                                                               <  1:164002/251‑1 (MQ=255)
                     aaCGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGGTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAACATTGTCGGATTACGTTGGGCTGTTAACGAGTAAACCGCGCTACGCCGCCGTGACGTTCGACGCGAGTGGGGAACAGGGGgcg                                                                                                                                                                                                                            >  2:633648/1‑251 (MQ=255)
                       cGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCa                                                                                                                                                                                                                          <  1:395251/251‑1 (MQ=255)
                            aaTCCGCCGCGAAAACGGCGAGCCGCGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACAGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGATCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCATAACCCGCGCTACGCCACCGTGACGATCGCCGCGAGTGGGGAACGGTGGGCGAAGGccc                                                                                                                                                                                                                     >  2:571744/1‑251 (MQ=255)
                                 gccgcGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGGAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGGGCGGCAACGGGGGCGCCGGCCCtcacg                                                                                                                                                                                                                >  2:193247/1‑247 (MQ=255)
                                     cgAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGAcg                                                                                                                                                                                                            <  1:546767/251‑1 (MQ=255)
                                                         ataACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATcc                                                                                                                                                                                         <  1:257439/250‑1 (MQ=255)
                                                           aaCCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTc                                                                                                                                                                                       >  1:526468/1‑250 (MQ=255)
                                                           aaCCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGCGCTATGCCACCGATCCTcc                                                                                                                                                                                      >  2:230581/1‑250 (MQ=255)
                                                            aCCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCAc                                                                                                                                                                                     >  1:671268/1‑251 (MQ=255)
                                                             ccTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTccct                                                                                                                                                                                    >  1:250871/1‑248 (MQ=255)
                                                               tGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTccctct                                                                                                                                                                                  >  2:333474/1‑246 (MQ=255)
                                                                gttcggtttCAAAGCTTCGGGAAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTAAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCGACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTa                                                                                                                                                                                   <  2:553603/241‑1 (MQ=255)
                                                                gTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTa                                                                                                                                                                                   <  1:59398/249‑1 (MQ=255)
                                                                gTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATg                                                                                                                                                                                 >  2:189999/1‑251 (MQ=255)
                                                                          aaaGCCTCTGGCAACTGGACAGGGCCAGTTACTGAAATCACCCCGACTGAATATGAAAACGGCGAAGCGAAGAAAGAAAAAGCGAAGTTTCGCGTCTAAAGCTAGTATCTGGAAGCCTTGTCGGATTACGTTGGGATGTTAACGCGTCACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGAGCGGAAAGGGGGGCGCCGGCCCTACAGGACGCGGGCAATGCCACCGATCCTCACTATACCCCCAAAc                                                                                                                                                                        >  2:403926/1‑250 (MQ=255)
                                                                          aaaGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACt                                                                                                                                                                       >  1:209083/1‑251 (MQ=255)
                                                                                         tGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCa                                                                                                                                                                                      <  1:526597/221‑1 (MQ=255)
                                                                                                  ccAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGCATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGCTTCAGCAGATGa                                                                                                                                                >  1:84325/1‑250 (MQ=255)
                                                                                                        aCTGAAATCACCACGACTGAATATGAAAACGGCGCAGCGAAGACAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGAATACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTCTGCCCGCAAACTCACCAACATGATTCAGCAGATGAACTCGAt                                                                                                                                         >  2:82851/1‑251 (MQ=255)
                                                                                                        aCTGAAACCACCACGGTCGAATATGAAAACGGCGAATCGAAGAAAGTAAAAGCGAAGTTTCGCGTTTACATCTCGTATCTGGAACCCTTGTCGGATTCCGTTTGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGCCGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGa                                                                                                                                          <  2:246766/250‑1 (MQ=255)
                                                                                                         cggaaacccccccggttGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCTTCTACATCTCGTATCTGGAAGCCTTGGCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATa                                                                                                                                        <  2:547810/235‑1 (MQ=255)
                                                                                                          tGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATa                                                                                                                                        >  1:489790/1‑250 (MQ=255)
                                                                                                           gAAATCACCACGACTGAAAATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGTGTTTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGGTGGGCTGTTACCGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATTAAATCGATa                                                                                                                                        <  2:680632/249‑1 (MQ=255)
                                                                                                           gAAAACACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTGCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGGTGGGCTGTTAACGCGTAACCCGGGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATaa                                                                                                                                       <  1:396992/250‑1 (MQ=255)
                                                                                                            aaaTCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGc                                                                                                                                     >  1:185552/1‑251 (MQ=255)
                                                                                                                caccacGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTAGCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGAc                                                                                                                                  <  1:406705/250‑1 (MQ=255)
                                                                                                                      gACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTGGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAg                                                                                                                                                                                                                <  1:570248/166‑1 (MQ=255)
                                                                                                                      gACTGAATATGAAAAAGGCGCAGCGAAGAAAGTAACAGCGAAGTTTCGAGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGAGGGCGCAGGACCTTCAg                                                                                                                                                                                                                >  2:570248/1‑166 (MQ=255)
                                                                                                                             cacgaaaactggcaatcgaagAAAGACAAAGTGAAGTGTGGCGTCTACAGCCGGGAGCTGGAACCCTTGACGTATTACGTGGGCCTGTTAGCGCGTAACCGTGCCTACTCCGCCGCGACGACCGCCGCGTGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTAGGCCACCGATCCTCACTATGCCCGCAAACTCGCAAACATGATGCAGCAGAGGAAAGCGCTAAGCCACAAGGTGAGcggaa                                                                                                                     <  2:301341/234‑5 (MQ=255)
                                                                                                                                gAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCt                                                                                                                  >  1:481879/1‑250 (MQ=255)
                                                                                                                                 aaaaCGGCGAAGCGAACAAAGTAAAAGCGAAGTTGCGCGTTGACAGCTCGTTACTGGAAGCCTAGTCGTATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACATGGTGAGCAAATCCTAc                                                                                                                <  2:150061/251‑1 (MQ=255)
                                                                                                                                            gcgaagAAATTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACAt                                                                                                     <  2:489790/251‑1 (MQ=255)
                                                                                                                                               aagaaAGTAAAAGCGAAGGTTCGCGTCTACAGCTCGTATCGGGAAGCCTTGACGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCGACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGa                                                                                                  <  1:315208/251‑1 (MQ=255)
                                                                                                                                                agaaAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTACCCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGAt                                                                                                 <  1:63409/251‑1 (MQ=255)
                                                                                                                                                agaaAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGACCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGAGACAAGGTGAGCAAAACCTACAGTATGAACATTGAt                                                                                                 >  2:478712/1‑251 (MQ=255)
                                                                                                                                                   aaGTAAAAGCTAAGTTGCGAGTCTGCAGCCCGAATCAGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCGACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCGCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAt                                                                                              <  2:255959/251‑1 (MQ=255)
                                                                                                                                                   aaGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAt                                                                                              <  1:193247/251‑1 (MQ=255)
                                                                                                                                                      aaaaatggatgTTTCGCGTCTCCAGCTCTTATCGGGAAGCCTTGTCGGTTTACGTTGGGCTGTTAACTCGTAACCCGCGCTCCGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCt                                                                                            <  2:569945/240‑1 (MQ=255)
                                                                                                                                                       aaaaGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGt                                                                                          >  1:133728/1‑251 (MQ=255)
                                                                                                                                                        aaaGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGGGCAAAACCTACAGTATGAACCTTTATAATCTGtt                                                                                         >  2:561995/1‑251 (MQ=255)
                                                                                                                                                        aaaGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGtt                                                                                         <  1:614261/251‑1 (MQ=255)
                                                                                                                                                         aaGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTAAGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTc                                                                                        >  1:528039/1‑251 (MQ=255)
                                                                                                                                                              aaGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAAt                                                                                   <  1:559829/251‑1 (MQ=255)
                                                                                                                                                                gTTTCGCGTTTACAGCCCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTACCCCGCGCTCCGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATaa                                                                                 <  2:526468/251‑1 (MQ=255)
                                                                                                                                                                    cgcgTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCa                                                                             <  1:566448/251‑1 (MQ=255)
                                                                                                                                                                        tCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTc                                                                         >  2:519613/1‑251 (MQ=255)
                                                                                                                                                                           aCAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTcc                                                                        <  1:230581/249‑1 (MQ=255)
                                                                                                                                                                             aGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTCCAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCcggcg                                                                    >  1:453967/1‑251 (MQ=255)
                                                                                                                                                                                 cGTATCTGGAAGCCTGGTCGGCTTACGGTGGGCTGTTAACGAGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGGCGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTc                                                                <  1:251021/251‑1 (MQ=255)
                                                                                                                                                                                    aTCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGc                                                              <  1:554875/250‑1 (MQ=255)
                                                                                                                                                                                        ggAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTg                                                            <  1:403926/248‑1 (MQ=255)
                                                                                                                                                                                        ggAAGCCTTGTCGGATTACGTTGGGCTATTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGcc                                                          <  1:226963/250‑1 (MQ=255)
                                                                                                                                                                                         gAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCg                                                         >  1:313207/1‑250 (MQ=255)
                                                                                                                                                                                           aGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGATATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGCCCCTGCCGat                                                       >  1:295903/1‑250 (MQ=255)
                                                                                                                                                                                                      aTTACGTTGGGCTGTTAACGCGTAACCCGCTATCCGCCGCCGTGACGATCGCCGCGAGTGCGGAACATGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCACACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACaat                                                                                                                               <  2:406533/167‑2 (MQ=255)
                                                                                                                                                                                                      aTTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAt                                                                                              <  1:433488/200‑1 (MQ=255)
                                                                                                                                                                                                      aTTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAAt                                                                                              >  2:433488/1‑200 (MQ=255)
                                                                                                                                                                                                      aTTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGACGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAg                                                                                                                               >  1:406533/1‑167 (MQ=255)
                                                                                                                                                                                                                gCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTAc                                                                                                                >  2:129543/1‑172 (MQ=255)
                                                                                                                                                                                                                gCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTAc                                                                                                                <  1:129543/172‑1 (MQ=255)
                                                                                                                                                                                                                   gTTAACGCGTAACCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCCGATGAAATCGATAAGCGCCAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGGTAATACTCTGTAATTGAAGGCTtat                               >  2:628151/1‑250 (MQ=255)
                                                                                                                                                                                                                             aaCCCGCGCTACGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTcc                    >  1:169544/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aCGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCt          <  1:551426/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aCGCCGCCGTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAAGACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCt          <  1:75406/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cgccgTGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGAt      <  1:109683/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 tGACGATCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATtaat  >  1:13508/1‑250 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
CACTGGAATCTGGTTGGGGGCAACGGCAAATCCGCCGCGAAAACGGCGAGCCGAGCTATAACCTGTTTGGTGTCAAAGCCTCTGGCAACTGGAAAGGGCCAGTTACTGAAATCACCACGACTGAATATGAAAACGGCGAAGCGAAGAAAGTAAAAGCGAAGTTTCGCGTCTACAGCTCGTATCTGGAAGCCTTGTCGGATTACGTTGGGCTGTTAACGCGTAACCCGCGCTACGCCGCCGTGACGACCGCCGCGAGTGCGGAACAGGGGGCGCAGGCCCTACAGGACGCGGGCTATGCCACCGATCCTCACTATGCCCGCAAACTCACCAACATGATTCAGCAGATGAAATCGATAAGCGACAAGGTGAGCAAAACCTACAGTATGAACATTGATAATCTGTTCTGAATAACTCAAGTCCGGCGGGTCGCTGCCGATAATACTCTGTAATTGAAGGCTTATAAGGAACCTCCATGTCCAGCTTGATTAAT  >  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: