Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,510,609 C→T P327L (CCG→CTG)  yfeO → putative transport protein YfeO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,510,6090CT86.7% 31.0 / NA 17P327L (CCG→CTG) yfeOputative transport protein YfeO
Reads supporting (aligned to +/- strand):  ref base C (2/0);  major base T (11/2);  minor base G (0/2);  total (13/4)
Fisher's exact test for biased strand distribution p-value = 5.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGG‑TAATTAAAC‑TTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGAT  >  NC_000913/2510373‑2510843
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
cgcgATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCCGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGG‑GAATTAAAC‑TTGCCCCCCTGGTGGTTTCTGCCGCCCGTGGCTTTTGCGGTGGGCGCATCTTCCCGGGCGGGTTTGTCg                                                                                                                                                                                                                              >  1:51454/1‑251 (MQ=255)
     tgatgCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGAATTCTGGGGGTTATTGGTGGCACAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGGGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGGGAAATTAAA‑‑TTGCCGCCCCTGTCGTTGGTTCCGCCAGTGGGTTTCGCGGGGGGCGAATCTTCCTTGCCTTGTTTGTCgggggg                                                                                                                                                                                                                         >  2:227463/1‑247 (MQ=255)
                     aaaTCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTATCGGGATTCTGGGGGTTATTGGTGGACCAGTTTCGATGTTTAAAGGGCTGGATGAGGTGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCCGCGATTAATTTTTGCTGGCGG‑TAATTAAAA‑CTGACCCCCTGGTCGGTTCTGCCGGCCGTGGGTTTCGCGGTGGGCGTATCTTCCTGGCAGGGTTTTTCGCCGTGGCATTTggggttatgg                                                                                                                                                                                                         >  2:273892/1‑244 (MQ=255)
                      aaTCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGG‑TAATTAAAC‑TTGCCGCCCCGGGCGTTGCTGCCCCCCGTGGCTTTCGCGGTGGGCGGATCTTCCCGGCAGTGTTTTTCGGCGTGGCCTTAGGGTTTAtgct                                                                                                                                                                                                        >  1:267272/1‑251 (MQ=255)
                              gcgcgggcggggTATTGGCGGATTTATTCGCGGTATTCTGGGGGGTATTGGTGGACCAGTGTCGGTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGGGGCAAATCAGGCTGTCAGCTCCAGCGATTACTTTTGGCTGTCGG‑TAATTAAAC‑TTGCCGCATTGGTCGTTCCTGCGGCCAGTGGCTTGCGCGGTGGGCGAATCTTCCGGGCAGTGTTTGTCGGCGTGGCATAAGGGTTGATGCTGCATGAg                                                                                                                                                                                                 <  2:312890/241‑1 (MQ=255)
                                 cgtgtggggTATTGGCGTATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTTCTTTTTGCTGGCGG‑TAATTAAAC‑TTGCCGCCTTGGTCGTGGCGGCCGCCAGGGGCTATCGCGGTTGGCGAAGCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGATGATGCTGCATGAGCAc                                                                                                                                                                                              <  2:310047/244‑1 (MQ=255)
                                                                  tcgggggtTTAGTGGTGGACCTGTTGCGCTTTTTAAGGGGCGGGAGGAGATGCACCAGAGGGGGGCAAATCAGGTTTCCAGCACCAGGGTTTTCTTTTTGCGGCCGG‑TAATTAACC‑TTCCCGCCTTGGTCGTTGCTGCCGCCAGTGTCTTGCGGGGGGGGCAAATCTTCCTGGCGTTGTTTTTCGGCGTGGCATTAGGTTTAATGCGGAACGGCCCCGTTCCCGCCGTACCAGCGGCAATAACCGTTTCt                                                                                                                                                             <  1:236675/247‑1 (MQ=255)
                                                                                                                             tggtggCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGG‑TAATTAAAC‑TTGCCGCCCTGGGCGTTGCTGCCGCCAGTGGCTTTCGCGGGGGGCGGATCTTCCTGGCAGTGTTTGTTGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTTCCAGCCGCCATAACCGGTTCTTGCGCTTTTCTCGGCATCGTGCGGGTGGGAACACGCGATGGCGGGTTAAGTCTTTTTAtt                                                                                                 >  2:437399/1‑250 (MQ=255)
                                                                                                                                tggCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGG‑TAATTAAAC‑TTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGGGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGGGGGAACACGCGATGGCTGGTTAAGTCTTTTTATggcg                                                                                              >  1:21733/1‑251 (MQ=255)
                                                                                                                                            cTTTCAGCACCAGCGATTACTTTTTGCTGGCGG‑TAATTAAAC‑TTGCCGCACCGGTCGTTGCTGCCGCAAGTGGCTTTCGCGGTGGGCGAATATTCCTGGCAGTCTTTGTCGGCGTGGCATTTGGGTTGATGCTGCCTGAGCACGTTCCCCCCGTACCAGCGGCAATAACAGTTTCTTCGGGTATTCTCGGCATCGTGCTGGTGGTAACAAGCGATGGCTGGTTAAGTGCTTTTATGGAGGCAGGCGtttgt                                                                                  >  1:349975/1‑248 (MQ=255)
                                                                                                                                              ttcaccaccggggATTTCTTTTGGCTGGTG‑‑TAAATAAACCTTGCCGCCTGGGTGGTTGCTGCCGCAGGTGTTTTGCGCGGGGGGCAACTCTCCCGGGCAGTTTTTGCCGGCGTGGCATTGGGGTTGATGCTGCATGACCCCTTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTa                                                                                  <  1:576692/237‑1 (MQ=255)
                                                                                                                                                                          cGG‑TAATTAAAC‑TTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGACGTCGTTGGACCCAATACCAAATTGCTACCgcacctcttt                                                    >  1:504214/1‑243 (MQ=255)
                                                                                                                                                                                 ttAAAC‑TTGCCGCCCTGGTCGTTGCTGCCGCCCGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGATGTTCACACTACCACAATGCTACCgccgattcttattggc                                              >  2:519697/1‑237 (MQ=255)
                                                                                                                                                                                               ccTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTTGTTGTACCCCATACCCAATTGCTCCCGCCGctcagtgtttacatggtttctgcat                                 >  2:596316/1‑229 (MQ=255)
                                                                                                                                                                                               ccTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACCTTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCAt                                 >  1:245581/1‑251 (MQ=255)
                                                                                                                                                                                                                   gTGGCTTTCGCGGTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGttatttggaggga               >  2:202252/1‑241 (MQ=255)
                                                                                                                                                                                                                              ggTGGGCGAATCTTCCTGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCCCATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCgatgatgat  >  1:213572/1‑251 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                          
CGCGATGATGCATCAAATGAAAAATCCGGTGCTCGTGCTGGGTATTGGCGGATTTATTCTCGGTATTCTGGGGGTTATTGGTGGACCAGTTTCGCTGTTTAAAGGGCTGGATGAGATGCAGCAGATGGTGGCAAATCAGGCTTTCAGCACCAGCGATTACTTTTTGCTGGCGG‑TAATTAAAC‑TTGCCGCCCTGGTCGTTGCTGCCGCCAGTGGCTTTCGCGGTGGGCGAATCTTCCCGGCAGTGTTTGTCGGCGTGGCATTAGGGTTGATGCTGCATGAGCACGTTCCCGCCGTACCAGCGGCAATAACCGTTTCTTGCGCTATTCTCGGCATCGTGCTGGTGGTAACACGCGATGGCTGGTTAAGTCTTTTTATGGCGGCAGTCGTTGTACCCAATACCACATTGCTACCGCTGCTCTGTATCGTCATGCTTCCGGCATGGCTGTTATTAGCAGGTAAGCCGATGATGAT  >  NC_000913/2510373‑2510843

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