Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,465,154 C→T L25F (CTC→TTC)  gspM → Type II secretion system protein GspM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,465,1540CT92.3% 27.0 / ‑6.1 13L25F (CTC→TTC) gspMType II secretion system protein GspM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (5/7);  minor base A (1/0);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

GGAGGCGAACTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCA  >  NC_000913/3464917‑3465358
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                
ggAGGCGAACTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGTTCGGCGtttttt                                                                                                                                                                                                      <  2:93703/250‑1 (MQ=255)
                cAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGCTCTCTCGTTCGGCGTTTTTTTTTTGCTTgggggggt                                                                                                                                                                                     >  1:205025/1‑246 (MQ=255)
                                                    aaGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGCTCTCTCGTTCTGGGTTTTTTGTTGGGTTGGGGGGATAAAAACGAAAGATACCTACATTGCAGAACATCAggcc                                                                                                                                                >  1:100912/1‑248 (MQ=255)
                                                           aaGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGTTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACgcgc                                                                                                                                         <  1:144212/251‑1 (MQ=255)
                                                                       atcgTTTCTTTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGTTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTaaga                                                                                                                           <  2:117267/251‑1 (MQ=255)
                                                                                      tAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGTTCGTTGCGGCACTGGCGGATCTCTCCATCCGCGTTTTTTTTTTGGTTTGGGGGGTTAAAACCATAAATCACTACCTTGCACGGAATTAGTTCCCCCCGCAAAAAATTAAAAAGGaaattcaaa                                                                                                               >  2:170915/1‑240 (MQ=255)
                                                                                                                aaaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCCTTCGGCGGTTTTTTTTTGCTTGGGGGGGTAAAACCCATAGATACCTACCTTTCCGAGCATCCGGCCCCCCCGCACAAAATTTAAAAGGGACTCCAAATGAAGCAGGATCCAGCGAGCcccccaggt                                                                                  >  1:44339/1‑246 (MQ=255)
                                                                                                                 aaaCAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGTTCGGGGTTTTTTTTTTGGTTGGGGGGGTAAAAACGACAGATAACTCCATTGCAGAGCATCAGTCCCACGCGCCAAAAAATAAAAAAGGCATCAAAAGGATGCAGGATCAAGCGAGAACCCCTgggt                                                                                 >  1:87323/1‑249 (MQ=255)
                                                                                                                                                       aGGAACTCATGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGTTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAat                                           <  1:221796/251‑1 (MQ=255)
                                                                                                                                                                                aTGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGTTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTg                 <  1:235192/251‑1 (MQ=255)
                                                                                                                                                                                             aaaTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGTTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGTATGCAGGAt                                                                                                    <  2:199414/156‑1 (MQ=255)
                                                                                                                                                                                               aTCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGTTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGGCATCCTGCGCTTACGCAGCCCATCAAAAAGAGCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCa  <  2:152633/251‑1 (MQ=255)
                                                                                                                                                                                                tCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGGTCTCTCCTTCGGCGTTTTTTTTTGGGTTTGGGGGGTAAAACCGATAGATACCTACATTGTAGAGCATCAGGCCCCCGCGCAAAAAAATAAGAAGGAAATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCCTCAAAAATATCCTTCTCGAGGGAAGAAAAAcagagaaatttgttatcac  >  2:145440/1‑233 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                
GGAGGCGAACTATAGCCAGAAGCAAAAGTCGCTGACCTTAAAAATGAGTGCTAAGAGTGAAGCTAATATCGATCGTTTCTGTGAGTTAACACAGTCCTGGCTTCCGATGGAAAAAACAGAAAAAGATCCGGTCAGTGGTGTATGGACAGTAAGGAACTCAGGGAAATGATTAAATCATGGTGGGCAGAAAAATCCACATCGGAAAAACAGATCGTTGCGGCACTGGCGGTTCTCTCGCTCGGCGTTTTTTGTTGGCTTGGGGTGATAAAACCGATAGATACCTACATTGCAGAGCATCAGTCCCACGCGCAAAAAATTAAGAAGGACATCAAATGGATGCAGGATCAAGCGAGCACCCATGGTTTACTCGGTCATCCTGCGCTTACGCAGCCCATCAAAAATATCCTTCTCGAGGAAGCAAAACGTGAGAATTTGGCTATCA  >  NC_000913/3464917‑3465358

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