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	95.5%	70.9 / ‑2.8	22	noncoding (80/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base G (11/10);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.86e-01

GTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTT  >  NC_000913/4606197‑4606454                                                                                                                                |                                                                                                                                      gTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTg                                                                                                                         <  1:123650/140‑1 (MQ=255)       tGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAAcccc                                                                                                                   >  1:376134/1‑140 (MQ=255)         cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCa                                                                                                                 <  1:380684/140‑1 (MQ=255)                ggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGt                                                                                                          <  2:258207/140‑1 (MQ=255)                 gggATTTGAACCCCCGCCACCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                         <  2:29155/140‑1 (MQ=25)                 gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                         >  2:276888/1‑140 (MQ=255)                    aCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAgg                                                                                                      <  1:276888/140‑1 (MQ=255)                         aaCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                                 <  2:376134/140‑1 (MQ=255)                                                     ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAt                                                                     >  1:524300/1‑140 (MQ=255)                                                        gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                  >  2:343566/1‑140 (MQ=255)                                                              aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                            <  2:524300/140‑1 (MQ=255)                                                                cgcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                                                          >  1:337918/1‑140 (MQ=255)                                                                    tCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCa                                                      >  2:441635/1‑140 (MQ=255)                                                                               cGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAt                                            >  2:163469/1‑140 (MQ=255)                                                                                                        tataAACGTGGTTTTTGGGGGCGG‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg                  <  1:496876/140‑1 (MQ=21)                                                                                                         atCAACGTGGTTTTT‑GGTGCGAGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg                  <  1:163469/140‑1 (MQ=255)                                                                                                                tGGTTTTT‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTc          >  1:4387/1‑140 (MQ=255)                                                                                                                     ttt‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcgc              <  2:143679/131‑1 (MQ=255)                                                                                                                     ttt‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcgc              >  1:143679/1‑131 (MQ=255)                                                                                                                      tt‑GGTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTgg    >  1:340346/1‑140 (MQ=255)                                                                                                                       t‑GGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCTCAATTCTTTAAAAGAAATGCTACCACGGAGGCGCATTCTAgctgtc  >  2:496876/1‑135 (MQ=17)                                                                                                                         ggTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg               <  1:523648/127‑1 (MQ=255)                                                                                                                         ggTGCGAG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg               >  2:523648/1‑127 (MQ=255)                                                                                                                                |                                                                                                                                      GTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTT‑GGTGCGA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTT  >  NC_000913/4606197‑4606454

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