Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,346	T→G	10.2%	T13P (ACG→CCG)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,346	0	T	G	10.2%	102.1 / 2.7	41	T13P (ACG→CCG)	trmE	GTPase
Reads supporting (aligned to +/- strand):  ref base T (19/17);  new base G (2/3);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAG  >  NC_002947/6202‑6481                                                                                                                                                 |                                                                                                                                         ccAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGa                                                                                                                                     >  2:1039721/1‑150 (MQ=255)     cACCATCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCggtgtg                                                                                                                                 <  2:134683/150‑6 (MQ=255)      aCCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATggcg                                                                                                                                >  1:796292/1‑150 (MQ=255)      aCCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCggtgtgg                                                                                                                                <  2:529411/150‑7 (MQ=255)          gCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGa                                                                                                                            <  2:926070/150‑1 (MQ=255)               tCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGtt                                                                                                                       <  2:231887/150‑1 (MQ=255)                  tcgtcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTTAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCa                                                                                                                    >  2:729227/1‑150 (MQ=255)                  atctcgCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACCCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTCCa                                                                                                                    <  1:985723/147‑1 (MQ=255)                       gcgGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGGTGCTATGGCGGCGATGGTTCCAcagat                                                                                                               >  2:985723/1‑146 (MQ=255)                           aaTGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAgtg                                                                                                           >  2:401929/1‑150 (MQ=255)                                 ccGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATg                                                                                                     >  1:287148/1‑150 (MQ=255)                                      aTGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGa                                                                                                <  1:930899/150‑1 (MQ=255)                                       tGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCggtgtg                                                                                                                                 <  1:88384/116‑6 (MQ=255)                                       tGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCCAGGGCCGGGGCGGTgtg                                                                                                                                 >  2:88384/1‑114 (MQ=255)                                        gggCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGaaa                                                                                              <  1:435654/150‑1 (MQ=255)                                          gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                            >  1:433019/1‑150 (MQ=255)                                              gCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCctct                                                                                        >  1:90108/1‑150 (MQ=255)                                                       tCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACaaaa                                                                               >  1:580762/1‑150 (MQ=255)                                                               cGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGa                                                                       >  2:484132/1‑150 (MQ=255)                                                                         tCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAAc                                                             >  2:803613/1‑150 (MQ=255)                                                                            aCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                          <  1:1039721/150‑1 (MQ=255)                                                                            aCAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACAACAAAAACGACAGATAGCAAAAACGcc                                                          >  1:241358/1‑150 (MQ=255)                                                                             cAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGccc                                                         <  1:803613/150‑1 (MQ=255)                                                                                  cggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGGAATCCTCGACGACAAAAACAACAGGTAGCAAAAACGCCCCACTa                                                    >  2:124874/1‑150 (MQ=255)                                                                                   ggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                              <  2:226922/91‑1 (MQ=255)                                                                                   ggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                              >  1:226922/1‑91 (MQ=255)                                                                                   ggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg                                                   <  1:881100/150‑1 (MQ=255)                                                                                    gccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACAACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑g                                                 <  2:241358/150‑1 (MQ=255)                                                                                     ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑gg                                                <  1:729227/150‑1 (MQ=255)                                                                                      cTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑ggg                                               <  2:433019/150‑1 (MQ=255)                                                                                            aGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGttt                                         <  2:580762/150‑1 (MQ=255)                                                                                              ccAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCC‑CTAGTGGGGCGttttt                                       <  1:204918/150‑1 (MQ=255)                                                                                                agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTa                                     <  1:401929/150‑1 (MQ=255)                                                                                                       cGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAg                              <  1:541379/150‑1 (MQ=255)                                                                                                            aaCCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGGTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGCAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCCCGAT‑GGGGCGTGTTTTATTCACAGGCTAg                         >  1:726639/1‑150 (MQ=255)                                                                                                                 gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTAt                    >  1:927649/1‑150 (MQ=255)                                                                                                                          cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc           >  2:349293/1‑150 (MQ=255)                                                                                                                              tcGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc       >  1:842722/2‑150 (MQ=255)                                                                                                                              ccGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc       <  2:90108/150‑1 (MQ=255)                                                                                                                                   ccGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCCAAAACGCACCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcag  >  2:151263/1‑150 (MQ=255)                                                                                                                                   ccGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcag  >  1:871256/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                         CCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAG  >  NC_002947/6202‑6481

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