Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 71,548 C→T G66G (GGC→GGT)  yabI → DedA family protein YabI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,5480CT83.3% 46.0 / ‑6.0 19G66G (GGC→GGT) yabIDedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (2/13);  minor base G (3/0);  total (6/13)
Fisher's exact test for biased strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGAT‑TTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAA  >  NC_000913/71365‑71766
                                                                                                                                                                                       |                                                                                                                                                                                                                           
ggAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTg                                                                                                                                                           <  1:166219/249‑1 (MQ=255)
ggAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGt                                                                                                                                                         <  1:44481/251‑1 (MQ=255)
              cccAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGGTTGGGAGCGGCGGGTTTAGGTTTTGGCCCGCCCGGCTGGCCGGGATTTTCGGGTTTTTGCTGGGCGAATGGGTCTCCCTTTTGGTGGGGTGGGGGTTTTAAAAA‑CCGTTGCCTCGCTGGTCAttttctaagaa                                                                                                                                            >  1:130313/1‑242 (MQ=255)
                             aTTCATTGGTGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCattact                                                                                                                             <  2:262906/250‑1 (MQ=255)
                                     aTGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAAc                                                                                                                    <  1:75519/251‑1 (MQ=255)
                                       ggCGGTGGTGGTGGTGGCCTTTCAGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGGGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGGTTCTGGCACGGCTGGCTGGCAGGGGTTGTCGGGTGGCTGCTGGGCGACAGGAT‑CTCCTTCTGGGTGGGCTGGGGTTTTTAAAAACCGTTGCATTGCTGGGCATTTTTGAAGAAAAAAAAAAGAATACTCGATaaaaatg                                                                                                                  >  1:211933/1‑248 (MQ=255)
                                           gtggtgTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGGCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACa                                                                                                              <  1:164887/251‑1 (MQ=255)
                                                  tggtggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGtt                                                                                                        <  2:89535/250‑1 (MQ=255)
                                                   ggtggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGAATCTACCAGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGCTTGGCAGCGGCGAGTTAAGTTTCTGGAACGCCTGGCTGGCAGGGGATGTCGGGTGTCTGCTGGGCGACTGGGT‑CTCCTTCTGGCTGGGCTGGCGGTTTAAAAAGCCGCTGCATCGACGGTCATTTCTGGAGAAAAAACAAGCCTTACTCGGTAAAACAGAACATGAGTTGc                                                                                                      >  1:44097/1‑251 (MQ=255)
                                                              tGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACGCGATAAAACTGAACATGCGTTGCATCAACACAg                                                                                            <  2:106320/250‑1 (MQ=255)
                                                                          cgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCAtt                                                                               <  1:184893/251‑1 (MQ=255)
                                                                                  ggTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGGTGGGGGCGCTGATTGGCAGCGGCGAGTTTAGTTTCTGGGACGCCTGGCTGGCCGGGGTTGTCGGTTGGCTGCTGGGCGGCCGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTTAAAAGACGTTGCATCGCTGGGCATTTTTGTAGAAAAAACAAGCATTACTCGCTAAAACTGAACCTGCGTTGCATCAACACAGCCTGTTCACCATTCTtgtcgg                                                                       >  2:111874/1‑251 (MQ=255)
                                                                                      tGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGgtcggtcg                                                                    <  1:132457/250‑1 (MQ=255)
                                                                                        gttCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGTATTGTCGGTTGTCTGCTGGTCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAGACTGAACATGCGTTGCATCAACACAGCATGTTCTCCATTCTGgtcggtcgctt                                                                 <  2:80310/250‑4 (MQ=255)
                                                                                                                            gcgcTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGc                             <  1:23455/251‑1 (MQ=255)
                                                                                                                                 gATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGTCGGGaa                        <  1:188347/251‑1 (MQ=255)
                                                                                                                                 gATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGaa                        <  2:90662/251‑1 (MQ=255)
                                                                                                                                           ggcgAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCggggaagttggggct              >  2:2690/1‑239 (MQ=255)
                                                                                                                                                        tCTGGCACGCCTGGCTGGCAGGGATTGTCGGATGTCTGCTGGGCGACTGGAT‑CTCCTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGCCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTCGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGACAATGGAGGCGGGAATGCTGGATCTGCCGGTGGcaaa  >  1:158514/1‑247 (MQ=255)
                                                                                                                                                                                       |                                                                                                                                                                                                                           
GGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGAT‑TTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAA  >  NC_000913/71365‑71766

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