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	92.7%	126.2 / 3.8	41	noncoding (80/87 nt)	leuP	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base . (1/2);  new base G (17/21);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.29e-01

TAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCC  >  NC_000913/4606069‑4606341                                                                                                                                             |                                                                                                                                        tAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑g                                                                                                                                       <  2:256539/140‑1 (MQ=255) tAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑g                                                                                                                                       <  1:771733/140‑1 (MQ=255)        aaTTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑gggggggg                                                                                                                                >  2:55536/1‑140 (MQ=255)           tGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑gggggggggac                                                                                                                             >  2:361990/1‑137 (MQ=255)               gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGa                                                                                                                          <  1:361990/140‑1 (MQ=255)               gCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑ggggggggggacttg                                                                                                                         >  1:380963/1‑133 (MQ=255)                cGAGGGGGGGGACTTGACCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑gggggggggactt                                                                                                                           <  2:733240/137‑6 (MQ=255)                cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGaa                                                                                                                         <  2:203595/140‑1 (MQ=255)                cGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑gggggggggactt                                                                                                                           >  1:733240/1‑132 (MQ=255)                     ggggggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAAccccc                                                                                                                    <  2:1035431/140‑1 (MQ=255)                         ggACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTTCCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑GGGGGGGGACTTGAACCCCCACGTc                                                                                                               <  2:837532/140‑1 (MQ=255)                            cTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAACTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCg                                                                                                             >  1:413411/1‑140 (MQ=255)                             ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAATCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGt                                                                                                            <  1:921709/140‑1 (MQ=255)                             ttGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGt                                                                                                            >  1:206980/1‑140 (MQ=255)                             ttGAACCCCCACGTCCGTAAGGACACAAACACTTGAAGCTCGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGt                                                                                                            <  2:922560/140‑1 (MQ=37)                                aaCCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAACTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                         <  2:413411/140‑1 (MQ=255)                                    cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                           >  2:973781/1‑134 (MQ=255)                                    cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTa                                                                                                           <  1:973781/134‑1 (MQ=255)                                    cccACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAaca                                                                                                     >  1:760498/1‑140 (MQ=255)                                       aCGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTa                                                                                                  >  2:104911/1‑140 (MQ=255)                                        cGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTaa                                                                                                 <  2:657597/140‑1 (MQ=255)                                                     acTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAg                                                                                    >  1:983116/1‑140 (MQ=255)                                                      cTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAA‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgc                                                                                   <  1:104911/140‑1 (MQ=255)                                                       tAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAgcg                                                                                  <  2:760498/140‑1 (MQ=255)                                                             cTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTa                                                                            >  1:132327/1‑140 (MQ=255)                                                                 aGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCaa                                                                        <  2:983116/140‑1 (MQ=255)                                                                       cgcgTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGc                                                                  <  2:132327/140‑1 (MQ=255)                                                                           tCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGCCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                              >  2:235773/1‑140 (MQ=255)                                                                                  aTTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGccc                                                       >  2:25994/1‑140 (MQ=255)                                                                                    tCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                                                     >  1:633164/1‑140 (MQ=255)                                                                                     ccGCCACCTTCGCCCAGTCATTCTGATTTTTTTAAGATCCACACGTTTGGTGCG‑AGGGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAg                                                     <  2:380963/140‑1 (MQ=17)                                                                                     ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                    <  1:235815/140‑1 (MQ=255)                                                                                     ccGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGt                                                    >  2:235815/1‑140 (MQ=255)                                                                                      cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                              <  1:679719/129‑1 (MQ=255)                                                                                      cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAcc                                                              >  2:679719/1‑129 (MQ=255)                                                                                      cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                   >  2:359507/1‑140 (MQ=255)                                                                                      cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                   >  1:383345/1‑140 (MQ=255)                                                                                      cGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                   <  1:25994/140‑1 (MQ=255)                                                                                      cGCCACCTTCACACAGACAACTTACTTTTTTTGATATCGCATCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTg                                                   <  2:861077/140‑1 (MQ=18)                                                                                        ccACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCg                                                 <  1:55536/140‑1 (MQ=255)                                                                                          gccTTCGAACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAg                                               <  2:595118/139‑1 (MQ=255)                                                                                                      tCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGt                                   >  2:856341/1‑140 (MQ=255)                                                                                                         tCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGt                                >  2:398721/1‑140 (MQ=255)                                                                                                               actttttGAAATAGCGACGATTGGTGCG‑A‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGt                         <  2:734342/138‑1 (MQ=18)                                                                                                                        tatCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggg                 <  1:359507/140‑1 (MQ=255)                                                                                                                          tCGCCTCGTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAgggggg               <  1:398721/140‑1 (MQ=255)                                                                                                                                 gTTTGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTg        <  1:235773/140‑1 (MQ=255)                                                                                                                                  tttGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGa       <  2:206980/140‑1 (MQ=255)                                                                                                                                  tttGGTGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGa       >  1:809930/1‑140 (MQ=255)                                                                                                                                       tGCG‑AG‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGGCCTGAAaccc  >  1:207771/1‑140 (MQ=21)                                                                                                                                             |                                                                                                                                        TAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCG‑A‑‑GGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCCCAGTGCGAGCAATATCAACGTGGTTTTTGGTGCGAGGGGGGGGACTTGAACCCC  >  NC_000913/4606069‑4606341

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