Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,947,465	Δ1 bp	14.6%	intergenic (+2/‑32)	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,465	0	T	.	14.6%	116.0 / 14.5	41	intergenic (+2/‑32)	metZ/metW	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base T (10/25);  new base . (3/3);  total (13/28)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGT  >  NC_000913/2947317‑2947599                                                                                                                                                     |                                                                                                                                        ccAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAt                                                                                                                                         >  8:40627/1‑149 (MQ=255) ccAAATCCGTGAGCGGGGTAAAACTCCGGGGGCGTGAAAAACCGAGCGGAGTATAGAGCGCTTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAAt                                                                                                                                         <  4:262804/149‑1 (MQ=17)    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                                                                                                                        >  4:186236/1‑141 (MQ=255)                        aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAG‑t                                                                                                                        <  3:186236/141‑1 (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                                               >  1:175492/1‑100 (MQ=255)                                                                                                                                           ccGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCa                                               <  2:175492/100‑1 (MQ=255)                                                                                                                                                   aaTTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt  <  7:160659/137‑1 (MQ=255)                                                                                                                                                   aaTTAAAATTTGATGAAG‑TAAAGCAGTAATGTGACGCGGTGTGGAGCATCCTGGTAGCTCGTAGGGCTCATAACCCTAAGATATTCGGTTCAAATCATGAACCAGAAACAAATAAAATTTGATGAAGTAAAAGCatt  >  8:160659/1‑135 (MQ=255)                                                                                                                                                     |                                                                                                                                        CCAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAG‑TAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGT  >  NC_000913/2947317‑2947599

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