Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,707,816	A→G	11.1%	V353A (GTG→GCG)	actP‑III ←	acetate permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,707,816	0	A	G	11.1%	89.3 / 3.0	36	V353A (GTG→GCG)	actP‑III	acetate permease
Reads supporting (aligned to +/- strand):  ref base A (19/13);  new base G (2/2);  total (21/15)
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.41e-01

AAGGCAATGTTCTGCGACTCGAACATCAGGCCCAGGATCACCGCCAGCAGGCCGATGAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGCCAGGTGCACGGCGACCATGTTGCCGCCACC  >  NC_002947/3707672‑3707952                                                                                                                                                 |                                                                                                                                          aaggCAATGTTCTGCGACTCGAACATCAGGCCCAGGATCACCGCCAGCAGGCCGATGAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACAccgc                                                                                                                                       <  3:116792/149‑1 (MQ=255)                  cTCGAACATCAGGCCCAGGATCACCGCCAGCAGGCCGATGAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAgcg                                                                                                                      >  7:27118/1‑149 (MQ=255)                                ccAGGATCACCGCCAGCAGGCCGATGAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGAc                                                                                                         >  4:96671/1‑148 (MQ=255)                                ccAGGATCACCGCCAGCAGGCCGATGAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGAc                                                                                                         >  2:29880/1‑148 (MQ=255)                                ccAGGATCACCGCCAGCAGGCCCATGAGCAACGTCGCAATGCGAGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCAAGGCCAGCGAc                                                                                                         >  2:161483/1‑148 (MQ=255)                                 cAGGATCACCGCCAGCAGGCCGATGAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACca                                                                                                       >  1:84530/1‑149 (MQ=255)                                         ccGCCAGCAGGCCGATGAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAg                                                                                                >  6:192694/1‑148 (MQ=255)                                          cGCCAGCAGGCCGATGAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACa                                                                                                                         <  4:49334/122‑1 (MQ=255)                                                         gAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACGCCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCa                                                                              >  2:51837/1‑149 (MQ=255)                                                         gAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACGCCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCa                                                                              >  2:55602/1‑149 (MQ=255)                                                               cGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCt                                                                         >  1:207183/1‑148 (MQ=255)                                                                 tGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCt                                                                         >  5:167932/1‑146 (MQ=255)                                                                    cAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGAt                                                                    >  1:28477/1‑148 (MQ=255)                                                                     aaTGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAAtggtg                                                                                         <  1:62709/127‑1 (MQ=255)                                                                     aaTGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAAtggtg                                                                                         >  2:62709/1‑127 (MQ=255)                                                                             acacGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAg                                                           <  2:84530/148‑1 (MQ=255)                                                                                 gcgcATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGaa                                                        >  1:152368/1‑147 (MQ=255)                                                                                      tttCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTg                                                 <  2:207183/149‑1 (MQ=255)                                                                                         ccTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACGCCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGcc                                               <  1:55602/148‑1 (MQ=255)                                                                                         ccTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACGCCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGcc                                               <  1:51837/148‑1 (MQ=255)                                                                                           tctTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTgaga                                                                     >  6:91901/1‑124 (MQ=255)                                                                                           tctTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTgaga                                                                     <  5:91901/124‑1 (MQ=255)                                                                                                cTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGGCCCACGGCTGAGATGAAGCCGAGGAACAGGTTGcc                                               <  1:172511/142‑1 (MQ=255)                                                                                                cTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGcc                                               >  2:172511/1‑141 (MQ=255)                                                                                                 tCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGtt                                                  <  2:213247/137‑1 (MQ=255)                                                                                                 tCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGtt                                                  >  1:213247/1‑137 (MQ=255)                                                                                                     tCGCCTTGCCTTTGCGCAGCACGCAGGCATACTGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCttg                                  <  2:218719/149‑3 (MQ=255)                                                                                                       gccttgccttTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTggg                                 <  8:27118/148‑1 (MQ=255)                                                                                                        ccttgccttTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGc                                >  2:77870/1‑148 (MQ=255)                                                                                                            gccttTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGCCAg                             <  5:77816/147‑1 (MQ=255)                                                                                                                   cgcATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGCCAGGTGCACGGc                    >  5:145960/1‑149 (MQ=255)                                                                                                                    gcATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGCCAGGTGCACGGc                    >  8:240151/1‑148 (MQ=255)                                                                                                                          cGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGCCAGGTGCACGGCGACCATg             <  2:152368/149‑1 (MQ=255)                                                                                                                               gCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGCCAGGTGCACGGCGACCATGTTgcc        <  3:96671/149‑1 (MQ=255)                                                                                                                                   aCAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGCCAGGTGCACGGCGACCATGTTGCCGCCa    >  5:130459/1‑149 (MQ=255)                                                                                                                                      ggTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGCCAGGTGCACGGCGACCATGTTGCCGCCAcc  >  7:226623/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                          AAGGCAATGTTCTGCGACTCGAACATCAGGCCCAGGATCACCGCCAGCAGGCCGATGAGCAACGTGGCAATGCGTGACACGCGCATTTCCTCTTGCTCGGTCGCCTTGCCTTTGCGCATCACGCAGGCATACAGGTCGTGGGACACCGCCGATGCGCCGGACAGCGCCAGGCCAGCGACCACAGCCAGAATGGTGGCGAAGG‑CCACGGCTGAGATGAAGCCGAGGAACAGGTTGCCACCCACGGCCTGGGCCAGGTGCACGGCGACCATGTTGCCGCCACC  >  NC_002947/3707672‑3707952

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