Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,208,990 A→C 20.1% T342T (ACA→ACC)  yehP → VMA domain putative YehL ATPase stimulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,208,9900AC20.1% 114.1 / 3.0 56T342T (ACA→ACC) yehPVMA domain putative YehL ATPase stimulator
Reads supporting (aligned to +/- strand):  ref base A (33/9);  new base C (2/11);  total (35/21)
Fisher's exact test for biased strand distribution p-value = 7.20e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCG  >  NC_000913/2208760‑2209211
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  
gCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTa                                                                                                                                                                                                                           >  2:17564/1‑239 (MQ=255)
           agagTTATTAATGAAAGAACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGGGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATa                                                                                                                                                                                                                >  2:60999/1‑239 (MQ=255)
                      tgcaaGTACAGTTGGGCGGCGGGATCAATATCGCCAGTGCCGGGGAGTATGGTCGGCAACTTATTGAACCACCAGCGAAAATCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGCCGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGgcg                                                                                                                                                                                                     <  2:837352/236‑1 (MQ=255)
                         aaGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGGGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGGCGGGAATAAAAGTGCTGGGCTGGGCACCGCACGTAGG‑CACCGCAACCCC‑TTGATATAACCGCGAAACGAGCCAGgcgcgg                                                                                                                                                                                                  >  1:953112/1‑237 (MQ=255)
                                    ggcggcGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCATAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑AACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAAGCGCTGGTTAATGttg                                                                                                                                                                                        >  2:88210/1‑236 (MQ=255)
                                       ggcggGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑Ggcg                                                                                                                                                                                    <  2:328843/238‑1 (MQ=255)
                                           ggACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGGGCTGGGACTGGCAGCGCTCGATAG‑CCCCGCAACACC‑TTGCTATGACCGCGATACGGCCCCGGCGCGGGTGAAGGTC‑Ggcgccaaa                                                                                                                                                                               >  1:1018113/1‑235 (MQ=255)
                                                ctcctcccccgTGCCGTGGGTTTTGGTAGGAAATTTATTGAAAAATACGTGAAAAGCGTAATTCTCCTGTGGAGGGATTTTTAGGAGGGGGGTTCACCATCTTTCCTGCCGCATCGGGTGAAAAAGT‑GTGT‑CCAGGGCGGCATCAAAGTGCTGGGCCTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAgc                                                                                                                                                                          <  1:904575/229‑1 (MQ=255)
                                                      gCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGGCGC‑CAAATAgccgc                                                                                                                                                                       >  2:192744/1‑236 (MQ=255)
                                                      gCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCAt                                                                                                                                                                    >  2:297257/1‑239 (MQ=255)
                                                            gCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCCCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCGGGTTAATGTC‑GGCGCACAAAAACACCCCATGACGcc                                                                                                                                                              >  1:712730/1‑239 (MQ=255)
                                                               ggggAGTATGGGCGGCAACTTATTGAACAACCAGCAAAAAGCGTCATTATCCGCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGGGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCggg                                                                                                                                                           <  2:234578/237‑1 (MQ=255)
                                                                tGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACCCC‑TTGCTAGTACCGCGATACGGCCCAGGCGCTGGTTAAGTTC‑CGGGCCCAAAAAGCCGCCAAGACGCCgggg                                                                                                                                                          >  1:5050/1‑238 (MQ=255)
                                                                     tATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCCCAACACC‑TTGCTATGACCCCGATACGGCCCAGGCGCGGTGTAAGGTC‑GGCGCAACAATACCCGCCacacccccgggccggcg                                                                                                                                                     >  1:1164/1‑223 (MQ=255)
                                                                     gttggcgggaaattttttaaacaccCAGCAAAAGGTGTAATTACCCGGGGGGGGGATTTTCAGGAAGGGGTTTCATCATCATTACTGACGCATCGGGTGAAAAAGG‑GTGC‑CCAGGGGGGCATAAAAGTGCGGGGACTGGCAGCGCTCGATAG‑CACCGCAACCCC‑TTGCTAGGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                                                                                                                                                      <  1:297257/219‑1 (MQ=255)
                                                                            ggCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCCGGCGCTGGTTAATGTC‑GGCGCCCAAATAGCCGCCATGACGCCGGGCGAGCTGGcacaa                                                                                                                                              >  2:251371/1‑236 (MQ=255)
                                                                             gCAACTTCTTGAACAACCAGCGAAACGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCACCGCTCGCTAG‑CACTGCAAAACC‑TTGTGAGGACCGCGAGACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAAAAGCCCCCCTGACGCCGGGCGGGCGGGcatcat                                                                                                                                             >  1:290272/1‑239 (MQ=255)
                                                                             gCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatcat                                                                                                                                             >  2:741323/1‑239 (MQ=255)
                                                                             gCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCCCAACACC‑TTGCTATGACCGCGCTACGGCCCCGTCCCTGGTTAATGTC‑GGCCCCCAAATAACCGCCCACACGCCGGGGGAggcggacaccg                                                                                                                                             >  1:532896/1‑229 (MQ=255)
                                                                                    aTTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCATCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATTTC‑GGCGCACAAATAGCGGCCATGACGCCGGGCGAGCTGGCATCATGGCttac                                                                                                                                      >  2:879507/1‑237 (MQ=255)
                                                                                            acCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GTCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTTGGGAGAAtt                                                                                                                              >  2:966025/1‑238 (MQ=255)
                                                                                                         acgttttccctcttgggggttttttACGAGGGGGGTTCAGCATCATTCCTGCAGACTCAGTGGAAAAGGT‑GTGC‑CCAGGGCGGCATCAATGTGCTGGGCCTGGCAGCGCTCGATAG‑AACCGCAACACC‑TTGCTAGGACCGCGATACGGCCCGGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGa                                                                                                                   <  1:339555/217‑1 (MQ=255)
                                                                                                                        aGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAATCTTCAGCCATGAATTACCTACggcggtc                                                                                                   >  2:339555/1‑232 (MQ=255)
                                                                                                                             tttttGCGAGGGGGGTTCATCCTCATTCGTGA‑GGAACCGAGGAAAAATTGGGGC‑CCAGGGGGGCTTCAAAGGGCGGGGACTTGCAGGGCTAGATAG‑CACCGCACCACC‑TGGCTATGACCGCGATACGGCCCAGGTGTTGGTTAATTTC‑GGCGCACAAATGGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAattatt                                                                                             <  1:741323/239‑1 (MQ=255)
                                                                                                                             ttttgAGGAGGGGGTTTCCTCATCATTATGGCCGATCCGGGGGAAAAAGG‑TTGT‑CCGGGGCGTCTACAAAGCGCTGGGACTGGCAGCCTTGGATGC‑CACCG‑AACTCCGTTGCTATGACCGCGAGACGGCCCAGGCTCGGGTTAAGGTA‑AGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAattatt                                                                                             <  1:947280/239‑1 (MQ=255)
                                                                                                                                       gggggTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGGATTCACTTCGTCCGGGATTATTTGAACTGac                                                                                    >  2:705838/1‑238 (MQ=255)
                                                                                                                                       gggggTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGACTTATTAGAACTGac                                                                                    >  2:928223/1‑238 (MQ=255)
                                                                                                                                                    atttctggccgcagagggTAAAAAGGT‑TTGT‑CCGGGGCGGCAAAAAAGTGCTGGGCTTGGCAGCGCTCGATAG‑CACCGCAACCAC‑TTCCTAGGACCGCGATACGGCCCAGGCGGTGGTTAATGTC‑GGCCCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGa                                                                      <  1:17564/223‑1 (MQ=255)
                                                                                                                                                         cTGGCGCATCAGGTGAAAAAGT‑GTGT‑CCAGGGCGGCATCAAAGTGCTGGGACGGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATCCGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCg                                                                 <  2:1164/239‑1 (MQ=255)
                                                                                                                                                         cTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACCCCGCAGGCCCTGACGGCg                                                                 >  2:370509/1‑239 (MQ=255)
                                                                                                                                                          ctgaggcctaGGGGAAAAAGG‑GTGTCCCAGGGGGGCAAAAAGGTGCTGGGTGTGGAAG‑GCTCCAAAGACACCAAAAACCC‑TGGTGATGCCCGGAATGCGCCCCGGGCGCTGGTTAATGCC‑GGGCCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCg                                                                 <  1:705838/230‑1 (MQ=255)
                                                                                                                                                             cGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGCATTCACTCCGCCCGCAATTATTAGAAATGACACCGCAGGCCCAGGCGGGGTTaa                                                             >  1:970823/1‑239 (MQ=255)
                                                                                                                                                                ataaGGTGAAAAAGT‑GTGC‑CCAGGGCGGCATAAAAGTGCTGGGACGGGAAGGGCTCGATAG‑CCCCGCAACCCC‑TTGCTAGGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGa                                                                      <  2:142868/224‑1 (MQ=255)
                                                                                                                                                                 tCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCCCTACGTCCGGAATTATTAAAAATGACACCGCAGGCCCTGCCGGCGTTAA‑GCaa                                                        >  2:7053/1‑239 (MQ=255)
                                                                                                                                                                      taaaaacGT‑GTGT‑CCAGAGCGCAATCAAAGTGCTGGGACTGTCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGCC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑c                                                    <  2:953112/235‑1 (MQ=255)
                                                                                                                                                                         aaaaGG‑GTGT‑CAAGGGGGGCATCAAGGTGCTGGGACTGGCAGGCCTCGATAG‑CCCCGCACCCCC‑TTGCTATGACCGGGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGtt                                              <  2:5050/239‑1 (MQ=255)
                                                                                                                                                                                gtgt‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCCTGAATTCACTACGTCCGGAATTATTAGAACTGACACCGACAGCCCTTCAGGCGTTAA‑GGCATG‑CCGGGTTTTGta                                         >  2:369973/1‑238 (MQ=255)
                                                                                                                                                                                   t‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAg                                                                                                                         >  2:116638/1‑158 (MQ=255)
                                                                                                                                                                                     tgtgtaaaaggaacaggTGCGGGGACGGGCAGCGCTCGATAG‑CCCCGAAACCCC‑TTGCTATGCCCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAgcg                                    <  2:315591/223‑1 (MQ=255)
                                                                                                                                                                                     ccAGAGCGGACTAAAAGGGCGGGGACTGGCAGCGCTCGATAG‑CCCCGCAACCCC‑TTGCTATGCCCGCGATCCGGCCCAGGCGCTGGTTAATGTC‑GGCGCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAgcg                                    <  2:1018113/239‑1 (MQ=255)
                                                                                                                                                                                      cAGGGCGGCATCAAACTGGTGGGATTGGCACGGCTCGGTAG‑CACCGCAACCCA‑TTGATATGACCGCGATACGGCCCAGGCGCTGGTTAATGGC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAgcg                                    <  2:290272/238‑1 (MQ=255)
                                                                                                                                                                                         aggtggttaaaaGAGCCGGGACTGGCAGCGCTCAATAG‑CACCGAAACCCT‑TTGCTATGACCGAGATGCGGCCCAGGCGCTGGTTAATGTC‑GGCGCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAGCGCAg                                 <  2:990092/229‑1 (MQ=255)
                                                                                                                                                                                          gCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATACACGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAAACATCCGCCATTGGTTACATACGCCCGGAATTGTTTGGACTTGCACCCCCAGCCCCGCAGGCATTAA‑AAAAAGCCC‑GGTTTTTTTaaa                                      >  1:385565/1‑232 (MQ=255)
                                                                                                                                                                                                        cTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGCAATCTTCAGTCCTGAATTCACTAAGTCCCGAATTATTTTAACTGAACCCCCAGGCCCTAACGCCGTTAA‑GACATG‑CC‑GGTTTTTTTAAGGCCAGCCTTAAGGAAATGGaa                 >  1:243220/1‑239 (MQ=255)
                                                                                                                                                                                                          gggACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCa                                                         >  2:338006/1‑199 (MQ=255)
                                                                                                                                                                                                                ggcggcGCTCGATAG‑CCCCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCaa         <  2:712730/239‑1 (MQ=255)
                                                                                                                                                                                                                ggCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAATCTTCAGTCCTGAATTCACTACGTCCCGAATTATTAGAACTGACACCCCCGGGCCTGACGGCCTTAA‑ACAATG‑CC‑GGGTTTTTTAAGCGCAGGCCTaaaggtggtctttaaaaaaa         >  1:315591/1‑221 (MQ=255)
                                                                                                                                                                                                                ggCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAATCTTCAGCCATGAATTCACAACGTCCGGAATTATTTAAACTGTAACCCCCAGCCCCGGTGCGGGTAA‑TCAATA‑CCGGGTTTTTTTAAG‑GCAGGCTTAAAGAACAGGGAAAGGGGaa         >  1:963665/1‑239 (MQ=255)
                                                                                                                                                                                                                 ggcaggcccgaaaagcccggCAACCCC‑TTGCTAGGACCGCGAACCGGCCCGGGCGCTGGTTAATGTC‑GGCGCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCa          <  2:243220/219‑1 (MQ=255)
                                                                                                                                                                                                                  cAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGCATTATTTGAACTGGCACCCCCGGCCCTGACGGAGTTAA‑GCAATG‑CC‑GGGTTTTTTAAAACCCGCATTAAGAAACGGGaaaaaaggaaag       >  1:55479/1‑231 (MQ=255)
                                                                                                                                                                                                                  cAGCGCTCGATAG‑AACCGCAACACC‑TTGATATGACCGCGATACGGCACAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATTAAGCCGGGCGAGCTGGCATCATGGCTTTCGGCGAATCTTCAGTCCTGAATTCACTCCGCCCGGAAATCTTAGAACCGAAACCCCCGGCCCTGACGGCGCTAACGCAA‑G‑CC‑CGTTTTGTTTAGGGCGCACTCTAGGAACAGGAAAAAAGCaaag       >  1:190628/1‑237 (MQ=255)
                                                                                                                                                                                                                     cgcTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GAAATG‑CC‑GGTTTTTTTAAGTGCCGCCTTAAGGAACTGGAAAATGGCAACGTc     >  2:906291/1‑238 (MQ=255)
                                                                                                                                                                                                                      gcTCGATAG‑CACCGCAACCCC‑TTGCTATGCCCGCCATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTcc    <  2:545872/238‑1 (MQ=255)
                                                                                                                                                                                                                       cTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAATCTTCAGGCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAGAGCAGTCTTAAGGAACTGGAAAATGGGCACGTGCCg  >  2:413215/1‑239 (MQ=255)
                                                                                                                                                                                                                       cTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCCTGAATTCACTACGTCCGGAATTTTTAGAACTGACACCCCAGGCCCCGCCGGCGTTAA‑GCAATT‑CC‑GGGTTTGTTAACGGCCGTCTTAAGGAACTGGCACAAGGCACCGCCCCg  >  2:782581/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑g                                                           >  2:103939/1‑171 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  
GCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGT‑GTGT‑CCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAG‑CACCGCAACACC‑TTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTC‑GGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAA‑GCAATG‑CC‑GGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCG  >  NC_000913/2208760‑2209211

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