Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 71,551 2 bp→TC coding (201‑202/765 nt) yabI → DedA family protein YabI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,5510CT89.5% 50.0 / NA 19C67C (TGC→TGTyabIDedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (4/13);  minor base G (2/0);  total (6/13)
Fisher's exact test for biased strand distribution p-value = 8.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_00091371,5520TC94.7% 51.8 / ‑5.7 19L68L (TTG→CTG) yabIDedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (5/13);  total (6/13)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

GGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAA  >  NC_000913/71365‑71766
                                                                                                                                                                                          ||                                                                                                                                                                                                                       
ggAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTg                                                                                                                                                           <  1:166219/249‑1 (MQ=255)
ggAACACTTTATTACCCAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGt                                                                                                                                                         <  1:44481/251‑1 (MQ=255)
              cccAATCTACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGGTTGGGAGCGGCGGGTTTAGGTTTTGGCCCGCCCGGCTGGCCGGGATTTTCGGGTTTTTGCTGGGCGAATGGGTCTCCCTTTTGGTGGGGTGGGGGTTTTAAAAACCGTTGCCTCGCTGGTCAttttctaagaa                                                                                                                                            >  1:130313/1‑242 (MQ=255)
                             aTTCATTGGTGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCattact                                                                                                                             <  2:262906/250‑1 (MQ=255)
                                     aTGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAAc                                                                                                                    <  1:75519/251‑1 (MQ=255)
                                       ggCGGTGGTGGTGGTGGCCTTTCAGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGGGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGGTTCTGGCACGGCTGGCTGGCAGGGGTTGTCGGGTGGCTGCTGGGCGACAGGATCTCCTTCTGGGTGGGCTGGGGTTTTTAAAAACCGTTGCATTGCTGGGCATTTTTGAAGAAAAAAAAAAGAATACTCGATaaaaatg                                                                                                                  >  1:211933/1‑248 (MQ=255)
                                           gtggtgTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGGCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACa                                                                                                              <  1:164887/251‑1 (MQ=255)
                                                  tggtggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGtt                                                                                                        <  2:89535/250‑1 (MQ=255)
                                                   ggtggCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGAATCTACCAGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGCTTGGCAGCGGCGAGTTAAGTTTCTGGAACGCCTGGCTGGCAGGGGATGTCGGGTGTCTGCTGGGCGACTGGGTCTCCTTCTGGCTGGGCTGGCGGTTTAAAAAGCCGCTGCATCGACGGTCATTTCTGGAGAAAAAACAAGCCTTACTCGGTAAAACAGAACATGAGTTGc                                                                                                      >  1:44097/1‑251 (MQ=255)
                                                              tGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACGCGATAAAACTGAACATGCGTTGCATCAACACAg                                                                                            <  2:106320/250‑1 (MQ=255)
                                                                          cgctggTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCAtt                                                                               <  1:184893/251‑1 (MQ=255)
                                                                                  ggTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGGTGGGGGCGCTGATTGGCAGCGGCGAGTTTAGTTTCTGGGACGCCTGGCTGGCCGGGGTTGTCGGTTGGCTGCTGGGCGGCCGGATCTCCTTCTGGCTGGGCTGGCGTTTTTAAAAGACGTTGCATCGCTGGGCATTTTTGTAGAAAAAACAAGCATTACTCGCTAAAACTGAACCTGCGTTGCATCAACACAGCCTGTTCACCATTCTtgtcgg                                                                       >  2:111874/1‑251 (MQ=255)
                                                                                      tGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGgtcggtcg                                                                    <  1:132457/250‑1 (MQ=255)
                                                                                        gttCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGTATTGTCGGTTGTCTGCTGGTCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAGACTGAACATGCGTTGCATCAACACAGCATGTTCTCCATTCTGgtcggtcgctt                                                                 <  2:80310/250‑4 (MQ=255)
                                                                                                                            gcgcTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGc                             <  1:23455/251‑1 (MQ=255)
                                                                                                                                 gATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGTCGGGaa                        <  1:188347/251‑1 (MQ=255)
                                                                                                                                 gATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGaa                        <  2:90662/251‑1 (MQ=255)
                                                                                                                                           ggcgAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTGTCGGTTGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGCTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCggggaagttggggct              >  2:2690/1‑239 (MQ=255)
                                                                                                                                                        tCTGGCACGCCTGGCTGGCAGGGATTGTCGGATGTCTGCTGGGCGACTGGATCTCCTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGCCATTTCTGAAGAAAAACAAAGCATTACTCGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTCGTCGGTCGCTTTGTTGGCCCGACGCGTCCGCTGGTGACAATGGAGGCGGGAATGCTGGATCTGCCGGTGGcaaa  >  1:158514/1‑247 (MQ=255)
                                                                                                                                                                                          ||                                                                                                                                                                                                                       
GGAACACTTTATTACCCAATCCACCGTGTATTCATTGATGGCGGTGGTGTTGGTGGCCTTTCTGGAGTCGCTGGCGCTGGTCGGTTTGATTCTACCCGGTACGGTGCTGATGGCGGGGCTGGGAGCGCTGATTGGCAGCGGCGAGTTAAGTTTCTGGCACGCCTGGCTGGCAGGGATTATTGGCTGCTTGATGGGCGACTGGATTTCTTTCTGGCTGGGTTGGCGTTTTAAAAAGCCGTTGCATCGCTGGTCATTTCTGAAGAAAAACAAAGCACTACTTGATAAAACTGAACATGCGTTGCATCAACACAGCATGTTCACCATTCTGGTCGGTCGTTTTGTTGGCCCGACGCGTCCGCTGGTGCCAATGGTGGCGGGAATGCTGGATCTGCCGGTGGCTAA  >  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: