Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,330,931 G→A E144K (GAA→AAA)  proP → osmolyte:H(+) symporter ProP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,330,9310GA90.0% 19.6 / ‑3.5 10E144K (GAA→AAA) proPosmolyte:H(+) symporter ProP
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (7/2);  total (7/3)
Fisher's exact test for biased strand distribution p-value = 3.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CACTTTCTCCGTTCCCTTTCTGATTCGACCGCTTGGCGGACTCTTCTTTGGTATGTTGGGCGATAAATATGGTCGCCAGAAGATCCTCGCTATCACTATTGTGATTATGTCGATCAGTACGTTCTGTATTGGCTTAATACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCGTGGTGGTGTTAATTTCGACCATTGTCGGCGAAGCGAACTTCCTCGATTGGGGCTGGCGTATTCCGTTCTTTATCGCTCTGC  >  NC_000913/4330711‑4331126
                                                                                                                                                                                                                            |                                                                                                                                                                                                    
naCTTTCTCCGTTCCCTTTCTGATTCGACCGCTTGGCGGACTCTTCTTTGGTATGTTGGGCAATAAATATGGTCGCCAGAAGATCCTCACTATCACTATTGTGATTATGTCGATCAGTACGTTCTGTATTGGCTTAATACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTAAATATACCGGGGCGTCGATATTTGTTGCgg                                                                                                                                                                       >  1:261303/2‑251 (MQ=255)
    ttCTCCGTTCCCTTTCTGATTCGTCCGCTTGGCGCGCTCTTCTTTGGTATGTTGGGCAATAAATAAGGTCGCCAGAAGAACCTCAATATAAATAATGAGATTATGACGATCAGAACGTTCTGTATTGGCTTAATACCGTCCTACGACACGATTGGTATTTGGGAACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTAGGCGGTAAATATACCGGGGCGGCGATATTTGTTGCGGAAt                                                                                                                                                                    >  1:129501/1‑250 (MQ=255)
                          gACCGCTGGGCGGACTCGTCTTTGGTATGTTGGGCAATAAATATGGTCGCGAGAAGATCCTCTCTATCACGATTGTGATTATGTCGATCAGGACGTTCTGTATTGGCTTAATACCGTCCTACGACACGATTGGTATTTGGGGCCCGATTCTGCGGTTTATCTGTAAGGTGTCACGAGGTTTCGCGGTCGGCGTGTAATATAGCGGGGCGTGGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGc                                                                                                                                             <  1:215001/251‑1 (MQ=255)
                                                                      ggTCGCCAGAAGATCCTCACTATCACTATTGTGATTATGTCGATCAGTACATTCTGTATTGGCTTAATACCGTCCTACGACACGATTGGTATTTGGGCAACGATCCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGTCGGTAACTATACCGGGGCGGCGATATTTGTTGCGGAATACTCCTCTGACCGTAAACGTGGCTTTATTGGCAGCTGGCTGGACTTcagagcggaataaaccaagtc                                                                                                 >  1:228668/1‑231 (MQ=255)
                                                                                                              ngATCAGTACGTTCTGTATTGGCTTAATACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTAAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATg                                                                                                                                       >  1:85620/2‑173 (MQ=255)
                                                                                                               gATCAGTACGTTCTGTATTGGCTTAATACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTAAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCGTGGTGGTGTTAATTTCAACCATTGTCg                                                        >  1:311178/1‑251 (MQ=255)
                                                                                                                                                                attggggCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTAAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCGGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGTTGCGGGCGTGGTGGTGTTAATTTCGACCATTGTCGGCGAAGCGAACTTCCTCGATTg                                  <  1:237950/220‑1 (MQ=255)
                                                                                                                                                                     ggCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTAAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTGGTGCTGGGTGCGGGCGTGGTGGTGTTAATTTCGACCATTGTCGGCGAAGCGAACTTCCTCGATTGGGGCTGGCGTATTCCGTTCTTTATCGCTCTGc  <  1:156557/251‑1 (MQ=255)
                                                                                                                                                                         ccGATTCTGCTGTTGATCTGTAAGATGGCACAAGGGTTCTCGGTCGGCGGTAAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCACTGACCGTAAACGTGGCTTTATGGGCCgctggct                                                                                                                            >  1:297285/1‑125 (MQ=255)
                                                                                                                                                                                    gTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTAAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCt                                                                                              >  1:209337/1‑144 (MQ=255)
                                                                                                                                                                                                                            |                                                                                                                                                                                                    
CACTTTCTCCGTTCCCTTTCTGATTCGACCGCTTGGCGGACTCTTCTTTGGTATGTTGGGCGATAAATATGGTCGCCAGAAGATCCTCGCTATCACTATTGTGATTATGTCGATCAGTACGTTCTGTATTGGCTTAATACCGTCCTACGACACGATTGGTATTTGGGCACCGATTCTGCTGTTGATCTGTAAGATGGCACAAGGTTTCTCGGTCGGCGGTGAATATACCGGGGCGTCGATATTTGTTGCGGAATACTCCCCTGACCGTAAACGTGGCTTTATGGGCAGCTGGCTGGACTTCGGTTCTATTGCCGGGTTTGTGCTGGGTGCGGGCGTGGTGGTGTTAATTTCGACCATTGTCGGCGAAGCGAACTTCCTCGATTGGGGCTGGCGTATTCCGTTCTTTATCGCTCTGC  >  NC_000913/4330711‑4331126

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