Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,286,493 T→G 11.8% N154K (AAT→AAG 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,4930TG11.8% 78.2 / 2.5 33N154K (AAT→AAGPP_2888ECF family RNA polymerase sigma factor PrtI
Reads supporting (aligned to +/- strand):  ref base T (12/16);  new base G (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 = 1.00e+00

GCCGACCTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTGCAG  >  NC_002947/3286344‑3286633
                                                                                                                                                     |                                                                                                                                            
gCCGACCTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAAt                                                                                                                                              >  2:385666/1‑150 (MQ=255)
     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                                                            <  1:153714/133‑1 (MQ=255)
                                                                                                   ttGGGCATCCCCATCGGCACTGTGACGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCgatga                                                            >  2:153714/1‑133 (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)
                                                                                                                                                     |                                                                                                                                            
GCCGACCTGCAGGCCTTTGCCCGGCTGACGGCGGAGCAGCGTGCCTTGTTGTTGCTGGTCAGCATCGAAGGCTTGAGTTACAAGGAGGCCGCCCAGGCCTTGGGCATCCCCATCGGCACTGTGATGTCGCGCCTGTCGCGCGCGCGCAATGCCTTGCGCGAACTGACCGAGGGCAACCCCCAACCCCCGGCCTTGCGGAGACTGAAATGACGCGGCTGATCCCCAGCGATGATGAATTGCACGCTTTCGTCGATGAACGCCTGGAACCTGTGCGCCGGGCCGAGGTGCAG  >  NC_002947/3286344‑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: