Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,393,552	G→A	F172F (TTC→TTT)	queG ←	epoxyqueuosine reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,393,552	0	G	A	100.0%	109.0 / NA	31	F172F (TTC→TTT)	queG	epoxyqueuosine reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/21);  total (10/21)

AACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGT  >  NC_000913/4393332‑4393794                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   aaCCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATa                                                                                                                                                                                                                      <  1:93562/251‑1 (MQ=255) aaCCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATa                                                                                                                                                                                                                      <  1:76128/251‑1 (MQ=255)      tAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAagtg                                                                                                                                                                                                                  >  1:245064/1‑250 (MQ=255)        aaTACGGTTTCCCATTAACGGTCGCACCTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATACACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAagtgatt                                                                                                                                                                                                               <  2:10142/251‑3 (MQ=255)                  cccATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGGGAGATAAGAGATACAGCGGCGAGCAGCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGt                                                                                                                                                                                                     <  2:201163/251‑1 (MQ=255)                              cGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGAAGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGCGCTTACCTGTCCAGCCGAGccc                                                                                                                                                                                         >  2:30856/1‑251 (MQ=255)                                                               aTGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAgg                                                                                                                                                        >  1:116600/1‑251 (MQ=255)                                                                           gagaTACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGAc                                                                                                                                            <  2:192866/251‑1 (MQ=255)                                                                                                 cGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGa                                                                                                                       >  2:127679/1‑250 (MQ=255)                                                                                                        tGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGcc                                                                                                                                                                                               <  2:243812/171‑1 (MQ=255)                                                                                                        tGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGcc                                                                                                                                                                                               >  1:243812/1‑171 (MQ=255)                                                                                                             cgacgaTGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGc                                                                                                          <  1:60846/251‑1 (MQ=255)                                                                                                             cgacgaTGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGACCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTg                                                                                                           <  2:176102/250‑1 (MQ=255)                                                                                                                cgaTGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGa                                                                                                       <  1:161620/251‑1 (MQ=255)                                                                                                                       aCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTAAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCAtct                                                                                                <  2:176491/251‑1 (MQ=255)                                                                                                                            tCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGccc                                                                                           >  1:75909/1‑251 (MQ=255)                                                                                                                                  aaaTCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCttt                                                                                     <  1:56434/251‑1 (MQ=255)                                                                                                                                  aaaTCGACATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAATAACTAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCAAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAAAAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAAGGTTGCTGAATCATCTCGCCCAGCttt                                                                                     <  2:113777/251‑1 (MQ=255)                                                                                                                                           tACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCgg                                                                            >  2:232324/1‑251 (MQ=255)                                                                                                                                            aCAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGTGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGt                                                                           <  2:126808/251‑1 (MQ=255)                                                                                                                                             cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGtt                                                                          <  1:127679/251‑1 (MQ=255)                                                                                                                                                           gCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTtgt                                                            <  1:43928/251‑1 (MQ=255)                                                                                                                                                             cGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTtgtg                                                           >  1:178630/1‑250 (MQ=255)                                                                                                                                                               cATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATa                                                        <  1:147765/251‑1 (MQ=255)                                                                                                                                                                              ggatGATCCACGGGCATCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGAATGAGGATAATTGAGTGCTTACCTGTCCAGCCGAGCCCAGCGTGTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGTCAAACCGGCTAAAATGCAGCGAAACACAATGTTGCTGAATCATCGCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGACAGTCACGGCGCAgcg                                          <  2:100460/247‑1 (MQ=255)                                                                                                                                                                                 tGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCAt                                       <  2:123641/250‑1 (MQ=255)                                                                                                                                                                                                     atCGACCAGCAATTCGCCTAAAAAAAAGAACGAAACGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATCACGGCTAACATAGCCGCGtt                  >  1:232334/1‑251 (MQ=255)                                                                                                                                                                                                          ccAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTgggg             >  1:222151/1‑251 (MQ=255)                                                                                                                                                                                                             gCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGttt          <  1:113426/251‑1 (MQ=255)                                                                                                                                                                                                                    gCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCGCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCAGCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCg   <  2:17362/251‑1 (MQ=255)                                                                                                                                                                                                                     ccTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGt  <  1:241553/251‑1 (MQ=255)                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   AACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGT  >  NC_000913/4393332‑4393794

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