Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,351	A→G	16.4%	I11T (ATC→ACC)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,351	0	A	G	16.4%	58.2 / 4.3	28	I11T (ATC→ACC)	trmE	GTPase
Reads supporting (aligned to +/- strand):  ref base A (14/9);  new base G (2/2);  total (16/12)
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.29e-01

TTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTC  >  NC_002947/6215‑6500                                                                                                                                         |                                                                                                                                                      ttCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGt                                                                                                                                          <  2:247655/150‑1 (MQ=255)      cgtcgCGAAATGGACAGTGATGGGCACGCGTCGGGTGCAGGGTGCGGCTGCTGACCAACAGGCCGGCCTTGGCAGTCAGAGGGCCGGATCCCCTGACAATTCCCTCTCCTCCTCGGCCCCGGCCCGTGGCGTTTGCTGCGATGGCTTctc                                                                                                                                     <  2:228/150‑3 (MQ=255)                  gACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCa                                                                                                                         >  2:189821/1‑150 (MQ=255)                     cGTAAGGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGGTGAACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGc                                                                                                                      <  1:613317/150‑1 (MQ=255)                           gggCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                                                <  2:699259/150‑1 (MQ=255)                                  ccGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTa                                                                                                         <  2:452713/150‑1 (MQ=255)                                            aGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAAc                                                                                               <  1:826515/150‑1 (MQ=255)                                                                             gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCTTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGAt                                                                                        <  1:83427/124‑1 (MQ=255)                                                                             gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCTTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGAt                                                                                        >  2:83427/1‑124 (MQ=255)                                                                                   agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTa                                                        <  1:189821/150‑1 (MQ=255)                                                                                       ggCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTca                                                    <  1:166540/150‑1 (MQ=255)                                                                                                 ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCg                                          >  1:183204/1‑150 (MQ=255)                                                                                                   tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTa                                        >  2:693240/1‑150 (MQ=255)                                                                                                    gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGGTGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAACCGGCCCACTAGGGGGCGCTTTTTATTCACAGGCTcgcggct                                       >  1:684912/1‑147 (MQ=255)                                                                                                       aaTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGACGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAg                                    >  1:419265/1‑150 (MQ=255)                                                                                                          tCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                                 <  2:183204/150‑1 (MQ=255)                                                                                                            ccACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGt                               >  1:205856/1‑150 (MQ=255)                                                                                                            ccACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGt                               >  2:748097/1‑150 (MQ=255)                                                                                                             cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc                              >  2:261766/1‑150 (MQ=255)                                                                                                              ccatcGCCACGGCCCAGGGCCGGGGCGGTGGTGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaaaaa                                                                                <  2:683952/95‑1 (MQ=37)                                                                                                              ccatcGCCACGGCCCAGGGCCGGGGCGGTGGTGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCaaaaa                                                                                >  1:683952/5‑99 (MQ=37)                                                                                                              acacCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCa                             >  1:803994/1‑150 (MQ=255)                                                                                                               catcGCCACGGCCCAGGGCCGGGGCGGTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCaa                            >  2:27198/4‑150 (MQ=255)                                                                                                                atcGCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                           >  1:19240/3‑150 (MQ=255)                                                                                                                acagCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                           >  2:594534/4‑150 (MQ=255)                                                                                                                  cGCCACGGCCCAGGGCCGGGGCGGTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagca                         <  1:27198/150‑1 (MQ=255)                                                                                                                      ccGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcag                     >  2:786364/1‑150 (MQ=255)                                                                                                                       cGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcagc                    >  1:191338/1‑150 (MQ=255)                                                                                                                                         aTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTc  <  1:786364/150‑1 (MQ=255)                                                                                                                                         |                                                                                                                                                      TTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTC  >  NC_002947/6215‑6500

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: ‑