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	94.4%	116.4 / 0.5	36	noncoding (80/87 nt)	leuQ	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base G (17/17);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

CCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAG  >  NC_000913/4606193‑4606449                                                                                                                                   |                                                                                                                                  ccTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑gggg‑ggggga                                                                                                                       >  1:1185140/1‑138 (MQ=255)      tttGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGa                                                                                                                   >  2:607768/1‑140 (MQ=255)          gTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAAccc                                                                                                               >  1:487943/1‑140 (MQ=255)           tGCGAGGGGGGGGACTTGAACCCCCCCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAAcccc                                                                                                              <  2:474683/140‑1 (MQ=255)           gggggggggggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAAcccc                                                                                                              <  1:1001879/135‑1 (MQ=255)             cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCa                                                                                                            <  1:922043/140‑1 (MQ=255)             cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCa                                                                                                            <  1:619758/140‑1 (MQ=255)               aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATTTCAACGTGGTTTTTGGAGCGAG‑GGGG‑GGGGACTTGGAcccccccc                                                                                                          >  2:960044/1‑137 (MQ=21)               aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACg                                                                                                          >  1:927326/1‑140 (MQ=255)               aGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACg                                                                                                          <  1:702622/140‑1 (MQ=255)                    ggggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGt                                                                                                     <  1:607768/140‑1 (MQ=255)                     gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCa                                                                                                            <  1:368019/132‑1 (MQ=255)                     gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCa                                                                                                            >  2:368019/1‑132 (MQ=255)                     gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTa                                                                                                    >  2:21796/1‑140 (MQ=255)                             aaCCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACt                                                                                            <  2:1185140/140‑1 (MQ=255)                               cccccACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTaa                                                                                          <  1:292641/140‑1 (MQ=255)                                            aaGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                             <  2:487943/140‑1 (MQ=255)                                                acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                                                           >  2:914588/1‑137 (MQ=255)                                                acacTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAgcg                                                                           <  1:914588/137‑1 (MQ=255)                                                      acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                   >  2:429725/1‑140 (MQ=255)                                                            gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                             >  1:841257/1‑140 (MQ=255)                                                            gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                             >  2:215458/1‑140 (MQ=255)                                                              aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGTGAGTAATATCAACGTGGTTTTTGGTGCGA‑‑GGGGCGGGGACTTGATCCTCCGCGTCCGTCATGAAACCAACACCTGAAGCTGGCGCGTCTACCTATTCCGc                                                           >  2:66150/1‑140 (MQ=25)                                                              aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                           <  1:429725/140‑1 (MQ=255)                                                                    cgcgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                                                     >  1:1032065/1‑140 (MQ=255)                                                                      cgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                   <  1:21796/140‑1 (MQ=255)                                                                      cgTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCg                                                   <  2:841257/140‑1 (MQ=255)                                                                        tCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCa                                                 >  1:339605/1‑140 (MQ=255)                                                                         cTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCAt                                                <  1:215458/140‑1 (MQ=255)                                                                                   cGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAAtt                                      >  1:292252/1‑140 (MQ=255)                                                                                     ccACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCt                                    >  1:625230/1‑140 (MQ=255)                                                                                                     cGAGCAATATCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTAc                    <  2:927326/140‑1 (MQ=255)                                                                                                             atCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCTCGGAgg            >  1:1046618/1‑140 (MQ=255)                                                                                                             atCAACGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAgg            >  2:502831/1‑140 (MQ=255)                                                                                                                 aCGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCa        <  1:960044/140‑1 (MQ=255)                                                                                                                  cGTGGTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAt       >  1:1268849/1‑140 (MQ=255)                                                                                                                    tGGTTTTTGGTGCGAGGGGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCAtt      >  1:872972/1‑140 (MQ=255)                                                                                                                      gTTTTTGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa   >  1:979457/1‑140 (MQ=255)                                                                                                                       tttttGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg  <  2:339605/140‑1 (MQ=255)                                                                                                                       tttttGGTGCGAG‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAg  <  1:66150/140‑1 (MQ=255)                                                                                                                                   |                                                                                                                                  CCTCGTTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGA‑‑GGGG‑GGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAG  >  NC_000913/4606193‑4606449

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