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

TGATATCGCCTCGTTTGGTGCGAGG‑GGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTT  >  NC_000913/4606185‑4606455                                                                                                                                              |                                                                                                                                      tGATATCGCCTCGTTTGGTGCGCGGAGGGCGGGACTTTAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑a                                                                                                                                        <  1:208057/140‑1 (MQ=21)                     ggggg‑ggg‑gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCa                                                                                                                  <  2:64791/137‑1 (MQ=255)                               gggACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                          >  1:364420/1‑140 (MQ=255)                                 gACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                        <  1:370374/140‑1 (MQ=255)                                    ttGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAAccccc                                                                                                                   <  2:421638/126‑1 (MQ=255)                                    ttGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAAccccc                                                                                                                   >  1:421638/1‑126 (MQ=255)                                              aCGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                                                                           >  2:557429/1‑140 (MQ=255)                                                     tAAGAACGCTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTa                                                                                    <  2:420718/140‑1 (MQ=255)                                                            acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTACGTAAGGACAGTAACAACTGAAGCTAGCGCGTc                                                                             >  2:318332/1‑140 (MQ=37)                                                                acacCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                         >  2:431736/1‑140 (MQ=255)                                                                      gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAACGACACTAGCACATGAAGCTAGGGCGTCTACCAATTcc                                                                   >  1:214435/1‑140 (MQ=25)                                                                      gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTCGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                                   >  1:680436/1‑140 (MQ=255)                                                                          cTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                               <  1:557429/140‑1 (MQ=255)                                                                            aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGGCTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                               <  2:108351/138‑1 (MQ=255)                                                                            aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                               >  2:670167/1‑138 (MQ=255)                                                                            aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                               <  1:670167/138‑1 (MQ=255)                                                                            aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                             >  1:415523/1‑140 (MQ=255)                                                                            aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCa                                                               >  1:108351/1‑138 (MQ=255)                                                                            aGCGCGTCTACCAATTCCGCCACCTTCGCCCAGAGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGCCTTGAACCCCCACGTCCGTAAGGACACTAATACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                             >  2:535589/1‑140 (MQ=25)                                                                                   cTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCAt                                                      <  1:431736/140‑1 (MQ=255)                                                                                            ccGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCACCTGCCGGAAGAAAACTAACACCTGAAGCGAGCGCGTCTCCCACTTCCGCCACCTTCCGATACCATCatt                                             >  2:471168/1‑138 (MQ=18)                                                                                                 cccTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTa                                        <  1:619399/139‑1 (MQ=255)                                                                                                       gCCCAGTGCGAACAACCTCAATTTGTTTTTTGGTGCGCG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGa                                  <  1:231490/140‑1 (MQ=18)                                                                                                         ccAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAAt                                <  2:415523/140‑1 (MQ=255)                                                                                                            gTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGc                             <  1:318332/140‑1 (MQ=255)                                                                                                                      tatCAACGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAg                   >  1:685422/1‑140 (MQ=255)                                                                                                                            cGTGGTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCTTCAATTCTTAAAAAGAATTGCCACCACGGAGGCGCAt             >  1:328551/1‑140 (MQ=255)                                                                                                                                gTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCACTTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTa         >  1:254835/1‑140 (MQ=255)                                                                                                                                gTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg                <  1:656113/133‑1 (MQ=255)                                                                                                                                gTTTTTGGTGCG‑AGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGgcg                >  2:656113/1‑133 (MQ=255)                                                                                                                                      ggTGCG‑A‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGttt  >  2:568466/1‑140 (MQ=255)                                                                                                                                              |                                                                                                                                      TGATATCGCCTCGTTTGGTGCGAGG‑GGG‑GGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCG‑A‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCATACCATCAATTCTTAAAAAGAATTGCTACCACGGAGGCGCATTCTAGTGGTTT  >  NC_000913/4606185‑4606455

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