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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,606,207	1	.	G	92.9%	40.1 / ‑2.7	14	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (5/8);  total (5/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAATTGGTGCGA‑GGGGG‑GG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGG  >  NC_000913/4606075‑4606328                                                                                                                                           |                                                                                                                            tAATTGGTGCGAGGGGGG‑GG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑ggg‑ggg                                                                                                                     >  1:80017/1‑140 (MQ=255) aaaTTGGTGGGA‑GGGGG‑GGTGACTTGCACCCCGCCGGCCGTAAATGCACTAAATCCTGAAGATAGTGCGTCTACCAAATCCGCCATCTTCGAACAATGATCTAACAATTTTTGAAATCACCTCGTTTGGTGC‑G‑A‑‑ggg‑ggg                                                                                                                     <  2:122385/139‑1 (MQ=255)  aaTTGGTGCGAGGGGGG‑GG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTCGATCTCCCCTCGTTTGGTGC‑G‑A‑‑ggg‑gggg                                                                                                                    >  1:240630/1‑140 (MQ=18)   aTTGGTGCGA‑GGGGGTGG‑GACTAGACCCCCCACGTCCCTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑ggg‑ggggg                                                                                                                   <  1:192262/140‑1 (MQ=21)          cgagggggg‑gg‑gACTAGAACCCACATGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGa                                                                                                             <  1:281467/136‑1 (MQ=21)                       gagttcgaCCCTCGCGTCCGTACGGTCACTACAACTTGATGCTAGCGCGGCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTATTTGGATATCGCCTTGTTTGGTGCGG‑G‑‑GGG‑GGGGGACTTGAACCCCCACGTc                                                                                                  <  2:322653/133‑1 (MQ=21)                             aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAg                                                                                            <  2:126856/140‑1 (MQ=255)                                ccccACGTCCGTAAGGATACTATCACCTGAATCTAGCGCTTCTACCAAGTCCGCCGCCTTCGCACAGTGATCTTACTTTTTTTGATGTCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAac                                                                                         <  2:240630/140‑1 (MQ=37)                                     cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                    <  1:241686/140‑1 (MQ=255)                                     cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                    <  2:468283/140‑1 (MQ=255)                                                  acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGGGC‑G‑AG‑GGG‑GGGGGACTTGAACTCCCGCGTCCGTATGAACCCTAACCTCTGAAGCttg                                                                       >  2:288147/1‑138 (MQ=37)                                                  acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATATTACTTTTTATGATATCCTCTCGTTTGGCGC‑GAA‑‑GGC‑ATCGGACTTGAACGCCCGCATCCAGACGAACCCTAACACCTGAAGCTac                                                                       >  2:234486/1‑139 (MQ=16)                                                        acCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGAATTGAACCCCCACGTCCGTAAGAACAGTAACTCCTGAAGCTAGCGCGTc                                                                 >  2:222240/1‑140 (MQ=37)                                                         ccTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCt                                                                >  1:144938/1‑140 (MQ=255)                                                                    cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                     <  1:49297/140‑1 (MQ=255)                                                                        tCTACCAATTCCGCCACCTTCGCACAGTCCTCTTACTTTTTTTGATATCGCCTCGTCTGGTGC‑G‑A‑‑GGGCGGGGGCCTTTAACCCCCAGGTCCGTAAGCACAGTAGCACCGGAAGCTAGCGCGTCTACCCATTCCGCCAcc                                                 >  2:51038/1‑140 (MQ=21)                                                                                     ccACCTTCGCACAGTCACCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTAAACCCCCACGTCCGTATGACCACTAACACCTGAATCTAGCGCGTCTACCAAGTCCGGCACCTTCGCCCacggcg                                    >  2:271692/1‑135 (MQ=21)                                                                                            cGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCActa                             >  2:58343/1‑140 (MQ=255)                                                                                              cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATc                           <  2:144938/140‑1 (MQ=255)                                                                                                         tACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACGCCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttt                <  1:51038/140‑1 (MQ=255)                                                                                                             ttttttGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGtttttcgg            >  2:619158/1‑137 (MQ=255)                                                                                                             ttttttGATATCGCCTCGTTTGGTGC‑G‑AG‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt            >  1:406165/1‑140 (MQ=255)                                                                                                                        cGCCTCATTTGGTGG‑G‑GGGGGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCCTTTTCCAATTCCGCCACCTACGCCCATTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg  <  2:80017/140‑1 (MQ=18)                                                                                                                                           |                                                                                                                            TAATTGGTGCGA‑GGGGG‑GG‑GACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGC‑G‑A‑‑GGG‑GGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGG  >  NC_000913/4606075‑4606328

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