Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	2,817,046	Δ1 bp	noncoding (43/77 nt)	PP_t44 →	tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,817,045	0	A	.	95.5%	92.3 / ‑3.2	22	noncoding (42/77 nt)	PP_t44	tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (8/13);  total (8/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.58e-01

ATGAGTGGAAAATTTTTTGAGAAAAAGCTTTCATTTATCAGAGGGTTAGGGTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTT  >  NC_002947/2816902‑2817186                                                                                                                                                |                                                                                                                                              atgaGTGGAAAATTTTTTGAGAAAAAGCTTTCATTTATCAGAGCGGTAGGGTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAggg                                                                                                                                            <  1:291821/147‑1 (MQ=255)    aGTGGAAAATTTTTTGAGAAAAAGCTTTCATTTATCAGACGGTTAGGGTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCg                                                                                                                                       <  4:293452/148‑1 (MQ=255)      tGGAAAATTTTTTGAGAAAAAGCTTTCATTTATCAGAGGGTTAGGGTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAg                                                                                                                                     <  8:57329/148‑1 (MQ=255)                         aaGCTTTCATTTATCAGAGGGTTAGGGTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTc                                                                                                                   >  5:245919/1‑147 (MQ=255)                                  tttATCAGAGGGTTAGGGTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCaaa                                                                                                          >  3:241342/1‑147 (MQ=255)                                      tCAGAGGGTTAGGGTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAAcc                                                                                                     <  8:210044/148‑1 (MQ=255)                                         gagGGTTAGGGTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCt                                                                                                   <  6:245919/147‑1 (MQ=255)                                                 gggTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATc                                                                                           <  2:32886/147‑1 (MQ=255)                                                  ggTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGGCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCa                                                                                         <  7:129128/148‑1 (MQ=255)                                                      cATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCa                                                                                     <  2:179636/148‑1 (MQ=255)                                                               aCGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGt                                                                            <  4:205094/148‑1 (MQ=255)                                                                             aTTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGAtttttt                                                               >  1:13899/1‑147 (MQ=255)                                                                                                    aGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTctcg                                       >  3:15867/1‑148 (MQ=255)                                                                                                          ggCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGaa                                 >  1:49389/1‑148 (MQ=255)                                                                                                                gcgcAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAttttg                            <  8:296179/147‑1 (MQ=255)                                                                                                                   cAGCCCGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAttttgtgcc                        >  4:164735/1‑148 (MQ=255)                                                                                                                       ccGGTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAttttgtgccactt                    >  3:151317/1‑148 (MQ=255)                                                                                                                         ggTAGCGCACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAttttgtgccactttt                  <  3:164735/148‑1 (MQ=255)                                                                                                                             gcgcACTTGCATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAAttttgtgccactttttgtg              <  2:49389/148‑1 (MQ=255)                                                                                                                                     gCATGGGGTTC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTGTTCTGTCTCGCTCGAATTTTGTGCCACCTTTTGTGCCACTg        <  7:274377/146‑1 (MQ=255)                                                                                                                                      cATGGGGTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTTTGCCACTTTTTGTGCCACTGGGc     >  8:202869/1‑148 (MQ=255)                                                                                                                                          gggTGC‑AGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCtat  <  4:151317/147‑3 (MQ=255)                                                                                                                                                |                                                                                                                                              ATGAGTGGAAAATTTTTTGAGAAAAAGCTTTCATTTATCAGAGGGTTAGGGTACATTCTTCAACGTTCTCGGCAAGATTGAGAACATGCTTCAAAGCCCAGTCGGGGCGTAGCGCAGCCCGGTAGCGCACTTGCATGGGGTGCAAGGGGTCGAGTGTTCGAATCACTCCGTCCCGACCAAAAAACCCTAGAAAATCCAGTCACTTAACGGTGACTGGATTTTTTTATGCCTGATGGTTTTCTGTCTCGCTCGAATTTTGTGCCACTTTTTGTGCCACTGGGCTTT  >  NC_002947/2816902‑2817186

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