Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,902,057 C→A 38.2% W3L (TGG→TTG)  sugE ← transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,902,0570CA38.2% 16.9 / 4.5 13W3L (TGG→TTG) sugEtransporter
Reads supporting (aligned to +/- strand):  ref base C (3/5);  new base A (3/2);  total (6/7)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCAGGCCCAGGCTGATGACCATGGCGCCGACGGTAAGGACAGTGGGTAGCGGTCGAGTGAAACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCT  >  NC_002947/1901927‑1902188
                                                                                                                                  |                                                                                                                                   
cagcaggcccagGCTGATGACCATGGCGCCGACGGTAAGGACAGTGGGTAGCGGTCGAGTGAAACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGAc                                                                                                                                <  1:470168/136‑1 (MQ=255)
          agGCTGATGACCATGGCGCCGACGGTAAGGACAGTGGGTAGCGGTCGAGTGACACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCAAGGACATGAgtgggg                                                                                                                      >  2:382797/1‑133 (MQ=255)
          agGCTGATGACCATGGCGCCGACGGTAAGGACAGTGGGTAGCGGTCGAGTGAAACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCAAGGACATGAgtgtgt                                                                                                                      >  2:380355/1‑136 (MQ=255)
             cTGATGACCATGGCGCCGACGGTAAGGACAGTGGGTAGCGGTCGAGTGAAACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCAAGGACATGAGTGTGtctc                                                                                                                   >  2:429348/1‑136 (MQ=255)
                   aCCATGGCGCCGACGGTAAGGACAGTGGGTAGCGGTCGAGTGAAACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGt                                                                                                             >  1:403394/1‑136 (MQ=255)
                              gACGGTAAGGACAGTGGGTAGCGGTCGAGTGAAACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATgg                                                                                                  <  2:311838/136‑1 (MQ=255)
                                                           tGAAACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGcg                                                                     >  1:259183/1‑136 (MQ=255)
                                                           tGAAACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGcg                                                                     >  1:64499/1‑136 (MQ=255)
                                                           tGAAACCGTCGGTGTATTTGAGGCCGACGGCACAGCCGACCTCGAACAGGCCGGCGAAAAACATGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGcg                                                                     <  1:480818/136‑1 (MQ=255)
                                                                                                         cAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGc                       <  2:233603/136‑1 (MQ=255)
                                                                                                                              gatCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCt  <  1:38423/136‑1 (MQ=255)
                                                                                                                              gatCAAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCt  <  1:380355/136‑1 (MQ=255)
                                                                                                                              gatCAAGGACATGAGTGTGTCTCCATCGTGTGAATGAGGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCt  <  1:382797/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                   
CAGCAGGCCCAGGCTGATGACCATGGCGCCGACGGTAAGGACAGTGGGTAGCGGTCGAGTGAAACCGTCGGTGTATTTGAGGCCGACGGCCCAGCCGACCTCGAACAGGCCGGCGAAAAACAGGATGATCCAGGACATGAGTGTGTCTCCATCGTTTGAATGATGGGGCCGTCCCCGGAGTGGTCACGCGGTGCGTCGTGAGGTCGTCCTCACCGTGGGCACCATGGTGCCTGTCTGCGGCCTTGCGGGCAAGACCGCGCCT  >  NC_002947/1901927‑1902188

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