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	96.3%	77.1 / ‑3.1	27	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (12/14);  total (12/15)
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.03e-01

TAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCC  >  NC_000913/4606069‑4606340                                                                                                                                             |                                                                                                                                       tAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑g                                                                                                                                       <  2:504079/140‑1 (MQ=255)       tAATTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCGTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑gggggg                                                                                                                                  >  2:512905/1‑140 (MQ=255)         aTTGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑ggggggg‑g                                                                                                                               <  2:78720/140‑1 (MQ=255)          ttGGTGCGAGGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑ggggggg‑gg                                                                                                                              <  1:416484/140‑2 (MQ=255)                    ggggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGG‑GACTTGAAcccc                                                                                                                    <  2:688687/140‑1 (MQ=255)                         gggcattGAAAACCCCCCTCCGTAAGGACACCAACACCTTAAGCTAGCTCGTCTACCAATTCAGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATAGCCACGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACg                                                                                                                <  2:318200/135‑1 (MQ=25)                             cTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTCCTCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGGAGCCTCGAACCCCCACGGccc                                                                                                            >  2:206799/1‑139 (MQ=21)                              ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGt                                                                                                           <  1:91265/140‑1 (MQ=255)                                      ccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAacac                                                                                                   >  1:9862/1‑140 (MQ=255)                                       cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCACCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACt                                                                                                  <  1:143610/140‑1 (MQ=255)                                         cGTCCGTAAGGACACTAACACCAGAAGCTCGCGCGTCAACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                <  1:323341/140‑1 (MQ=38)                                             cGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                                                            >  2:637018/1‑140 (MQ=255)                                                                gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                 >  1:195129/1‑132 (MQ=255)                                                                gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                 <  2:195129/132‑1 (MQ=255)                                                                gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTAcc                                                                         >  2:822339/1‑140 (MQ=255)                                                                 aaGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTTGTTTGGTTCGAGGGGGGGG‑GGCTTGATCCCCCGCGCCCGGAAGTACACTAACTCCTAGTGCTAGCGCGTCTACCa                                                                        >  2:365166/1‑140 (MQ=25)                                                                        cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                 <  1:637018/140‑1 (MQ=255)                                                                            tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                             >  1:128949/1‑140 (MQ=255)                                                                            tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                             >  1:22873/1‑140 (MQ=255)                                                                                  aaTTCCGCCACCAACCCACAGTAATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcc                                                       <  1:271316/140‑1 (MQ=18)                                                                                      ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                   <  2:9862/140‑1 (MQ=255)                                                                                       cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                  <  2:22873/140‑1 (MQ=255)                                                                                       cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                  >  1:734724/1‑140 (MQ=255)                                                                                          cACCTTCGCACAGTCTTCTTACCTTTTTTCATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGa                                               <  1:365166/140‑1 (MQ=21)                                                                                                       tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                  >  2:595684/1‑140 (MQ=255)                                                                                                               cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                          >  2:528582/1‑140 (MQ=255)                                                                                                               cTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTg                          >  1:580630/1‑140 (MQ=255)                                                                                                                    tttGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                     <  1:822339/140‑1 (MQ=255)                                                                                                                    tttGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCg                     <  2:734724/140‑1 (MQ=255)                                                                                                                     ttGATATCGCCTCGTTTGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGAGAGCAATATCAACGTGGTTTTTGGTGCGa                    >  1:218468/1‑140 (MQ=255)                                                                                                                                   tttGGTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGa      <  1:528582/140‑1 (MQ=255)                                                                                                                                       gTGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAAccc  >  1:556004/1‑140 (MQ=255)                                                                                                                                        tGCGAGGGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACAcc                                                                                            <  2:471983/49‑1 (MQ=11)                                                                                                                                             |                                                                                                                                       TAAAGATAATTGGTGCGA‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGA‑GGGGGGG‑GACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCC  >  NC_000913/4606069‑4606340

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