Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	547,739	G→A	10.8%	D196N (GAC→AAC)	fusA →	elongation factor G 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	547,739	0	G	A	10.8%	96.1 / 3.8	37	D196N (GAC→AAC)	fusA	elongation factor G 1
Reads supporting (aligned to +/- strand):  ref base G (14/19);  new base A (2/2);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

CAAACTTCCTGCGCGTCGTAGGTCAGATCAAGAACCGTCTGGGTCACACTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACAAGTACCTTGAAGAA  >  NC_002947/547599‑547873                                                                                                                                             |                                                                                                                                       cAAACTTACTGCGCGGCGTAGGTCAGATCAAGAACCGTCTGGGTCACAACCCGGAACCTTTGCTGCGCCGCTACGGTTCGGAAGATAACCGCCAGGGTCAGGTTGACCTGAGCAAGATGAAGGCTATCTACTGGAACGACGACGACAAgg                                                                                                                               <  1:551320/150‑1 (MQ=255)           gcgcGTCGTAGGTCAGATCAAGAACCGTCTGGGTCACACTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTActaccta                                                                                                                     <  1:508058/150‑1 (MQ=255)                  gtAGGTCAGATCAAGAACCGTCTGGGTCACACTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAg                                                                                                              <  2:462335/150‑1 (MQ=255)                     ggTCAGATCAAGAACCGTCTGGGTCACACTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGgaa                                                                                                           <  1:399061/150‑1 (MQ=255)                         aGATCAAGAACCGTCTGGGTCACACTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGaaa                                                                                                       >  2:472414/1‑150 (MQ=255)                               aGAACCGTCTGGGTCACACTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTg                                                                                                 >  2:359379/1‑150 (MQ=255)                                gAACCGTCTGGGTCACACTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGc                                                                                                <  2:563116/150‑1 (MQ=255)                                                acTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAg                                                                                                              <  2:484102/120‑1 (MQ=255)                                                acTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAg                                                                                                              >  1:484102/1‑120 (MQ=255)                                                      gTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAAc                                                                          <  2:376380/150‑1 (MQ=255)                                                         ccTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGaa                                                                       <  2:293260/150‑1 (MQ=255)                                                           tGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATg                                                                     <  1:508128/150‑1 (MQ=255)                                                           tGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATg                                                                     <  1:555173/150‑1 (MQ=255)                                                           tGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACAACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATg                                                                     >  2:452976/1‑150 (MQ=255)                                                           tGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACAACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATg                                                                     >  2:452896/1‑150 (MQ=255)                                                             ttCAGCTGGCCCTCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGgc                                                                   >  2:375881/1‑150 (MQ=255)                                                             ttCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGgc                                                                   >  1:67689/1‑150 (MQ=255)                                                             ttCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGgc                                                                   >  2:244484/1‑150 (MQ=255)                                                                   tGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGcaac                                                              >  2:304837/1‑149 (MQ=255)                                                                   tGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGcaac                                                              <  1:304837/149‑1 (MQ=255)                                                                      ccATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGAGGACAAAGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGACATGGTCGAGCTGGCTACCAAATCGCGCAacaaca                                                          >  2:188962/1‑150 (MQ=255)                                                                       cATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACAt                                                         >  2:349693/1‑150 (MQ=255)                                                                            gTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTg                                                    >  2:528123/1‑150 (MQ=255)                                                                              tCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGaa                                                  <  1:359379/150‑1 (MQ=255)                                                                              tCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGaa                                                  <  1:528123/150‑1 (MQ=255)                                                                               cGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAGGATGAAGGCTATCTACTGGAACGACAACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTCGAAg                                                 <  1:452896/150‑1 (MQ=255)                                                                               cGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACAACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTCGAAg                                                 <  1:452976/150‑1 (MQ=255)                                                                                  aaGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAgccg                                              >  1:470722/1‑150 (MQ=255)                                                                                       aaCTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAgccgccgcc                                         >  1:193827/1‑150 (MQ=255)                                                                                                          aCCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCtga                      <  2:67689/150‑1 (MQ=255)                                                                                                                 cAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACaa               <  1:244484/150‑1 (MQ=255)                                                                                                                   aGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACAAGt             >  2:77289/1‑150 (MQ=255)                                                                                                                           gCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACAAGTACCTTgaa     <  1:349693/150‑1 (MQ=255)                                                                                                                             tATCTACTGGAACGACGACGACAAGGGTACTACCTACCGGGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACAAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACAAGTACCTTgaaga   <  2:181685/150‑1 (MQ=255)                                                                                                                              aTCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGAATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGCAGCCACCGAAGAGCTGATGAACAAGTACCTTgaagaa  <  1:375881/150‑1 (MQ=255)                                                                                                                                       aaCGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACAAGTACCTTgaa     >  1:59382/1‑138 (MQ=255)                                                                                                                                       aaCGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACAAGTACCTTgaa     <  2:59382/138‑1 (MQ=255)                                                                                                                                             gacgacAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACAAGTACCTTg       <  2:19761/130‑1 (MQ=255)                                                                                                                                             gacgacAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACAAGTACCTTg       >  1:19761/1‑130 (MQ=255)                                                                                                                                             |                                                                                                                                       CAAACTTCCTGCGCGTCGTAGGTCAGATCAAGAACCGTCTGGGTCACACTCCGGTTCCTGTTCAGCTGGCCATCGGTTCGGAAGATAACTTCCAGGGTCAGGTTGACCTGATCAAGATGAAGGCTATCTACTGGAACGACGACGACAAGGGTACTACCTACCGCGAGGAAGAAATTCCTGCGGATATGGTCGAGCTGGCTAACGAATGGCGCAACAACATGGTTGAAGCCGCCGCCGAAGCCACCGAAGAGCTGATGAACAAGTACCTTGAAGAA  >  NC_002947/547599‑547873

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