Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,133,301 +G coding (657/714 nt) yghS ← putative ATP‑binding protein YghS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,133,2971.G83.3% 11.8 / ‑2.3 6coding (661/714 nt)yghSputative ATP‑binding protein YghS
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (0/5);  total (0/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 = 1.00e+00

TGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATC  >  NC_000913/3133166‑3133423
                                                                                                                                    |                                                                                                                               
tGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTAACGCTTGCAGAATTTCATcggcg                                                                                                                                <  1:416130/133‑1 (MQ=255)
tGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCAGTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATcggcg                                                                                                                                <  1:298414/133‑1 (MQ=255)
    gcgcGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCggggg                                                                                                                            <  1:294271/133‑2 (MQ=255)
               cTATTAACAGGTTGATTTTGTAGTCCATCCATTCGTTAAATTATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGTGTCGCCCATc                                                                                                                 <  2:154272/133‑1 (MQ=255)
                                 tgTAGTGCATCCATTCGTTAAATCCTCCACTACCCATTTTGCTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGTGTCGCCCATCCAGTTCAAGGATCTTTgc                                                                                               <  2:139773/133‑1 (MQ=255)
                                      tGCATCCATTCGTTACATCATCCCCTATCCATATTGGTGATTAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCTCCAt                                                                                          <  2:158713/133‑1 (MQ=255)
                                         aTCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGa                                                                                       <  1:510884/133‑1 (MQ=255)
                                                 gTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGt                                                                               <  1:8209/133‑1 (MQ=255)
                                                                                                             gACGCTTGCAGAATTTCATCGGCGGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCa                   <  1:440401/133‑1 (MQ=255)
                                                                                                                              aTCGGCGGGGATGTCCCCCATCCATTTCTAGGGTCTTTGCGACATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTtcatc  <  2:196881/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                                               
TGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATGTCAGTTGCGGCGTAACGGCGATTTTTTCCTGTAACGCTGCCAGTTGGTGGTCAGGTTTACGCGCAAACGCGGTTTGTTCATC  >  NC_000913/3133166‑3133423

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