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	78.1%	75.6 / 18.5	32	noncoding (80/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (3/4);  new base G (15/10);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCG  >  NC_000913/4606207‑4606441                                                                                                                     |                                                                                                                          aGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCAc                                                                                                   >  1:317445/1‑139 (MQ=255)       gggCCTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCGAGCGCGTCTACCAGTTCCGCCACCTTCGCCCAGTGCGTGAAATATCAACCTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGt                                                                                             <  1:450632/139‑1 (MQ=21)       gggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTACGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GTCTTGAACCTCAACTTCCgg                                                                                             >  2:190038/1‑138 (MQ=18)              gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCa                                                                                                    >  2:227568/1‑125 (MQ=255)              gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCa                                                                                                    <  1:227568/125‑1 (MQ=255)              gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGt                                                                                             <  2:298656/132‑1 (MQ=255)              gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGt                                                                                             >  1:298656/1‑132 (MQ=255)                 cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGAGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTa                                                                                   >  2:986674/1‑139 (MQ=255)                         tCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTGTTTGGTGCGAG‑GGGGGGG‑GCCTTGAACCCCCCCGTCCGGAAGGACACTaataacaga                                                                           >  1:487818/1‑132 (MQ=25)                                    acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                                               >  2:932535/1‑139 (MQ=255)                                        acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc                                                           <  1:932535/139‑1 (MQ=255)                                            cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                                        >  1:418371/1‑139 (MQ=255)                                             tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGG‑GACTTGAACCCCCa                                                                                                    >  2:867439/1‑93 (MQ=255)                                             tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGG‑GACTTGAACCCCCa                                                                                                    <  1:867439/93‑1 (MQ=255)                                              gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTa                                                                                            >  1:200127/1‑101 (MQ=255)                                              gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTa                                                                                            <  2:200127/101‑1 (MQ=255)                                                aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                    <  1:986674/139‑1 (MQ=255)                                                    aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                 <  2:169882/122‑1 (MQ=255)                                                    aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                 >  1:169882/1‑122 (MQ=255)                                                      cgcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                                              >  2:374958/1‑139 (MQ=255)                                                         gTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGGAGACTTAAACCCCCACATCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                           <  2:418371/139‑1 (MQ=255)                                                          tCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                          >  2:78925/1‑139 (MQ=255)                                                                                     tGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATa                                       <  2:203510/115‑1 (MQ=255)                                                                                     tGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATa                                       >  1:203510/1‑115 (MQ=255)                                                                                     tGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGc               >  1:939162/1‑139 (MQ=255)                                                                                         aGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGG‑GACTTGAACCCCCACGTCCGt                                                                                             >  1:477323/1‑56 (MQ=255)                                                                                         aGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGG‑GACTTGAACCCCCACGTCCGt                                                                                             <  2:477323/56‑1 (MQ=255)                                                                                              tatCAACGTGGTTTTTGGTGCGAGGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACgg       >  2:260387/1‑139 (MQ=255)                                                                                                  agcGTGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg  <  1:369691/137‑1 (MQ=255)                                                                                                   aCGTGGTTTTTGGTGCGAGGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg  <  1:260387/139‑1 (MQ=255)                                                                                                      tGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCATAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg  >  2:1025062/1‑135 (MQ=255)                                                                                                      tGGTTTTTGGTGCGAG‑GGGGGGG‑GACTTGAACCCCCACGTCCATAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg  <  1:1025062/135‑1 (MQ=255)                                                                                                                     |                                                                                                                          AGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCG  >  NC_000913/4606207‑4606441

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