Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,818,158	T→A	13.7%	intergenic (‑23/+40)	argZ ← / ← argY	tRNA‑Arg/tRNA‑Arg

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,818,158	0	T	A	13.7%	87.2 / 5.1	35	intergenic (‑23/+40)	argZ/argY	tRNA‑Arg/tRNA‑Arg
Reads supporting (aligned to +/- strand):  ref base T (10/20);  new base A (2/3);  total (12/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 = 9.91e-01

TTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTA  >  NC_000913/2818012‑2818292                                                                                                                                                   |                                                                                                                                       ttCCGCTACCTCTTTGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa                                                                                                                                      <  1:347858/149‑1 (MQ=255)  tCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGa                                                                                                                                      >  5:327560/1‑148 (MQ=255)    cGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGata                                                                                                                                    <  7:109478/148‑1 (MQ=255)        aCCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGAtttttt                                                                                                                                >  2:322787/1‑146 (MQ=255)         ccTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGa                                                                                                                                <  7:168122/147‑1 (MQ=255)         ccTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGAttttttg                                                                                                                               <  2:326401/148‑4 (MQ=14)               tgttACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGAttttt                                                                                                                                 <  4:324012/137‑2 (MQ=14)               cGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGc                                                                                                                          <  8:134906/147‑1 (MQ=255)               cGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCAATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCg                                                                                                                           <  4:116006/146‑1 (MQ=255)                   aCCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGCCCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGctacta                                                                                                                     >  8:198390/1‑148 (MQ=255)                      ttCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc                                                                                                                 >  6:52512/1‑149 (MQ=255)                        cATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGcc                                                                                                                <  6:129988/148‑1 (MQ=255)                        cATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc                                                                                                                 <  5:52512/147‑1 (MQ=255)                        cATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGAtttttt                                                                                                                                <  5:313721/132‑3 (MQ=11)                        cATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTACTGCAGAtttttt                                                                                                                                >  6:313721/1‑130 (MQ=11)                             cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATc                                                                                                           <  7:127186/148‑1 (MQ=255)                             cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATc                                                                                                           <  2:369675/148‑1 (MQ=255)                             cTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATc                                                                                                           >  1:369675/1‑148 (MQ=255)                                aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCaa                                                                                                         <  8:256160/147‑1 (MQ=255)                                aaGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAg                                                                                                        <  7:164718/148‑1 (MQ=255)                                   aTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAag                                                                                                    <  5:86175/149‑1 (MQ=255)                                   aTTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGAAag                                                                                                    <  8:248845/149‑1 (MQ=255)                                    tttCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAaga                                                                                                   <  8:149116/149‑1 (MQ=255)                                         gagaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATgg                                                                                              <  7:198390/149‑1 (MQ=40)                                          agaATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGt                                                                                             <  6:327560/149‑1 (MQ=40)                                                 ggTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc                                                                                                                 <  6:373612/122‑1 (MQ=40)                                                 ggTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGc                                                                                                                 >  5:373612/1‑122 (MQ=40)                                                                       ccTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAAATTATTTTACTGCTGATATTTATAAACCTAATAAAGCCATATAAAGAAAGAGATGGTGAATACGGGAGGATTAGACGCGCCGAcgg                                                                >  6:190002/1‑147 (MQ=11)                                                                                    cgctCGTAGCCGAGTACTCTTTCAAGATCAGCTACTTAAGCATAGCGAAATTTATTTTACTGCTCATCTTGTAACCCGTCCTAACTCCATATCTAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTa                                                     <  5:329202/144‑1 (MQ=16)                                                                                                                       ggATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                  <  2:372719/115‑1 (MQ=255)                                                                                                                       ggATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGcc                                                  >  1:372719/1‑115 (MQ=255)                                                                                                                                gggAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGGAGTCTATCAAGCTGAGCGACGGAGGCAGCgggcgact         >  7:370285/1‑142 (MQ=25)                                                                                                                                     atttatttTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATACGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTAATCTATCCAGCTGAGCTACGGATGCATCGGGAAACttacttta  >  5:66465/1‑149 (MQ=255)                                                                                                                                        tatttTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCg                <  2:84188/132‑1 (MQ=255)                                                                                                                                        tatttTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCg                >  1:84188/1‑132 (MQ=255)                                                                                                                                                   |                                                                                                                                       TTCCGCTACCTCTTCGATACCTTCATTGCTGAAGATTTCGGAGAATATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAATTTATTTTACTGCTGATATTGATACCGCTACTAAAGCCATATCAAGTAAGAGATGGTGCATCCGGGAGGATTCGAACCTCCGACCGCTCGGTTCGTAGCCGAGTACTCTATCCAGCTGAGCTACGGATGCATCGGGAAACTTACTTTA  >  NC_000913/2818012‑2818292

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑