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→ACCyehPVMA 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

GCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCG  >  NC_000913/2208760‑2209211
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  
gCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTa                                                                                                                                                                                                                           >  2:17564/1‑239 (MQ=255)
           agagTTATTAATGAAAGAACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGGGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATa                                                                                                                                                                                                                >  2:60999/1‑239 (MQ=255)
                      tgcaaGTACAGTTGGGCGGCGGGATCAATATCGCCAGTGCCGGGGAGTATGGTCGGCAACTTATTGAACCACCAGCGAAAATCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGCCGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGgcg                                                                                                                                                                                                     <  2:837352/236‑1 (MQ=255)
                         aaGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGGGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGGCGGGAATAAAAGTGCTGGGCTGGGCACCGCACGTAGGCACCGCAACCCCTTGATATAACCGCGAAACGAGCCAGgcgcgg                                                                                                                                                                                                  >  1:953112/1‑237 (MQ=255)
                                    ggcggcGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCATAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGAACCGCAACACCTTGCTATGACCGCGATACGGCCCAAGCGCTGGTTAATGttg                                                                                                                                                                                        >  2:88210/1‑236 (MQ=255)
                                       ggcggGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGgcg                                                                                                                                                                                    <  2:328843/238‑1 (MQ=255)
                                           ggACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGGGCTGGGACTGGCAGCGCTCGATAGCCCCGCAACACCTTGCTATGACCGCGATACGGCCCCGGCGCGGGTGAAGGTCGgcgccaaa                                                                                                                                                                               >  1:1018113/1‑235 (MQ=255)
                                                ctcctcccccgTGCCGTGGGTTTTGGTAGGAAATTTATTGAAAAATACGTGAAAAGCGTAATTCTCCTGTGGAGGGATTTTTAGGAGGGGGGTTCACCATCTTTCCTGCCGCATCGGGTGAAAAAGTGTGTCCAGGGCGGCATCAAAGTGCTGGGCCTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAgc                                                                                                                                                                          <  1:904575/229‑1 (MQ=255)
                                                      gCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGGCGCCAAATAgccgc                                                                                                                                                                       >  2:192744/1‑236 (MQ=255)
                                                      gCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCAt                                                                                                                                                                    >  2:297257/1‑239 (MQ=255)
                                                            gCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCCCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCGGGTTAATGTCGGCGCACAAAAACACCCCATGACGcc                                                                                                                                                              >  1:712730/1‑239 (MQ=255)
                                                               ggggAGTATGGGCGGCAACTTATTGAACAACCAGCAAAAAGCGTCATTATCCGCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGGGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCggg                                                                                                                                                           <  2:234578/237‑1 (MQ=255)
                                                                tGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACCCCTTGCTAGTACCGCGATACGGCCCAGGCGCTGGTTAAGTTCCGGGCCCAAAAAGCCGCCAAGACGCCgggg                                                                                                                                                          >  1:5050/1‑238 (MQ=255)
                                                                     tATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCCCAACACCTTGCTATGACCCCGATACGGCCCAGGCGCGGTGTAAGGTCGGCGCAACAATACCCGCCacacccccgggccggcg                                                                                                                                                     >  1:1164/1‑223 (MQ=255)
                                                                     gttggcgggaaattttttaaacaccCAGCAAAAGGTGTAATTACCCGGGGGGGGGATTTTCAGGAAGGGGTTTCATCATCATTACTGACGCATCGGGTGAAAAAGGGTGCCCAGGGGGGCATAAAAGTGCGGGGACTGGCAGCGCTCGATAGCACCGCAACCCCTTGCTAGGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGc                                                                                                                                                      <  1:297257/219‑1 (MQ=255)
                                                                            ggCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCCGGCGCTGGTTAATGTCGGCGCCCAAATAGCCGCCATGACGCCGGGCGAGCTGGcacaa                                                                                                                                              >  2:251371/1‑236 (MQ=255)
                                                                             gCAACTTCTTGAACAACCAGCGAAACGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCACCGCTCGCTAGCACTGCAAAACCTTGTGAGGACCGCGAGACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAAAAGCCCCCCTGACGCCGGGCGGGCGGGcatcat                                                                                                                                             >  1:290272/1‑239 (MQ=255)
                                                                             gCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGcatcat                                                                                                                                             >  2:741323/1‑239 (MQ=255)
                                                                             gCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCCCAACACCTTGCTATGACCGCGCTACGGCCCCGTCCCTGGTTAATGTCGGCCCCCAAATAACCGCCCACACGCCGGGGGAggcggacaccg                                                                                                                                             >  1:532896/1‑229 (MQ=255)
                                                                                    aTTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCATCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATTTCGGCGCACAAATAGCGGCCATGACGCCGGGCGAGCTGGCATCATGGCttac                                                                                                                                      >  2:879507/1‑237 (MQ=255)
                                                                                            acCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGTCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTTGGGAGAAtt                                                                                                                              >  2:966025/1‑238 (MQ=255)
                                                                                                         acgttttccctcttgggggttttttACGAGGGGGGTTCAGCATCATTCCTGCAGACTCAGTGGAAAAGGTGTGCCCAGGGCGGCATCAATGTGCTGGGCCTGGCAGCGCTCGATAGAACCGCAACACCTTGCTAGGACCGCGATACGGCCCGGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGa                                                                                                                   <  1:339555/217‑1 (MQ=255)
                                                                                                                        aGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAATCTTCAGCCATGAATTACCTACggcggtc                                                                                                   >  2:339555/1‑232 (MQ=255)
                                                                                                                             tttttGCGAGGGGGGTTCATCCTCATTCGTGAGGAACCGAGGAAAAATTGGGGCCCAGGGGGGCTTCAAAGGGCGGGGACTTGCAGGGCTAGATAGCACCGCACCACCTGGCTATGACCGCGATACGGCCCAGGTGTTGGTTAATTTCGGCGCACAAATGGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAattatt                                                                                             <  1:741323/239‑1 (MQ=255)
                                                                                                                             ttttgAGGAGGGGGTTTCCTCATCATTATGGCCGATCCGGGGGAAAAAGGTTGTCCGGGGCGTCTACAAAGCGCTGGGACTGGCAGCCTTGGATGCCACCGAACTCCGTTGCTATGACCGCGAGACGGCCCAGGCTCGGGTTAAGGTAAGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAattatt                                                                                             <  1:947280/239‑1 (MQ=255)
                                                                                                                                       gggggTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGGATTCACTTCGTCCGGGATTATTTGAACTGac                                                                                    >  2:705838/1‑238 (MQ=255)
                                                                                                                                       gggggTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGACTTATTAGAACTGac                                                                                    >  2:928223/1‑238 (MQ=255)
                                                                                                                                                    atttctggccgcagagggTAAAAAGGTTTGTCCGGGGCGGCAAAAAAGTGCTGGGCTTGGCAGCGCTCGATAGCACCGCAACCACTTCCTAGGACCGCGATACGGCCCAGGCGGTGGTTAATGTCGGCCCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGa                                                                      <  1:17564/223‑1 (MQ=255)
                                                                                                                                                         cTGGCGCATCAGGTGAAAAAGTGTGTCCAGGGCGGCATCAAAGTGCTGGGACGGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATCCGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCg                                                                 <  2:1164/239‑1 (MQ=255)
                                                                                                                                                         cTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACCCCGCAGGCCCTGACGGCg                                                                 >  2:370509/1‑239 (MQ=255)
                                                                                                                                                          ctgaggcctaGGGGAAAAAGGGTGTCCCAGGGGGGCAAAAAGGTGCTGGGTGTGGAAGGCTCCAAAGACACCAAAAACCCTGGTGATGCCCGGAATGCGCCCCGGGCGCTGGTTAATGCCGGGCCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCg                                                                 <  1:705838/230‑1 (MQ=255)
                                                                                                                                                             cGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGCATTCACTCCGCCCGCAATTATTAGAAATGACACCGCAGGCCCAGGCGGGGTTaa                                                             >  1:970823/1‑239 (MQ=255)
                                                                                                                                                                ataaGGTGAAAAAGTGTGCCCAGGGCGGCATAAAAGTGCTGGGACGGGAAGGGCTCGATAGCCCCGCAACCCCTTGCTAGGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGa                                                                      <  2:142868/224‑1 (MQ=255)
                                                                                                                                                                 tCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCCCTACGTCCGGAATTATTAAAAATGACACCGCAGGCCCTGCCGGCGTTAAGCaa                                                        >  2:7053/1‑239 (MQ=255)
                                                                                                                                                                      taaaaacGTGTGTCCAGAGCGCAATCAAAGTGCTGGGACTGTCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGCCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGc                                                    <  2:953112/235‑1 (MQ=255)
                                                                                                                                                                         aaaaGGGTGTCAAGGGGGGCATCAAGGTGCTGGGACTGGCAGGCCTCGATAGCCCCGCACCCCCTTGCTATGACCGGGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGtt                                              <  2:5050/239‑1 (MQ=255)
                                                                                                                                                                                gtgtCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCCTGAATTCACTACGTCCGGAATTATTAGAACTGACACCGACAGCCCTTCAGGCGTTAAGGCATGCCGGGTTTTGta                                         >  2:369973/1‑238 (MQ=255)
                                                                                                                                                                                   tCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAg                                                                                                                         >  2:116638/1‑158 (MQ=255)
                                                                                                                                                                                     tgtgtaaaaggaacaggTGCGGGGACGGGCAGCGCTCGATAGCCCCGAAACCCCTTGCTATGCCCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAgcg                                    <  2:315591/223‑1 (MQ=255)
                                                                                                                                                                                     ccAGAGCGGACTAAAAGGGCGGGGACTGGCAGCGCTCGATAGCCCCGCAACCCCTTGCTATGCCCGCGATCCGGCCCAGGCGCTGGTTAATGTCGGCGCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAgcg                                    <  2:1018113/239‑1 (MQ=255)
                                                                                                                                                                                      cAGGGCGGCATCAAACTGGTGGGATTGGCACGGCTCGGTAGCACCGCAACCCATTGATATGACCGCGATACGGCCCAGGCGCTGGTTAATGGCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAgcg                                    <  2:290272/238‑1 (MQ=255)
                                                                                                                                                                                         aggtggttaaaaGAGCCGGGACTGGCAGCGCTCAATAGCACCGAAACCCTTTGCTATGACCGAGATGCGGCCCAGGCGCTGGTTAATGTCGGCGCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAg                                 <  2:990092/229‑1 (MQ=255)
                                                                                                                                                                                          gCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATACACGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAAACATCCGCCATTGGTTACATACGCCCGGAATTGTTTGGACTTGCACCCCCAGCCCCGCAGGCATTAAAAAAAGCCCGGTTTTTTTaaa                                      >  1:385565/1‑232 (MQ=255)
                                                                                                                                                                                                        cTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGCAATCTTCAGTCCTGAATTCACTAAGTCCCGAATTATTTTAACTGAACCCCCAGGCCCTAACGCCGTTAAGACATGCCGGTTTTTTTAAGGCCAGCCTTAAGGAAATGGaa                 >  1:243220/1‑239 (MQ=255)
                                                                                                                                                                                                          gggACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCa                                                         >  2:338006/1‑199 (MQ=255)
                                                                                                                                                                                                                ggcggcGCTCGATAGCCCCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCaa         <  2:712730/239‑1 (MQ=255)
                                                                                                                                                                                                                ggCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAATCTTCAGTCCTGAATTCACTACGTCCCGAATTATTAGAACTGACACCCCCGGGCCTGACGGCCTTAAACAATGCCGGGTTTTTTAAGCGCAGGCCTaaaggtggtctttaaaaaaa         >  1:315591/1‑221 (MQ=255)
                                                                                                                                                                                                                ggCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAATCTTCAGCCATGAATTCACAACGTCCGGAATTATTTAAACTGTAACCCCCAGCCCCGGTGCGGGTAATCAATACCGGGTTTTTTTAAGGCAGGCTTAAAGAACAGGGAAAGGGGaa         >  1:963665/1‑239 (MQ=255)
                                                                                                                                                                                                                 ggcaggcccgaaaagcccggCAACCCCTTGCTAGGACCGCGAACCGGCCCGGGCGCTGGTTAATGTCGGCGCAAAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCa          <  2:243220/219‑1 (MQ=255)
                                                                                                                                                                                                                  cAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGCATTATTTGAACTGGCACCCCCGGCCCTGACGGAGTTAAGCAATGCCGGGTTTTTTAAAACCCGCATTAAGAAACGGGaaaaaaggaaag       >  1:55479/1‑231 (MQ=255)
                                                                                                                                                                                                                  cAGCGCTCGATAGAACCGCAACACCTTGATATGACCGCGATACGGCACAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATTAAGCCGGGCGAGCTGGCATCATGGCTTTCGGCGAATCTTCAGTCCTGAATTCACTCCGCCCGGAAATCTTAGAACCGAAACCCCCGGCCCTGACGGCGCTAACGCAAGCCCGTTTTGTTTAGGGCGCACTCTAGGAACAGGAAAAAAGCaaag       >  1:190628/1‑237 (MQ=255)
                                                                                                                                                                                                                     cgcTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGAAATGCCGGTTTTTTTAAGTGCCGCCTTAAGGAACTGGAAAATGGCAACGTc     >  2:906291/1‑238 (MQ=255)
                                                                                                                                                                                                                      gcTCGATAGCACCGCAACCCCTTGCTATGCCCGCCATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTcc    <  2:545872/238‑1 (MQ=255)
                                                                                                                                                                                                                       cTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGGGAATCTTCAGGCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGAGCAGTCTTAAGGAACTGGAAAATGGGCACGTGCCg  >  2:413215/1‑239 (MQ=255)
                                                                                                                                                                                                                       cTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCCTGAATTCACTACGTCCGGAATTTTTAGAACTGACACCCCAGGCCCCGCCGGCGTTAAGCAATTCCGGGTTTGTTAACGGCCGTCTTAAGGAACTGGCACAAGGCACCGCCCCg  >  2:782581/1‑239 (MQ=255)
                                                                                                                                                                                                                                     cAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAg                                                           >  2:103939/1‑171 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                  
GCCGATCCGGTAGAGTTATTAATGAAAGTACAGTTGGGCGGCGGGACCAATATCGCCAGTGCCGTGGAGTATGGTCGGCAACTTATTGAACAACCAGCGAAAAGCGTCATTATCCTCGTGAGCGATTTTTACGAAGGGGGTTCATCATCATTACTGACGCATCAGGTGAAAAAGTGTGTCCAGAGCGGCATCAAAGTGCTGGGACTGGCAGCGCTCGATAGCACCGCAACACCTTGCTATGACCGCGATACGGCCCAGGCGCTGGTTAATGTCGGCGCACAAATAGCCGCCATGACGCCGGGCGAGCTGGCATCATGGCTTGCGGAGAATCTTCAGTCATGAATTCACTACGTCCGGAATTATTAGAACTGACACCGCAGGCCCTGACGGCGTTAAGCAATGCCGGTTTTGTTAAGCGCAGTCTTAAGGAACTGGAAAATGGCAACGTCCCG  >  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: