Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,691,173 A→T 5.7% noncoding (168/184 nt) glmY ← sRNA activator of glmS mRNA, glmZ processing antagonist

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,691,1730AT 5.7% 812.3 / 27.0 231noncoding (168/184 nt)glmYsRNA activator of glmS mRNA, glmZ processing antagonist
Reads supporting (aligned to +/- strand):  ref base A (101/115);  new base T (9/4);  total (111/120)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-01

Only 100 of 231 total aligned reads displayed.

AGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTC  >  NC_000913/2690902‑2691429
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                
aGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTCACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTAcca                                                                                                                                                                                                                                                               >  2:161231/1‑274 (MQ=255)
    gctggtgggTCATCGCTTGCCTGCCGCGGGCAACAATAATGGTAGTGAGGTGAACCACGGCCGCCTGATCGCGAAGCGCATTCAGGGTTTGCCAGGCTTGCCATGCAAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCt                                                                                                                                                                                                                                                           <  4:106851/268‑1 (MQ=255)
                     tCACTGGGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAgg                                                                                                                                                                                                                                           <  4:192190/274‑1 (MQ=255)
                                          aGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTAAGATACAACAACGCCGGGAGTTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTca                                                                                                                                                                                                                     >  2:153876/1‑274 (MQ=255)
                                                gTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTAGACGTTTCGCCCg                                                                                                                                                                                                                 >  2:57113/1‑273 (MQ=255)
                                                                 cgctggATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCAAGCACCAAAAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACg                                                                                                                                                                                               <  4:8780/269‑1 (MQ=255)
                                                                   ccTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGAAATTACCCGGATTTGTTATGGAAGAAGGCGGTGCCTAACTCGACGATTCACCCGAAGGGTGATATAGGTaagatg                                                                                                                                                                                            >  1:29182/1‑270 (MQ=255)
                                                                                             tGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCCGATATCAGAAGTTGGACGGCAGGCAcc                                                                                                                                                                  >  2:137391/1‑275 (MQ=255)
                                                                                                           tgccaGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACg                                                                                                                                                                                               <  2:64397/232‑1 (MQ=255)
                                                                                                                     aaCAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTAccc                                                                                                                                                                                                                                                               >  2:317386/1‑158 (MQ=255)
                                                                                                                                cagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCgg                                                                                                                                                                                                                                                             <  6:127009/149‑1 (MQ=255)
                                                                                                                                cagAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGTCGGCCGGCACCTTTTTGTGCGTCATTCGTAtt                                                                                                                                             >  4:58962/1‑261 (MQ=255)
                                                                                                                                  gAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCt                                                                                                                                                                                                                                                           >  8:926621/1‑149 (MQ=255)
                                                                                                                                     tCAGCAAAAATGCCAGCATTATCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTg                                                                                                                                                                                                                                                        <  7:1338310/149‑1 (MQ=255)
                                                                                                                                      cAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCTTTGt                                                                                                                                                                                                                                                       <  8:669592/149‑1 (MQ=255)
                                                                                                                                        gCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTa                                                                                                                                                                                                                                                     >  6:906886/1‑149 (MQ=255)
                                                                                                                                         cAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCCCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTAt                                                                                                                                                                                                                                                    <  8:690540/149‑1 (MQ=255)
                                                                                                                                          aaaaaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATg                                                                                                                                                                                                                                                   >  5:536393/1‑149 (MQ=255)
                                                                                                                                             aaTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCTTTGttgtatct                                                                                                                                                                                                                                                <  8:570945/149‑7 (MQ=39)
                                                                                                                                              aTGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAAt                                                                                                                                                                                                                                               <  6:724182/149‑1 (MQ=255)
                                                                                                                                                  cAGCATTACCAGTTGACGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTAc                                                                                                                                                                                                                                                                 <  7:611825/127‑1 (MQ=255)
                                                                                                                                                       ttACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTa                                                                                                                                                                                                                                                     <  4:160653/134‑1 (MQ=255)
                                                                                                                                                            agatGTCGTAATGACCGGGGAAAAACGGGCTAGCATTTCAAGGTGTTACTCTCGTCAGACGTGAATAGCCTGATGCTAACCGAGGGGAAGGTCAGATACCACAAAACCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTaa                                                                                                                                                                                                                                 <  8:899578/146‑1 (MQ=255)
                                                                                                                                                            aGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTaa                                                                                                                                                                                                                                 >  7:424587/1‑149 (MQ=255)
                                                                                                                                                            aGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTaa                                                                                                                                                                                                                                 >  8:154420/1‑149 (MQ=255)
                                                                                                                                                             gTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAAc                                                                                                                                                                                                                                >  6:964140/1‑149 (MQ=255)
                                                                                                                                                                  gggtagggggggggAAAAACGGGCCAGCGTTTCAAGGTGTGACTCTCGTCAGACGCGAAGAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCa                                                                                             <  4:301467/265‑1 (MQ=255)
                                                                                                                                                                  cGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGAc                                                                                                                                                                                                                           <  6:207882/149‑1 (MQ=255)
                                                                                                                                                                    tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTGAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCa                                                                                             <  2:299150/273‑1 (MQ=255)
                                                                                                                                                                    tAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAAc                                                                                                                                                                                                                                <  6:311326/142‑1 (MQ=255)
                                                                                                                                                                       tGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTc                                                                                                                                                                                                                      <  8:1013065/149‑1 (MQ=255)
                                                                                                                                                                          gCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGcc                                                                                                                                                                                                                   >  5:960126/1‑149 (MQ=255)
                                                                                                                                                                              ggAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGAt                                                                                                                                                                                                               <  7:926621/149‑1 (MQ=255)
                                                                                                                                                                               gAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATg                                                                                                                                                                                                              <  7:159801/149‑1 (MQ=255)
                                                                                                                                                                               gAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATg                                                                                                                                                                                                              >  7:574367/1‑149 (MQ=255)
                                                                                                                                                                                  aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAg                                                                                                                                                                                                                                            <  6:846071/116‑1 (MQ=255)
                                                                                                                                                                                  aaaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                           <  8:1129237/149‑1 (MQ=255)
                                                                                                                                                                                   aaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTg                                                                                                                                                                                                          <  5:112478/149‑1 (MQ=255)
                                                                                                                                                                                   aaCGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAAATCGACGTTTCGCCCGATGGTTg                                                                                                                                                                                                          >  7:826452/1‑149 (MQ=255)
                                                                                                                                                                                        gCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGTAATTCCCGGCttt                                                                                                                                                                                                                                                         >  6:627283/1‑97 (MQ=39)
                                                                                                                                                                                          cAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGc                                                                                                                                                                                                   <  5:858282/149‑1 (MQ=255)
                                                                                                                                                                                          cAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGc                                                                                                                                                                                                   <  5:413262/149‑1 (MQ=255)
                                                                                                                                                                                           aGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                           >  7:310243/1‑140 (MQ=255)
                                                                                                                                                                                           aGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                           <  8:310243/140‑1 (MQ=255)
                                                                                                                                                                                           aGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCt                                                                                                                                                                                                  >  7:876399/1‑149 (MQ=255)
                                                                                                                                                                                            gCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGTCTAACTCGACGTTTCGccc                                                                                                                                                                                                                  <  7:1150793/132‑1 (MQ=255)
                                                                                                                                                                                              gTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACg                                                                                                                                                                                               >  7:1103973/1‑149 (MQ=255)
                                                                                                                                                                                                      gtgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATc                                                                                                                                                                                       >  8:724750/1‑149 (MQ=255)
                                                                                                                                                                                                       tgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGAtata                                                                                                                                                                                                     <  7:856743/134‑1 (MQ=255)
                                                                                                                                                                                                       tgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCa                                                                                                                                                                                      >  5:295558/1‑149 (MQ=255)
                                                                                                                                                                                                       tgtTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCa                                                                                                                                                                                      <  8:5966/149‑1 (MQ=255)
                                                                                                                                                                                                            ctctCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGt                                                                                                                                                                                 >  8:435891/1‑149 (MQ=255)
                                                                                                                                                                                                            ctctCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGt                                                                                                                                                                                 <  7:1280696/149‑1 (MQ=255)
                                                                                                                                                                                                             tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAAc                                                                                                                                                                                                                                <  7:66998/101‑1 (MQ=255)
                                                                                                                                                                                                             tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGatat                                                                                                                                                                                                      <  6:1429407/127‑1 (MQ=255)
                                                                                                                                                                                                             tctcGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGtt                                                                                                                                                                                >  5:790924/1‑149 (MQ=255)
                                                                                                                                                                                                             tctcGTCAGACGCGAATAGCCTGACGCTAACAGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGTTGATATCAGAAGtt                                                                                                                                                                                >  7:393654/1‑149 (MQ=255)
                                                                                                                                                                                                                cGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGACATAGCTACGCTGATATCAGAAGTTGGa                                                                                                                                                                             <  6:536393/149‑1 (MQ=255)
                                                                                                                                                                                                                 gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATa                                                             <  4:249738/260‑1 (MQ=255)
                                                                                                                                                                                                                 gTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATa                                                             >  2:249738/1‑260 (MQ=255)
                                                                                                                                                                                                                    aGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggc                                                                                                                                                                         >  5:1399196/1‑149 (MQ=255)
                                                                                                                                                                                                                     gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGatat                                                                                                                                                                                        >  5:70233/1‑133 (MQ=255)
                                                                                                                                                                                                                     gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                                                                                                                                        <  7:724750/149‑1 (MQ=255)
                                                                                                                                                                                                                     gACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACggca                                                                                                                                                                        <  8:526617/149‑1 (MQ=255)
                                                                                                                                                                                                                             tAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCtt                                                                                                                                                                <  7:1452442/149‑1 (MQ=255)
                                                                                                                                                                                                                                ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCg                                                                                                                                                                                                                                         >  6:954659/1‑73 (MQ=255)
                                                                                                                                                                                                                                ccTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCttgtt                                                                                                                                                             <  7:266407/149‑1 (MQ=255)
                                                                                                                                                                                                                                  tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                                                                                                                                           >  7:210456/1‑149 (MQ=255)
                                                                                                                                                                                                                                  tGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                                                                                                                                           <  8:103301/149‑1 (MQ=255)
                                                                                                                                                                                                                                    aTGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGc                                                                                                                                                         >  5:975684/1‑149 (MQ=255)
                                                                                                                                                                                                                                     tGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCg                                                                                                                                                        >  5:1442277/1‑149 (MQ=255)
                                                                                                                                                                                                                                      gCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGt                                                                                                                                                       <  5:635666/149‑1 (MQ=255)
                                                                                                                                                                                                                                          aCCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCAtt                                                                                                                                                   <  6:683085/149‑1 (MQ=255)
                                                                                                                                                                                                                                           ccGAGGGGAAGTTCAGATTCAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTc                                                                                                                                                  <  7:219731/149‑1 (MQ=255)
                                                                                                                                                                                                                                           ccGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTg                    <  4:141514/275‑1 (MQ=255)
                                                                                                                                                                                                                                            cGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTtgttgt                                                                                                                                                           <  7:316390/139‑1 (MQ=255)
                                                                                                                                                                                                                                             gAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGtt                                                                                                                                                                                                           >  6:411614/1‑90 (MQ=255)
                                                                                                                                                                                                                                               ggggAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAt                                                                                                                                              <  8:166688/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 ggAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACGCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTAttt                                                                                                                                            <  8:826452/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     gTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGAc                                                                                                                                                                            >  7:448998/1‑113 (MQ=255)
                                                                                                                                                                                                                                                       tCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTa                                                                                                                                      >  5:1288305/1‑149 (MQ=255)
                                                                                                                                                                                                                                                        cAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAg                                                                                                                                     <  8:1346482/149‑1 (MQ=255)
                                                                                                                                                                                                                                                          gATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCa                                                                                                                                   <  8:1398655/149‑1 (MQ=255)
                                                                                                                                                                                                                                                              caacaaAGCCGGGAATTACCCGGGGTGGTTAAGGAAAAATGGGCGGCCTAACTAGACGTTTCGCCCGATGGTTTATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCTTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTctcctc  <  4:3310/274‑1 (MQ=255)
                                                                                                                                                                                                                                                               aacaaAGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGATACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                                                                                                                              >  8:25016/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aacaaAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTcc                                                                                                                              >  8:1280666/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                acaaAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTccc                                                                                                                             >  6:760118/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 caaAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCg                                                                                                                            <  6:437576/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaaGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGa                                                                                                                           >  5:36869/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aaaGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGa                                                                                                                           >  5:46112/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aaaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGa                                                                                                                           >  5:877049/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                  aaaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGa                                                                                                                           >  7:346498/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCAc                                                                                                                                  >  7:1360742/1‑141 (MQ=255)
                                                                                                                                                                                                                                                                   aaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa                                                                                                                          >  8:355396/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                   aaGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGaa                                                                                                                          <  8:210456/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aGCCGGGTAATTCCCGGCTTTGTTGTGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                                         >  5:937083/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aGCCGGGTAATTCCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                                         >  6:237905/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    aGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                                         <  8:280680/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                    aGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAg                                                                                                                         <  8:393508/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                            aaTTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGa                                                                                                                  <  5:1037705/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                
AGCGCCGCGTTGGTCATCGCTTCACTGCGCCGGGCATCAATAAGCGTAGTGCGGTTAACCAGCGCCGCCTGATCGCTAAGCGCATTCAGGCTTTGCCAGGCTTGCCATGCCAGCACCAACAGGGGCAGCAGAATCAGCAAAAATGCCAGCATTACCAGTTGTCGTAATGAGCGGGGAAAAACGGGCCAGCGTTTCAAGGTGTTACTCTCGTCAGACGCGAATAGCCTGATGCTAACCGAGGGGAAGTTCAGATACAACAAAGCCGGGAATTACCCGGCTTTGTTATGGAATAAGGCGGTGCCTAACTCGACGTTTCGCCCGATGGTTGATATAGCTACGCTGATATCAGAAGTTGGACGGCAGGCACCTTGTTGTGCGTCATTCGTATTTTATGTAGCACGTCCCGAAGGGGCTGACATAAGTCGGTGAATGAGCCACTGGTTACTATTATGCAGTAACTGTGCCAATAAAGAAAATAGTTTGGTAACGTACTGATTATACGTTGCTTTGAGGGGTTTATGTCTCCTC  >  NC_000913/2690902‑2691429

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: