Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,390,853	C→A	23.5%	C218* (TGC→TGA)	aspS →	aspartate‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,390,853	0	C	A	23.5%	29.8 / 3.7	17	C218* (TGC→TGA)	aspS	aspartate‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base A (2/2);  total (9/8)
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.64e-01

CCGTGCCACCCCGGAAGGCGCCCGCGACTACCTGGTGCCAAGCCGTACCCACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTG  >  NC_002947/1390706‑1390949                                                                                                                                                    |                                                                                                 ccGTGCCACCCCGGAAGGCGCCCGCGACTACCTGGTGCCAAGCCGTACCCACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGAtt                                                                                                >  2:491400/1‑150 (MQ=255)       cACCCCGGAAGGCGCCCGCGACTACTTGGTGCCAAGCCGTACCCACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGa                                                                                          <  1:271524/150‑1 (MQ=255)                            cTACCTGGTGCCAAGCCGTACCCACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTGCGTGCCGACCg                                                                     >  1:237449/1‑150 (MQ=255)                               ccTGGTGCCAAGCCGTACCCACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGATTCCGTGACGAAGACCTGCGTGCCGACCGTCa                                                                  <  1:491400/150‑1 (MQ=255)                                       cAAGCCGTACCCACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGATTCCGTGACGAAGACCTGCGTGCCGACCGTCAGCCGGAAt                                                          >  1:52544/1‑150 (MQ=255)                                           ccGTACCCACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGATTCCGTGACGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCAc                                                      <  2:52544/150‑1 (MQ=255)                                                 ccACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTGCGTGCCGACCGt                                                                    >  1:42598/1‑130 (MQ=255)                                                 ccACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTGCGTGCCGACCGt                                                                    <  2:42598/130‑1 (MQ=255)                                                                                    cacaGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGAc                                                                                         <  2:55467/74‑1 (MQ=255)                                                                                    cacaGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGAc                                                                                         >  1:55467/1‑74 (MQ=255)                                                                                        gCTGTTCAAGCAGCTTCTGATGGTCGCTGGCTTCCATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTTCTTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCAt         <  2:241068/150‑1 (MQ=255)                                                                                            ttCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGc     >  1:108823/1‑150 (MQ=255)                                                                                               aaGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTg  >  1:325396/1‑150 (MQ=255)                                                                                                                       ttCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTGCGTGCCGAcc                                                                      <  2:164540/58‑1 (MQ=255)                                                                                                                       ttCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTGCGTGCCGAcc                                                                      >  1:164540/1‑58 (MQ=255)                                                                                                                                                  tGCTTCCGTGACGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATg        >  1:442591/1‑93 (MQ=255)                                                                                                                                                  tGCTTCCGTGACGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATg        <  2:442591/93‑1 (MQ=255)                                                                                                                                                    |                                                                                                 CCGTGCCACCCCGGAAGGCGCCCGCGACTACCTGGTGCCAAGCCGTACCCACGCCGGTAGCTTCTTCGCCCTGCCGCAATCGCCACAGCTGTTCAAGCAGCTGCTGATGGTCGCTGGCTTCGATCGTTACTACCAGATCGCCAAGTGCTTCCGTGACGAAGACCTGCGTGCCGACCGTCAGCCGGAATTCACCCAGATCGACATCGAGACCAGTTTCCTCGATGAAAGCGAGATCATGGGCCTG  >  NC_002947/1390706‑1390949

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