Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,606,330	(G)8→9	noncoding (72/87 nt)	leuQ ←	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,322	1	.	G	94.4%	52.3 / ‑2.8	18	noncoding (80/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (13/4);  total (13/5)
Fisher's exact test for biased strand distribution p-value = 2.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

CGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGT  >  NC_000913/4606191‑4606450                                                                                                                                       |                                                                                                                                 cGCCTCATTTGGTGGGGGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCCTTTTCCAATTCCGCCACCTACGCCCATTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑gggggg                                                                                                                            <  2:80017/140‑1 (MQ=18)                  agggggggggACTTGAACCCCCCCGCCCGAAAGAACAATAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACg                                                                                                           <  1:222240/138‑1 (MQ=25)                         gggACTTGAACCCCCCCGTCCGTAACAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                     <  1:288147/140‑1 (MQ=255)                               tGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGAGTTAAACCCCAACGTCCGTAAGGaca                                                                                               >  1:170845/1‑140 (MQ=25)                                                      acTAACACCTGAACCTAACGCGTCTACCAATTCCGCCGCCTTCGCCCAGTGCGAGCAATAACAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                                                        <  2:406165/140‑1 (MQ=25)                                                             ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                                                 >  1:394255/1‑140 (MQ=255)                                                                gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGTCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                              >  2:581113/1‑140 (MQ=255)                                                                gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                              >  1:342414/1‑140 (MQ=255)                                                                      aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                        >  1:437581/1‑140 (MQ=255)                                                                                               tCGCCCAGTGCGAGCTATATCACCGGGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTaaaaa                               <  1:234486/140‑1 (MQ=18)                                                                                                             cAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACg                 <  1:619158/140‑1 (MQ=255)                                                                                                               atatCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGa               >  1:491721/1‑140 (MQ=255)                                                                                                               atatCAACGCGGTTTTTGGTGCGA‑GTGGGCGGACTTGTACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg              <  1:271692/140‑1 (MQ=255)                                                                                                                 atCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg             >  1:523527/1‑140 (MQ=255)                                                                                                                 atCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg             >  1:283476/1‑140 (MQ=255)                                                                                                                      cGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAt        >  1:418414/1‑140 (MQ=255)                                                                                                                      cGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAt        >  1:29158/1‑140 (MQ=255)                                                                                                                        tGGTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTGCGCATACCGTCAATTCTTAAAAAGAATTGCAACCACGGAGGCGCATTc      >  2:257708/1‑140 (MQ=21)                                                                                                                          gTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa    >  1:633502/1‑140 (MQ=255)                                                                                                                          gTTTTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa    >  1:659881/1‑140 (MQ=255)                                                                                                                            ttttGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGt  >  2:486080/1‑140 (MQ=255)                                                                                                                                       |                                                                                                                                 CGCCTCGTTTGGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGT  >  NC_000913/4606191‑4606450

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