Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,895,120 G→A D157N (GAT→AAT)  yieF → quinone reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,895,1200GA81.8% 13.1 / 0.0 11D157N (GAT→AAT) yieFquinone reductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (4/5);  total (5/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 = 4.41e-01

ACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAAC  >  NC_000913/3894905‑3895308
                                                                                                                                                                                                                       |                                                                                                                                                                                            
actCGGTACCGGGTGGTGTTAAAAATGGCATCGACGGGTGTTCCCGCCTGCCGGATCAGCCGCTGGCCGGTAAGCCGGTATTGATTCAGACCAGCTCAATAGGCGTGATTGGCGGCGCGCGCTGTCAGTATGTGCTGCGCCAGATGCTTGTTTTCCTCGATGCAATGGATATGAACAGGCCTGAATTTATGTGGGGTGTGATTCAGAACAAAGTTAATCCGCAAGCCGGAGAAGTGATTGATCAGGGTACg                                                                                                                                                           <  1:139518/251‑1 (MQ=255)
                gCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCATTGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTAATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCg                                                                                                                                           >  1:322344/1‑251 (MQ=255)
                                                                                                                         ttctcattttCACCTGCGCCAGATTTTGGTTTTCCTCGGTGCACTGGTGATGAACAACACGTAATTTATGGGCGCGGTGTTCCAGAACAAATTTAACCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGAATTTGGTGATTTTATTCAGCTAGTTAAGATCTAAATAAAAACCCCGCCAGCCATCATGCATGTCCGGTTTTTAACGCGCTATCGATTTTAg                                   <  1:26857/241‑1 (MQ=255)
                                                                                                                                                          ncTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTAATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGAGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAAc  >  1:59692/2‑250 (MQ=255)
                                                                                                                                                                       gTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCCGAACAAAGTTGATACGCAAACCGGCGAAGTGATTGATCAgg                                                                                                                                                                >  1:128949/1‑79 (MQ=255)
                                                                                                                                                                               cAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTAATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCg                                                                                                                                           <  1:260040/92‑1 (MQ=255)
                                                                                                                                                                                    cNGANTTTATNGGCGNCGTGATTCAGAACAANGTTGATCCGCAAACCGGATNNGNGATTGNTCAgg                                                                                                                                                                <  2:128949/66‑1 (MQ=255)
                                                                                                                                                                                                     gTGATTCAGAACAAAGTTAATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAAttn                                                                                                                                   <  1:40339/78‑2 (MQ=255)
                                                                                                                                                                                                            aGAACAAAGTTAATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGc                                                                                                          >  1:133862/1‑96 (MQ=255)
                                                                                                                                                                                                             naaCAAAGTTAATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGAt                                                                                                >  1:205362/2‑105 (MQ=255)
                                                                                                                                                                                                             gAACAAAGTTAATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGAt                                                                                                <  2:205362/105‑1 (MQ=255)
                                                                                                                                                                                                                       |                                                                                                                                                                                            
ACTCGGTACCGGGTGGGCTGAAAAATGCCATCGACTGGCTTTCCCGCCTGCCGGATCAACCGCTGGCCGGTAAACCGGTATTGATTCAGACCAGCTCAATGGGCGTGATTGGCGGCGCGCGCTGTCAGTATCACCTGCGCCAGATTCTGGTTTTCCTCGATGCAATGGTGATGAACAAGCCGGAATTTATGGGCGGCGTGATTCAGAACAAAGTTGATCCGCAAACCGGAGAAGTGATTGATCAGGGTACGCTGGACCACCTGACCGGGCAATTGACCGCATTTGGTGAGTTTATTCAGCGAGTTAAGATCTAAATAAAAAACCCGCCAGCAATCATGCATGGCGGGTTTTTAACGCGCTATCGATTTTAGTGAGCGTCGATAAAGACAATCTTCAGGATAAAC  >  NC_000913/3894905‑3895308

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