Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,286,496 C→A 12.1% A155A (GCC→GCA PP_2888 → ECF family RNA polymerase sigma factor PrtI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,286,4960CA12.1% 86.1 / 3.0 33A155A (GCC→GCAPP_2888ECF family RNA polymerase sigma factor PrtI
Reads supporting (aligned to +/- strand):  ref base C (12/17);  new base A (2/2);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.11e-01

CCTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTGCAG  >  NC_002947/3286349‑3286633
                                                                                                                                                   |                                                                                                                                         
ccTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCtt                                                                                                                                         <  2:443149/150‑1 (MQ=255)
           ttGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGa                                                                                                                              <  2:589703/150‑1 (MQ=255)
               ccGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGa                                                                                                                          >  2:193266/1‑150 (MQ=255)
                      aCGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGcaac                                                                                                                   <  1:536933/150‑1 (MQ=255)
                         gcggAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGcaacccc                                                                                                                <  2:539242/150‑1 (MQ=255)
                                agcGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGcaacccccaacccc                                                                                                         <  1:193266/150‑1 (MQ=255)
                                    tGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGc                                                                                                     <  1:267797/150‑1 (MQ=255)
                                                         cATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGAcg                                                                                >  1:50847/1‑150 (MQ=255)
                                                             gAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGc                                                                            >  2:347631/1‑150 (MQ=255)
                                                              aaGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCt                                                                           >  1:447649/1‑150 (MQ=255)
                                                               aGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGACTg                                                                          >  1:606113/1‑150 (MQ=255)
                                                                           cAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCgat                                                              <  1:347631/150‑1 (MQ=255)
                                                                             aggaggCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCgatga                                                            >  1:392494/1‑150 (MQ=255)
                                                                                 ggCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGACACTGaaa                                                                                     <  1:426606/121‑1 (MQ=255)
                                                                                 ggCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGACACTGaaa                                                                                     >  2:426606/1‑121 (MQ=255)
                                                                                  gccgccCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAAt                                                       <  2:261665/150‑1 (MQ=255)
                                                                                              ttGGGCATCCCCATCGGCACTGTGACGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCgatga                                                            >  2:153714/1‑133 (MQ=255)
                                                                                              ttGGGCATCCCCATCGGCACTGTGACGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCgatga                                                            <  1:153714/133‑1 (MQ=255)
                                                                                               tGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTtcgt                                          <  1:456661/150‑1 (MQ=255)
                                                                                                                aCTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAAc                         <  2:50847/150‑1 (MQ=255)
                                                                                                                  tgtgATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGgaga                                                                                           >  2:189833/1‑82 (MQ=255)
                                                                                                                  tgtgATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGgaga                                                                                           <  1:189833/82‑1 (MQ=255)
                                                                                                                  tgtgATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCt                       <  2:447649/150‑1 (MQ=255)
                                                                                                                    tgATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCtgt                     >  2:250008/1‑150 (MQ=255)
                                                                                                                             gcCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGc            >  1:67983/1‑150 (MQ=255)
                                                                                                                                  tCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGTGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAgg       <  1:385666/150‑1 (MQ=255)
                                                                                                                                  tCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAgg       >  2:277182/1‑150 (MQ=255)
                                                                                                                                   cgcgcgcgcgCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGt      <  2:36327/150‑1 (MQ=255)
                                                                                                                                   cgcgcgccccCAACGCCTTGCGCGCACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGt      <  2:36298/150‑1 (MQ=255)
                                                                                                                                    ccgatctcgcAAGGCATTGCGCGAACTGACCGAGGGCAACCCCCAAAACCCGTCCGTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTg     <  2:418008/143‑1 (MQ=255)
                                                                                                                                      gttcgcgcAAGGCATTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTg                             <  2:10624/121‑1 (MQ=255)
                                                                                                                                      gttcgcgcAAGGCATTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGAACCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTg                             >  1:10624/4‑124 (MQ=255)
                                                                                                                                       ttcgcgcAAGGCATTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCAAGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTGCAg  >  2:222973/3‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
CCTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTGCAG  >  NC_002947/3286349‑3286633

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