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	86.1%	124.4 / 14.6	43	noncoding (80/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (2/4);  new base G (17/20);  total (19/24)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCTCGTTTGGTGCGA‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAG  >  NC_000913/4606191‑4606449                                                                                                                                       |                                                                                                                                 cGCCTCGTTTGGTGCGAGGGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCCATATCAACGTGGTTTTTGGTGCGA‑gg‑ggggg                                                                                                                          >  2:203070/1‑140 (MQ=255)            gTGCGA‑GGGGTGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAAcc                                                                                                                >  1:367438/1‑140 (MQ=255)             tGCGAGGGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAAccc                                                                                                               >  2:308510/1‑140 (MQ=255)                ggggggg‑gggaACTTGAACCCCCGCGTCTGTAAGAACACTAACACCTAAAGCTAGCACGTATACCAATTCCGCCACCTTCGCCCAGTGCGAGAAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCa                                                                                                            <  2:867523/137‑1 (MQ=21)                gaggggg‑ggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCCCCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCa                                                                                                            <  1:351545/137‑1 (MQ=255)                gaggggg‑ggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCa                                                                                                            <  1:308510/137‑1 (MQ=255)                 aggggg‑ggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCAc                                                                                                           >  1:210039/3‑140 (MQ=255)                         gggACTTGAACCCCCTCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GGTGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                    <  2:56793/140‑1 (MQ=255)                         gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTa                                                                                                    <  2:244738/140‑1 (MQ=255)                         gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTa                                                                                                    <  2:367438/140‑1 (MQ=255)                           gACTTGAACCCCCTCGTCCGTAAGAACACTAAAGCCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                  <  2:607668/140‑1 (MQ=25)                             cTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGa                                                                                                >  1:681827/1‑140 (MQ=255)                                 aaCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                            <  2:210039/140‑1 (MQ=255)                                   cccccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCTCCACCTTCGCCCAGTACAAGCAATATCAACGAGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                          <  2:1081221/140‑1 (MQ=25)                                                    acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                         <  2:681827/140‑1 (MQ=255)                                                    acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                         <  1:203070/140‑1 (MQ=255)                                                      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                             >  2:807750/1‑134 (MQ=255)                                                      acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                             <  1:807750/134‑1 (MQ=255)                                                          acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCACTATCAACGTGGTTTTTGGAGCGAGGG‑GGGGGGTCTTGAACCCCCACGTCCGTAAGGTCACTAACACCTGAAGCAAGcgcgcttatc                                                                   >  2:122231/1‑134 (MQ=37)                                                               tGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                             >  1:525394/1‑140 (MQ=255)                                                                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                         >  2:980462/1‑126 (MQ=255)                                                                  aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGcgcg                                                                         <  1:980462/126‑1 (MQ=255)                                                                   gCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                         >  2:574656/1‑140 (MQ=255)                                                                   gCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                         <  1:208667/140‑1 (MQ=255)                                                                    gtAGCGCGTCTACCAATTCCGCCCGCTTCGCCCAGTTCGTGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                         <  2:601119/139‑1 (MQ=21)                                                                            tCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCa                                                 >  2:1082022/1‑140 (MQ=255)                                                                              tACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATa                                               >  1:686238/1‑140 (MQ=255)                                                                              tACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATa                                               >  1:408432/1‑140 (MQ=255)                                                                                    ttCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCaa                                        <  1:574656/140‑1 (MQ=255)                                                                                      ccGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAt                                       <  1:702046/140‑1 (MQ=255)                                                                                       cGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                        <  2:991938/122‑1 (MQ=255)                                                                                       cGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                        >  1:991938/1‑122 (MQ=255)                                                                                         ccACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCt                                    >  1:948842/1‑140 (MQ=255)                                                                                                 gCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGaa                           <  2:525394/140‑1 (MQ=255)                                                                                                   ccAGTGCGATCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCGCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAAt                          <  1:1033652/140‑1 (MQ=255)                                                                                                        gCGAGCAATATCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTa                     >  1:134964/1‑140 (MQ=255)                                                                                                                 atCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg            >  2:445985/1‑140 (MQ=255)                                                                                                                 atCAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg            >  1:191320/1‑140 (MQ=255)                                                                                                                   cAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg          >  2:160054/1‑140 (MQ=255)                                                                                                                   cAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg          <  1:218487/140‑1 (MQ=255)                                                                                                                   cAACGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg          <  1:160054/140‑1 (MQ=255)                                                                                                                    aaCGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcgc         <  2:686238/140‑1 (MQ=255)                                                                                                                    aaCGTGGTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCAGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcgc         <  2:408432/140‑1 (MQ=255)                                                                                                                          gTTTTTGGTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTTAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa   >  1:504424/1‑140 (MQ=255)                                                                                                                                ggTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCGACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg  <  2:346580/135‑1 (MQ=255)                                                                                                                                ggTGCGAGGG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCGACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg  >  1:346580/1‑135 (MQ=255)                                                                                                                                       |                                                                                                                                 CGCCTCGTTTGGTGCGA‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑GG‑GGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAG  >  NC_000913/4606191‑4606449

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