Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,351,095	C→G	11.4%	P11A (CCA→GCA)	arfB →	peptidyl‑tRNA hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,351,095	0	C	G	11.4%	101.7 / 5.2	35	P11A (CCA→GCA)	arfB	peptidyl‑tRNA hydrolase
Reads supporting (aligned to +/- strand):  ref base C (16/15);  new base G (2/2);  total (18/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCATCGTTGCCCGAGTTTTACAAGGCGCGGCTGCTGGCG  >  NC_002947/3350949‑3351238                                                                                                                                                   |                                                                                                                                                cccTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGc                                                                                                                                                 >  8:128704/1‑147 (MQ=255)        gAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCg                                                                                                                                        <  7:128704/149‑1 (MQ=255)             gAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGGCAGATGcacgtta                                                                                                                                    <  8:189471/148‑7 (MQ=255)              aGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCg                                                                                                                                  <  6:285759/149‑1 (MQ=255)                 gCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAAc                                                                                                                               <  4:308495/149‑1 (MQ=255)                           gTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGaa                                                                                                                                <  4:144813/138‑1 (MQ=255)                               gCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGGCAGATGc                                                                                                                                          <  8:199080/124‑1 (MQ=255)                               gCTTGCGCACTGATTCAAAGAAAACTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGGCAGATGc                                                                                                                                          >  7:199080/1‑124 (MQ=255)                                cTTGCGCCCTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCcgcg                                                                                                                >  5:301976/1‑149 (MQ=255)                                     gcACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACACCcgcgccttg                                                                                                           >  8:235456/1‑145 (MQ=255)                                            ttCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCgg                                                                                                    <  3:138466/149‑1 (MQ=255)                                              cAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGt                                                                                                  >  1:190371/1‑149 (MQ=255)                                                cgatctccTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTgg                                                                                                <  7:281739/143‑1 (MQ=255)                                                 aGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTggg                                                                                               <  8:19805/149‑1 (MQ=255)                                                       cTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATGTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTgg                                                                                                >  8:281739/1‑142 (MQ=255)                                                           aTTGTTCTATGGGGGAGCGGCACATTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTc                                                                                      <  5:22632/148‑1 (MQ=255)                                                                               cacaGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTa                                                                 >  5:126294/1‑149 (MQ=255)                                                                               cacaGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTa                                                                 <  6:301976/149‑1 (MQ=255)                                                                                      tgtgATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTAc                                                          >  6:38224/1‑149 (MQ=255)                                                                                      tgtgATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGCCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTAc                                                          >  6:102883/1‑149 (MQ=255)                                                                                        tgATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCgg                                                                   >  1:151680/1‑138 (MQ=255)                                                                                        tgATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCgg                                                                   <  2:151680/138‑1 (MQ=255)                                                                                        tgATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGc                                                        <  2:190371/149‑1 (MQ=255)                                                                                          aTAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCtt                                                      >  6:42937/1‑149 (MQ=255)                                                                                              tGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGa                                                   >  8:174027/1‑148 (MQ=255)                                                                                              tGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGa                                                   >  3:21081/1‑148 (MQ=255)                                                                                                       aCAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGcc                                         <  5:42937/149‑1 (MQ=255)                                                                                                         aaGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGAGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCt                                        >  7:263442/1‑148 (MQ=255)                                                                                                                   tCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCGTCGTTGccc                             <  3:48691/149‑1 (MQ=255)                                                                                                                             tATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCATCGTTGCCCGAGTTTTAc                    >  4:269715/1‑148 (MQ=255)                                                                                                                                  ccAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCATCGTTGCCCGAGTTTTACAAGgcg              >  3:123630/1‑149 (MQ=255)                                                                                                                                         tcggCATCTGGCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCATCGTTGCCCGAGTTTTACAAGGCGCGgctgc       >  8:138753/4‑149 (MQ=255)                                                                                                                                          gTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCATCGTTGCCCGAGTTTTACAAGGCGCGgctgct      <  6:126294/149‑1 (MQ=255)                                                                                                                                          gTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCAGCGTTGCCCGAGTTTTACAAGGCGCGgctgct      >  7:212496/1‑149 (MQ=255)                                                                                                                                             cATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCATCGTTGCCCGAGTTTTACAAGGCGCGGCTGCTGgc   <  4:123630/149‑1 (MQ=255)                                                                                                                                              aTCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCATCGTTGCCCGAGTTTTACAAGGCGCGGCTGCTGgcg  >  2:182676/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                CCCTGACGAGTGGAGGGCTCGCTTGGGTCGGCTTGCGCACTGATTCAAAGAACCCTGCATTGTGCTATGGGGGAGCGGCACAGTCTGTGATAATGGGCAATGACAAGGGCCCTATCATGCTGACTATCTCCAATAACGTGCATCTGCCAGATGCCGAAATCGAACTCACATACATCCGCGCGCAAGGCGCGGGTGGGCAGAATGTCAACAAAGTGTCCAGCGCGGTACACCTACGCTTCGACATCCCGGCCTCATCGTTGCCCGAGTTTTACAAGGCGCGGCTGCTGGCG  >  NC_002947/3350949‑3351238

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