Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,444,286 C→T 26.7% A296V (GCC→GTC)  yjfF → inner membrane putative ABC superfamily sugar transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,444,2860CT26.7% 13.5 / 5.3 15A296V (GCC→GTC) yjfFinner membrane putative ABC superfamily sugar transporter permease
Reads supporting (aligned to +/- strand):  ref base C (7/4);  new base T (2/2);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

AGTGGTAATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAG  >  CP009273/4444149‑4444423
                                                                                                                                         |                                                                                                                                         
ccgatctaTTGGCGGTACGCTTTTGAGCGGTGGCGTTGGGACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGATGAGTTCCTGGTGGACGAAAATCGTCATCg                                                                                                                                      <  2:234791/136‑1 (MQ=255)
ccgatctaTTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGTCATCg                                                                                                                                      <  2:302056/136‑1 (MQ=255)
    gTAATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTTTTGTTGCTATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGtgt                                                                                                                                  >  2:162504/1‑141 (MQ=255)
    gTAATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAGATCGCCAGCGGTATTttgtt                                                                                                                            >  1:14935/1‑149 (MQ=255)
        ttGGCGGTACGCTTTTGAGCGGTGGCGTTGGGACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGTCATCggagatcggaagagc                                                                                                                        >  1:262827/1‑136 (MQ=255)
        ttGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGTCATCggagatcggaagggc                                                                                                                        >  1:335101/1‑136 (MQ=255)
                    ttttGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGAGGAAAATCGCAATCGGTATTTTGTTGATTAttttt                                                                                                                  >  2:331706/1‑143 (MQ=255)
                                                          tttGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCa                                                                      <  1:343600/149‑1 (MQ=255)
                                                                gtggcgATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGCTTAACGGTGCTGTGGGAGAAGCGTCAGAGTTc                                                                >  1:220555/1‑149 (MQ=255)
                                                                            ggACTGATTCAGACTTACATAACCTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATGACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGt                                                          <  2:133951/143‑1 (MQ=255)
                                                                                                aaaCTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCCTTACAGCGTGGATTAACGGTGCTGTGCGAGAATCGCCTGAGTTCGCTACTGACAAGAGTCAAAATTgctgcccgat                                >  1:160833/1‑141 (MQ=255)
                                                                                                     ttGACGGCACGCTGAGTTCCTGGTGGACGAAAACCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATaaaaa                           <  1:19028/149‑1 (MQ=255)
                                                                                                                 tGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTTTAGCATTTCAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAGGAGTAAACATTGCGCACCAATGAAAACGCctc                     >  2:163751/1‑143 (MQ=255)
                                                                                                                          ggtggACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAAAAAAAACGCCTCCCCGCGTgg            >  2:122616/1‑143 (MQ=255)
                                                                                                                              gACGCAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGGGGATAATTGGTGCCTTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAg  <  1:177159/149‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                         
AGTGGTAATTGGCGGTACGCTTTTGAGCGGTGGCGTTGGAACGGTATTAGGGACGCTTTTTGGCGTGGCGATTCAGGGACTGATTCAGACTTACATAAACTTTGATGGCACGCTGAGTTCCTGGTGGACGAAAATCGCCATCGGTATTTTGTTGTTTATTTTTATAGCATTACAGCGTGGATTAACGGTGCTGTGGGAGAATCGTCAGAGTTCGCCAGTGACAAGAGTCAACATTGCGCAGCAATAAAAACGCCTCTCCGTGTGGAGAGGCGCAG  >  CP009273/4444149‑4444423

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