Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 518,196 G→A F28F (TTC→TTT cnoX ← chaperedoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913518,1960GA90.9% 21.7 / NA 11F28F (TTC→TTTcnoXchaperedoxin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (7/3);  minor base C (1/0);  total (8/3)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

CGGCCCCTGGAAGCCATCTACCGGTTGCCCGTTCTGGAACAGATACACGGTCGGAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGAT  >  NC_000913/517992‑518425
                                                                                                                                                                                                            |                                                                                                                                                                                                                                     
cGGCCCCTGGAAGCCATCTACCGGTTGCCCGTTCTGGAACAGATACACGGTCGGAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAAAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGt                                                                                                                                                                                         <  1:134706/251‑1 (MQ=255)
                 cTACCGGTTGCCCGTTCTGGAACAGCTACACGGTCGTCATCGCACGCAGACCAAACTGCGCGGCAATCATCTGATCCGAGTAGCAGTCCAGCTTCGACAGAATAAATTGCCCGTTGTACTGCGCCGCGCGGCTTTCCAGAATTGGGGGTAACTGCAAACAGTGCTGGCTCCGTTCAGACCAAAAATAAACCAGCACCGGAGTGGTCATCTACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTg                                                                                                                                                                        >  1:84429/1‑251 (MQ=255)
                                                            angcaGACCAACCTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAAAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGat                                                                                                                              >  1:132156/1‑250 (MQ=255)
                                                                                                   ncAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGCATTGGGGTTAACTGCAAAAAGTGCTGGCTACGTTCAGACCAAAAATAAACCAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTg                                                                                                                                                                        >  1:41479/2‑169 (MQ=255)
                                                                                                      ntCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGATGCTTTCCCGAATTGGGGTTAACTGCAAACAGTACTGGCTACGTTCAGACCAAAAATAAAACAGCACCGGCGTGGTCATCGACTGTTCCAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGTGGGCTTAAGCACGCGCTTCAACTCAcc                                                                                   >  1:101324/2‑251 (MQ=255)
                                                                                                                             gCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTCCGTTCAGACCAAAAATAAAACAGCACCGGAGTGGTCATCGCCTGTTCCAGAACTTGTTGCAGGTTAGATTCGTTAATGTTGACAATCTTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATCTTTTGTCCCTCAcg                                                            >  1:21775/1‑251 (MQ=255)
                                                                                                                               ccGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGTTACGTTGAGACCAACAATACAACAGCAACGGAGTGGTCATCGCCTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCAcgcg                                                          >  1:81284/1‑251 (MQ=255)
                                                                                                                                            gcgcgAGGCTTTCTAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAAAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAgg                                             >  1:132633/2‑251 (MQ=255)
                                                                                                                                                gAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAAAACAGCACCGGAGTGGTCATCGACTGTGCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGGTTAAGCGCGCGCTTCAACTCCCCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGgcgat                                         <  1:156208/251‑3 (MQ=255)
                                                                                                                                                    cTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAAAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAgcgcg                                                                                                 <  1:84071/191‑1 (MQ=255)
                                                                                                                                                                                       ncTACGTTCAGACCAAAAATAAAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGAt  >  1:100698/2‑251 (MQ=255)
                                                                                                                                                                                                            |                                                                                                                                                                                                                                     
CGGCCCCTGGAAGCCATCTACCGGTTGCCCGTTCTGGAACAGATACACGGTCGGAATCGCACGCAGACCAAACTGCGCGGCAATCATCTGCTCCGCGTCGCAGTCCAGCTTCGCCAGAATAAATTGCCCGTTGTACTGCGCCGCGAGGCTTTCCAGAATTGGGGTTAACTGCAAACAGTGCTGGCTACGTTCAGACCAAAAATAGAACAGCACCGGAGTGGTCATCGACTGTTCAAGAACCTGTTGCAGGTTAGATTCGTTAATGTTGACAATATTTTCTACGGACATGGAGTCGCTCTCTGTTGTCGATATTTTCTTTGACATGGGGGCTTAAGCGCGCGCTTCAACTCACCCCTGCAATATTTTGTCCATCACGCGCCCCGGCAGCAGGCGCTTAAGCACCATTACCGCCCAGGTCACCAGCGTCACCGGAT  >  NC_000913/517992‑518425

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