Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,287,294	C→G	noncoding (35/85 nt)	tyrV ←	tRNA‑Tyr

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,294	0	C	G	95.2%	58.9 / NA	22	noncoding (35/85 nt)	tyrV	tRNA‑Tyr
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (10/10);  minor base T (0/1);  total (11/11)
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.58e-01

GATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCG‑TTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATT  >  NC_000913/1287064‑1287536                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    gATGGTGGTGGGGGAAGGATTACTCAGCGCTCCGCGATTCGCCCTTCGGGTCGTTGCCTGTGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGATGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGCTTCGTTGGGACGTTCAGGGACTTTTGAAAGTGATGGTGGTTGGGGAAGGTTTCGAACCTTCGAAGTCGATGACGGCAGCTTTAGGGTCTGCTCCCTTTGGCCgt                                                                                                                                                                                                                                 >  2:90477/1‑250 (MQ=18)           ggggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCTGCTCGGGAAccc                                                                                                                                                                                                                       >  2:74841/1‑251 (MQ=255)              gAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGCTTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCccac                                                                                                                                                                                                                    >  2:71944/1‑251 (MQ=255)                 ggagtACTAAGCGCTGCCCCCGTCGCCCTTCGGCTCGTTGCCTGCGGCACTCCTCTCTCGCTGGCGCTCCAGTCGAACCTTGGTCTAAGCTTCTCATCCTTCCCCGCATGGTCAGAATATTTAATTGCGGATTCGTAGGAAATTTCAAGTACTTGTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTGCGAAGTCGATGACGGCAGATTTATAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAg                                                                                                                                                                                                                 <  2:113447/247‑1 (MQ=21)                             gcTGCGCGCTTCTCCCGTCGGGGCGTTGCCTGCGGCAACGGTCGCGCGCTGGCGATCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTGAATTGCGGAGTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCaaa                                                                                                                                                                                                     <  2:222293/251‑1 (MQ=25)                                       tCGCCCTTCGGGTCGTTGCCTGCGTGAACGCTCTCTCGCTGGCTCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAATTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTGTGGCCGCTCGGGATCCCCACCAGGGGTAATTCAAATTTTGAg                                                                                                                                                                                              <  1:86307/248‑1 (MQ=25)                                        cGCCCGTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTaa                                                                                                                                                                                          <  1:55347/251‑1 (MQ=255)                                                  gTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTCATGCTTGTGAt                                                                                                                                                                                >  1:64026/1‑251 (MQ=255)                                                         ccTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGAtggtggtg                                                                                                                                                                         <  1:70706/251‑1 (MQ=255)                                                                cAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCCAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGa                                                                                                                                                                    >  1:111288/1‑249 (MQ=255)                                                                             cTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAATTTTTGAGGTAATGCTTGAGATGGTGGAGGGGGAAGGATTATTCGtcgc                                                                                                                                                     >  1:37529/1‑251 (MQ=255)                                                                                                               cATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCG‑TTGCCTGTGGc                                                                                                                  <  1:40585/251‑1 (MQ=255)                                                                                                                      cgccgttttggCAGAATAGTTCATTGCGTAGTCGTTGGTAATTTCAGGGACTTTGGGAAGTGAGTGTGGTGGGGGAAGGAGTCGATCCTTCGAAGTCGAGGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACTAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCG‑TTGCCTGTGGCAACGTtc                                                                                                           <  2:5487/242‑1 (MQ=25)                                                                                                                          gCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCG‑TTGCCTGTGGCAACGTtctctc                                                                                                       >  1:48355/1‑251 (MQ=255)                                                                                                                           cATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGCAATGCTTGAGATGGGGGAGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCATCGGGCAG‑GTTCCTGTGGCAACGTTCTCTCg                                                                                                      >  1:62220/1‑251 (MQ=37)                                                                                                                                                           ggAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTAAAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTTATTGATTTTTTCAGGCCTTGCTTGAGATCGTGGTGGGGGAAGGATTATTTGTAGCGTCGCTCCTCACCCTTCGGGCCGGTTGA‑TGTGGCAACGTTCTCTCCCTTTCGTTCGAACCGAACCtgcttcgaagagt                                                                      >  2:111031/1‑239 (MQ=18)                                                                                                                                                                                      ggtggtggGGGCAGGATTCGAGCCTTCGAAGTCGATGACGGCAGAGTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCG‑TTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACt                                           <  2:183709/251‑1 (MQ=255)                                                                                                                                                                                                 aaGGAGTCGAACCTGCGAAGTCGAGGACGGAAGATTTAGAGCCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAAGTCAAATTTTGAGGTAGTGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCG‑TTGCCTGTGGCAACGTTCTCTCGCTTTCGCGCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAAtctc                                <  2:71725/251‑1 (MQ=255)                                                                                                                                                                                                               tCGAAGTCGATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCTCCCTTCGGGCCG‑TTTCCTGTGACTACGTTCTCTCGCTTTCGCTCGAATCGAACCTTTGTCGAGGATTCTCACCCTTCCCCGATGGGTTCAAACTTTAACAATCTCTCCGCAGTTCCca                   >  2:178234/1‑250 (MQ=37)                                                                                                                                                                                                                      ggATGACGGCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCG‑TTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCtg           <  2:48355/250‑1 (MQ=255)                                                                                                                                                                                                                              gCAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAGTGCTTGAGATGGTGGTGGGGGACGGATTATTCGTCGCTTCGCGCCTCACCCTTCGGGCCG‑TTGCCTGTGGCAACCTTCTCTCGCTTTCGCTCGAATCGAACCTTATTCGAAGGTCCTCACCCTTCCCCGATGCGTTCAAACTTTCACATTCTCACCCCAGTTACCACATCGCTGTGGTGAAt   >  2:172530/1‑251 (MQ=255)                                                                                                                                                                                                                               cAGATTTAGAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCG‑TTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAAtt  <  1:60910/251‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    GATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTAATTGCGGATTCGTTGGGAAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCAGGGGTAATTCAAATTTTGAGGTAATGCTTGAGATGGTGGTGGGGGAAGGATTATTCGTCGCTTCGCTCCTCACCCTTCGGGCCG‑TTGCCTGTGGCAACGTTCTCTCGCTTTCGCTCGAATCGAACCTTAGTCGAAGGTTCTCACCCTTCCCCGATGAGTGCAAACTTTCACAATCTCACCGAAGTTACCACATCGCTGTGGTGAATT  >  NC_000913/1287064‑1287536

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