Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 623,386 G→A G363S (GGC→AGC)  entS → enterobactin exporter EntS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913623,3860GA85.7% 13.8 / NA 7G363S (GGC→AGC) entSenterobactin exporter EntS
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (2/4);  minor base C (1/0);  total (3/4)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACGCTGGGATCGTTCCTCGCCATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCAT  >  NC_000913/623180‑623589
                                                                                                                                                                                                              |                                                                                                                                                                                                           
ccacggtgggatctggacTCGCCATTGTGCTGTTGGGCCTGATGCCGATGTGGATTTGAGTCGTGGGTTGTCTTGCGCGGTTCGTCTGGTTGAGTGCGGTCAGCTCGGTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATTGTAGGGCGGATTAACGGTTTGTGGACGGCTCAGAACGTGACGGGCGATTCCATAAGCGCGGCGATGCTGGGTGGTTTGGGCGCGATGATGACACCGGtt                                                                                                                                                                 <  1:123966/234‑1 (MQ=255)
         gATCGTTCCTCGCCATTGGTCTGTTTGGACTGATGCCGATGTGGATTTTCGGCGGAGTTAGCCTGGCGCTGTTCGGCTGGGTGAGTGCGGTCAGATCGTTGATGCAATACACAATGCTGCAAACGCAAACCCCGGAAGAGATGTTAGGGCGGATTAAAGGATTGTGTACGGAGCAGAACGAGACGGGCGATGCCATACGCTCGGCGATGCTGGGGGGTTTGGGCGAGATGATGACACCGGGTGCTTCCGc                                                                                                                                                         >  1:123875/1‑250 (MQ=255)
                                     ccTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCCATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGACCGTGACGGGCGATGCCATAAGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGAtt                                                                                                                             >  1:68304/1‑250 (MQ=255)
                                           gCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAAGCGCGGCGCTGCTGGGTGGGTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGc                                                                                                                       <  1:64447/250‑1 (MQ=255)
                                                                                                          gTTGCTGCAATACACATTGCGGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAAGCGCGGCGCTGCGGGGGGGTTTGGGCGCGATGATGACACCGTTGGCTGCCGCAAGCGGGGGCGGTTTTGGTTTGTTGATTATCGGGGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAg                                                       <  1:155998/251‑1 (MQ=255)
                                                                                                                                                               ggATTAACGGTTTGTGGACGGCGCAGAACGTGGCGGGCGATGCCATAAGCGCGGCTCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGGGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCAt  <  1:68119/251‑1 (MQ=255)
                                                                                                                                                                               gACGGCGCAGAACGTGACGGGCGATGCCATAAGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTTTTGCTGNGGCnnnnnggnnnnnnnnnn                                                                                   >  2:99792/1‑143 (MQ=255)
                                                                                                                                                                                                              |                                                                                                                                                                                                           
CCACGCTGGGATCGTTCCTCGCCATTGGTCTGTTTGGCCTGATGCCGATGTGGATTTTAGGCGTGGTTTGTCTGGCGCTGTTCGGCTGGTTGAGTGCGGTCAGCTCGTTGCTGCAATACACAATGCTGCAAACGCAAACCCCGGAAGCGATGTTAGGGCGGATTAACGGTTTGTGGACGGCGCAGAACGTGACGGGCGATGCCATAGGCGCGGCGCTGCTGGGTGGTTTGGGCGCGATGATGACACCGGTTGCTTCCGCAAGCGCGAGCGGTTTTGGTTTGTTGATTATCGGCGTGTTGTTATTGCTGGTGCTGGTGGAGTTGCGACATTTTCGCCAGACGCCGCCGCAGGTGACAGCGTCCGACAGTTAATGCTTAAAACAGCGCCTTAAGCCTATCCAGCACTTGCAT  >  NC_000913/623180‑623589

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