Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,947,467	A→C	12.4%	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.4%	92.2 / 4.8	40	intergenic (+4/‑30)	metZ/metW	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (9/25);  new base C (2/3);  total (12/28)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGG  >  NC_000913/2947320‑2947612                                                                                                                                                    |                                                                                                                                                   aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTa                                                                                                                                                    >  7:372897/1‑148 (MQ=255)   tCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaaa                                                                                                                                                 <  7:437545/149‑1 (MQ=255)    ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                                <  5:177858/149‑1 (MQ=255)     cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAt                                                                                                                                                <  2:197181/148‑1 (MQ=255)      gTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAAtt                                                                                                                                               >  3:134792/1‑148 (MQ=255)         aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                           <  8:52165/149‑1 (MQ=255)         aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                                                                            <  3:342160/148‑1 (MQ=255)         aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtga                                                                                                                                            <  6:215633/148‑1 (MQ=255)          tCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCTCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgat                                                                                                                                           <  1:159308/148‑1 (MQ=255)            ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                                        <  4:182429/149‑1 (MQ=255)             gggTAAAAATAAGGTTGCATGAAAACGATAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑t                                                                                                                                     <  4:134792/149‑1 (MQ=18)              ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑t                                                                                                                                     <  1:168478/149‑1 (MQ=255)                tAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑Taa                                                                                                                                   <  1:184663/149‑1 (MQ=255)                    aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCa                                                                                                                               <  7:40627/149‑1 (MQ=255)                    aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCa                                                                                                                               <  4:400165/149‑1 (MQ=255)                     aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑t                                                                                                                                     <  3:186236/141‑1 (MQ=255)                     aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑t                                                                                                                                     >  4:186236/1‑141 (MQ=255)                       aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTa                                                                                                                            <  8:372897/149‑1 (MQ=255)                               aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGAc                                                                                                                     >  4:202708/1‑148 (MQ=255)                               aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGAc                                                                                                                     <  3:202708/148‑1 (MQ=255)                                  aaaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCgg                                                                                                                 <  2:101658/149‑1 (MQ=255)                                   aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCggg                                                                                                                <  1:256963/149‑1 (MQ=255)                                    aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAACTCCGTCCCACGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCgggg                                                                                                               <  4:13028/149‑1 (MQ=255)                                     aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGG‑CCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTg                                                                                                             <  1:6601/149‑1 (MQ=255)                                      cgcgAGCGGAGTATAGTGCCCATCCCCGGACGCTGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGCCCCCCGCAACCAATTAAAATTTGATGAAG‑TAACGCAGTACTGGGACGCGGGGTg                                                                                                             <  2:219703/149‑1 (MQ=18)                                           gCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTT‑AAGAAGTTAAAGCAGTACGGTGAAGCggggagatgg                                                                                                         >  6:174367/1‑142 (MQ=17)                                                gTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCGCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACTGTGACGCGGGGTGGAGCAGCCTg                                                                                                   <  8:63333/149‑1 (MQ=255)                                                      tGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgat                                                                                                                                           >  5:20044/1‑103 (MQ=255)                                                      tGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgat                                                                                                                                           <  6:20044/103‑1 (MQ=255)                                                           aTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TGAAATTtga                                                                                                                                            >  8:389695/1‑97 (MQ=255)                                                            tCCACGGACGCGGGGGGGAGCAGCCTGGTAGCTCTTCGGGCTCATATCCCGAAGGTCGTCGGTTCAAATCCGGCCCCGGCATGCAATTAAAATTTGATGAAG‑TAAAGCAGTACTGTGACGCGGGGTGGa                                                                                                           <  7:416742/129‑1 (MQ=17)                                                                          gTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTAATTATT‑AAAAGCAGTACGGTGAc                                                                                                                     >  4:269523/1‑105 (MQ=11)                                                                                   ccTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTg                                                                                                             <  7:42109/104‑1 (MQ=255)                                                                                   ccTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTg                                                                                                             >  8:42109/1‑104 (MQ=255)                                                                                     tGGTATCTCGTCGGGCTCATAACCCGAAGGTCGTAGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACg                                                                                                                          >  6:17340/1‑89 (MQ=255)                                                                                     tGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACg                                                                                                                          <  5:17340/89‑1 (MQ=255)                                                                                                            ccGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCggg                                                                                      >  8:68258/1‑102 (MQ=255)                                                                                                                                        ccGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCa                                                            <  2:175492/100‑1 (MQ=255)                                                                                                                                        ccGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCa                                                            >  1:175492/1‑100 (MQ=255)                                                                                                                                                aaTTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt               <  7:160659/137‑1 (MQ=255)                                                                                                                                                aaTTAAAATTTGATGAAG‑TAAAGCAGTAATGTGACGCGGTGTGGAGCATCCTGGTAGCTCGTAGGGCTCATAACCCTAAGATATTCGGTTCAAATCATGAACCAGAAACAAATAAAATTTGATGAAGTAAAAGCatt               >  8:160659/1‑135 (MQ=255)                                                                                                                                                   tAAAATTTGCAGCAG‑GAAAGCACTCCGCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCggg  <  7:77157/147‑1 (MQ=14)                                                                                                                                                    |                                                                                                                                                   AATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGG  >  NC_000913/2947320‑2947612

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