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	95.6%	139.7 / 0.3	45	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (23/20);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.93e-01

AGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCC  >  NC_000913/4606067‑4606340                                                                                                                                              |                                                                                                                                         aGTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGa                                                                                                                                           <  2:392753/140‑1 (MQ=255)  gTAAAGATAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑g                                                                                                                                        <  1:642189/140‑1 (MQ=255)         tAATTGGTGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑gggg‑gtgg                                                                                                                                >  1:80195/1‑140 (MQ=255)                tGCGA‑GGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGa                                                                                                                          >  1:722360/1‑140 (MQ=255)                  gcgagggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGaa                                                                                                                         <  2:412426/136‑1 (MQ=255)                   cgagggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAAc                                                                                                                        >  1:250240/4‑140 (MQ=255)                     agggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAAccc                                                                                                                      >  2:817142/2‑140 (MQ=255)                       ggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAAccccc                                                                                                                    >  1:674909/1‑140 (MQ=255)                               ttGAACCCCCGCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGt                                                                                                            <  1:685497/140‑1 (MQ=255)                               ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGt                                                                                                            <  2:674909/140‑1 (MQ=255)                               ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGt                                                                                                            <  1:658663/140‑1 (MQ=255)                                  aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAg                                                                                                         >  1:548621/1‑140 (MQ=255)                                  aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAg                                                                                                         <  1:50811/140‑1 (MQ=255)                                       ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAacac                                                                                                    >  2:489004/1‑140 (MQ=255)                                        cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                                   <  1:782801/140‑1 (MQ=255)                                        cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACt                                                                                                   >  1:187986/1‑140 (MQ=255)                                          cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                 <  2:250240/140‑1 (MQ=255)                                          cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                 <  2:722360/140‑1 (MQ=255)                                            tCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAaca                                                                                               >  2:991043/1‑140 (MQ=255)                                                tAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGATTTGAACCCCCACGCCCGAAATAACACTAACACCtt                                                                                           >  2:1076506/1‑139 (MQ=37)                                                        cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                                                                   <  2:187986/140‑1 (MQ=255)                                                            cacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGt                                                                               >  2:435226/1‑140 (MQ=255)                                                               cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                                                            >  2:264734/1‑140 (MQ=255)                                                                 gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                          >  1:247930/1‑140 (MQ=255)                                                                   aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                        <  1:435226/140‑1 (MQ=255)                                                                     cTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                                                      <  1:264734/140‑1 (MQ=255)                                                                       aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGCCTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACAACTGAAGCTAGCGCACGTAGCAATTcc                                                                    >  2:351545/1‑140 (MQ=25)                                                                       aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTATTATTACTTTAAATGATATCCTCCCGTTTAGGGCGAG‑GGGG‑GGGGACTTCAACTCCCACGCACTAAATAACACTCGCACCTGGAGCAATTGCTGCTACCATTTcc                                                                    >  2:131753/1‑140 (MQ=255)                                                                         cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTccc                                                                   >  1:607668/1‑139 (MQ=255)                                                                         cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCg                                                                   >  1:56793/1‑140 (MQ=255)                                                                         cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  1:991043/140‑1 (MQ=255)                                                                         cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  1:489004/140‑1 (MQ=255)                                                                             tCTACCAATTCCGCCACCTTCGCACAGTCATCTTAC‑TTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCt                                                             >  1:601119/1‑140 (MQ=255)                                                                                  cAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGc                                                         <  2:548621/140‑1 (MQ=255)                                                                                   aaTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  2:247930/140‑1 (MQ=255)                                                                                   aaTACCGCCGCATTCGAACAGTCATCTTACTTTAGTTGATATCGAAACGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  2:436917/140‑1 (MQ=17)                                                                                      tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                                                     >  2:208667/1‑140 (MQ=255)                                                                                        cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        <  1:755479/136‑1 (MQ=255)                                                                                        cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                        >  2:755479/1‑136 (MQ=255)                                                                                           cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGCCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                   >  1:950965/1‑137 (MQ=255)                                                                                           cACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGGCGGGGACTTGAACCCCCACGCCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                   <  2:950965/137‑1 (MQ=255)                                                                                                   cacaGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATc                                        <  1:817142/140‑1 (MQ=255)                                                                                                        tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACga                                   >  2:661220/1‑139 (MQ=255)                                                                                                           tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                     <  1:1067297/135‑1 (MQ=255)                                                                                                           tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAAc                                     >  2:1067297/1‑135 (MQ=255)                                                                                                           tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGCGTGAACAATATAAACGAGGt                                >  2:702046/1‑140 (MQ=18)                                                                                                                cTTTTTTTGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                           >  2:618773/1‑140 (MQ=255)                                                                                                                     tttGATATCGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                      <  1:618773/140‑1 (MQ=255)                                                                                                                             cGCCTCGTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCCATATCAACGTGGTTTTTGGTGCGA‑ggggggg             >  2:203070/1‑140 (MQ=255)                                                                                                                                        gTGCGA‑‑GGGGTGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAAcc   >  1:367438/1‑140 (MQ=255)                                                                                                                                         tGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGACTTGAAccc  >  2:308510/1‑140 (MQ=255)                                                                                                                                              |                                                                                                                                         AGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGGGGGGGACTTGAACCC  >  NC_000913/4606067‑4606340

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: ‑