Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,947,467	A→C	12.7%	intergenic (+4/‑30)	metZ → / → metW	tRNA‑Met/tRNA‑Met

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,467	0	A	C	12.7%	76.3 / 2.0	30	intergenic (+4/‑30)	metZ/metW	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (8/18);  new base C (2/2);  total (10/20)
Fisher's exact test for biased strand distribution p-value = 5.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCC  >  NC_000913/2947321‑2947566                                                                                                                                                   |                                                                                                    aTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                                   <  7:390857/149‑1 (MQ=255)  tCCGTGATCTGGTTAAAAAAAAGGTTGCATGAAATCGATTTCGGAGTATAGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                  <  8:415630/149‑2 (MQ=18)    cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAtt                                                                                                <  6:233979/149‑1 (MQ=255)        aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                            <  3:224965/149‑1 (MQ=255)        aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                            <  2:90380/149‑1 (MQ=255)         tCGGTGTAAAAAAAAGGTTGCATGAAAACGCGAGCGTAGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                           <  2:163930/149‑1 (MQ=255)         tCGGGGTAAAAACAAGGTGGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCGGCGGGCTCATCACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                           <  1:286030/149‑1 (MQ=25)         tCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatg                                                                                           <  1:150436/149‑1 (MQ=255)            gggTAAAAAAAAGGTTGCATGAAATCTCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTAGGGCTCTTAACCCGAAGGTCGTCGGTTCAATTCCGGCCCCCGCAACCAATTAAAATTTGATGAAg                                                                                        <  7:376852/149‑1 (MQ=21)            gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAg                                                                                        <  1:116799/149‑1 (MQ=255)                  aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAACCTCGTAACACCTCCGGCTAATAAACCAAAGGTCGTCCGTTAAAATACAGCCCCCGAAACCAATTAAAAATAAAAAAAGTAAAg                                                                                   >  2:279143/1‑148 (MQ=255)                   aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                               <  8:144321/135‑1 (MQ=255)                   aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                               >  7:144321/1‑135 (MQ=255)                        ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                            <  4:284530/149‑1 (MQ=255)                                   aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAATGCAGTACGGTGACGCgggt                                                                 >  5:202676/1‑148 (MQ=255)                                    aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                <  2:421724/149‑1 (MQ=255)                                    aCGCGAGCGGAGTATAGTGCGCATCCACGGAAGCGGGGTGGAGCAGCCTTGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCTGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                <  4:229347/149‑1 (MQ=255)                                        gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgatg                                                                                           >  1:137172/1‑117 (MQ=255)                                        gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgatg                                                                                           <  2:137172/117‑1 (MQ=255)                                         cgCGGAGTATAGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCTTAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagc                                                           <  1:279143/148‑1 (MQ=255)                                         aGCGGAGTATAGGGCTAATTCACGGACGCTGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATATCCCGAAGGTGTTCGGTTCAAATCATGCCCCCTCAACCAAATAAAATTTGATGAAGATAAGCAGTACCGTGACGCGGGGTGGagc                                                           <  6:63857/149‑1 (MQ=14)                                                       cgcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTaaaa                                                                                   <  1:392974/110‑2 (MQ=2)                                                       cgcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTaaaa                                                                                   >  2:392974/1‑109 (MQ=2)                                                                                   cTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                             <  2:226049/73‑1 (MQ=255)                                                                                   cTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                             >  1:226049/1‑73 (MQ=255)                                                                                    tGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                             >  6:29763/1‑104 (MQ=255)                                                                                    tGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGa                                                             <  5:29763/104‑1 (MQ=255)                                                                                                                            aaaTCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGc                                                >  6:350833/1‑77 (MQ=255)                                                                                                                                       ccGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                               >  5:353677/1‑51 (MQ=255)                                                                                                                                       ccGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                               <  6:353677/51‑1 (MQ=255)                                                                                                                                               aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccccc  >  7:281440/1‑104 (MQ=40)                                                                                                                                               aaTTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccccc  <  8:281440/104‑1 (MQ=40)                                                                                                                                                   |                                                                                                    ATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCC  >  NC_000913/2947321‑2947566

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