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	86.5%	107.6 / 11.0	37	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base G (16/16);  total (18/19)
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

ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑GGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACT  >  NC_000913/4606074‑4606333                                                                                                                                        |                                                                                                                                 aTAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑gggggg                                                                                                                            >  1:272990/1‑140 (MQ=255)   aaTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑gggggggg                                                                                                                          >  1:907952/1‑140 (MQ=255)    aTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑ggggggggg                                                                                                                         >  2:910793/1‑139 (MQ=255)        gTGCGAGGGGGGGGACTTGAACCCCCACTTTCGTAAGGGCACTAACCCCTGAAGCTAGCGCGGCTACCAATTCCGCCGCCTTCGCACAGTCATCTTACTTTTTTTGATAACGCCTCGTTTGGTGCg‑g‑gggggggggactt                                                                                                                     <  2:907952/140‑6 (MQ=21)        gTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑gggggggggactt                                                                                                                     <  1:582641/140‑6 (MQ=255)          gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGa                                                                                                                  <  2:816926/140‑1 (MQ=255)          gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGa                                                                                                                  >  2:761826/1‑140 (MQ=255)          gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGa                                                                                                                  <  1:647547/140‑1 (MQ=255)           cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGaa                                                                                                                 <  1:31207/140‑1 (MQ=255)               gggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAAcccc                                                                                                             <  1:761826/140‑1 (MQ=255)                 gggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCa                                                                                                           >  2:320409/1‑140 (MQ=255)                  ggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCAc                                                                                                          <  1:320409/140‑1 (MQ=255)                        ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGt                                                                                                    >  1:882991/1‑140 (MQ=255)                          gAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGCG‑GGGGGGGG‑A‑CTTGAACCCCCACGTCCGTaa                                                                                                  >  2:652108/1‑140 (MQ=255)                               cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTc                                                                                                       <  2:519465/130‑1 (MQ=255)                               cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTc                                                                                                       >  1:519465/1‑130 (MQ=255)                                 cACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACt                                                                                           <  1:544228/140‑1 (MQ=255)                                  aCGTCCGTAAGGACACTAACACCTGAAGCTATCTCGCCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑GGGGGGGG‑AGCTTTAACCCCCACGTCCGTAAGAAGACTa                                                                                          >  2:155499/1‑140 (MQ=37)                                   cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCTCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                         <  1:155499/140‑1 (MQ=255)                                                acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                            <  1:241689/140‑1 (MQ=255)                                                  tAACCCCTGAACCTAGTGTGTCTCCCAATTCGGCAGCCTTCGCACACTTATCTTACTTTTTTTGATAACGACTCGTTTAGTGCG‑ACGGGGGGGG‑A‑CTTGAACACCCACCTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                          <  2:146366/140‑1 (MQ=18)                                                          gAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCTTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCTCGTCTAcc                                                                  >  1:146366/1‑140 (MQ=255)                                                                aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                            >  1:725218/1‑140 (MQ=255)                                                                aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                            >  2:528193/1‑140 (MQ=255)                                                                aGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CGTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTcc                                                            >  1:768702/1‑140 (MQ=255)                                                                  cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                          <  1:652108/140‑1 (MQ=255)                                                                  cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                          <  2:882991/140‑1 (MQ=255)                                                                  cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑GGGGGGGGTA‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                          <  1:172556/140‑1 (MQ=255)                                                                     gTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                       >  2:182997/1‑140 (MQ=255)                                                                                 cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                           >  2:583421/1‑140 (MQ=255)                                                                                 cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CGTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                           <  2:768702/140‑1 (MQ=255)                                                                                 cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑GGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGc                                          >  2:134371/1‑140 (MQ=255)                                                                                   ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                         >  2:814457/1‑140 (MQ=255)                                                                                   ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑GGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                               >  2:703850/1‑133 (MQ=255)                                                                                   ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑GGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                               <  1:703850/133‑1 (MQ=255)                                                                                     aCCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                       <  1:814457/140‑1 (MQ=255)                                                                                                 tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa                                    >  2:122395/1‑131 (MQ=255)                                                                                                 tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCaa                                    <  1:122395/131‑1 (MQ=255)                                                                                                    tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                        >  1:660855/1‑140 (MQ=255)                                                                                                             ttttGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGc               >  1:804059/1‑140 (MQ=255)                                                                                                              tttGATATCGCCTCGTTTGGTGCG‑A‑GGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa             <  1:134371/140‑1 (MQ=255)                                                                                                               ttGATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGa             <  2:660855/140‑1 (MQ=255)                                                                                                                 gATATCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgg           <  1:528193/140‑1 (MQ=255)                                                                                                                   tatCGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg         >  1:693906/1‑140 (MQ=255)                                                                                                                      cGCCTCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGAGCGAGCAATATCACCGTGGTTTTTGGTGCGAggggggg      >  2:132429/1‑140 (MQ=255)                                                                                                                          tCGTTTGGTGCG‑AGGGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAAAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggggggac  <  2:725218/140‑4 (MQ=255)                                                                                                                                        |                                                                                                                                 ATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑GGGGGGGG‑A‑CTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACT  >  NC_000913/4606074‑4606333

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