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%	35.7 / NA	13	F172F (TTC→TTT)	queG	epoxyqueuosine reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/8);  total (5/8)

GATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTGCTG  >  NC_000913/4393312‑4393798                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       gATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCACGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCCTCCAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAaaaaga                                                                                                                                                                                                                                                    >  1:246240/1‑245 (MQ=255)         ggCAGTCATCGCAACAGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTcg                                                                                                                                                                                                                                      >  1:364146/1‑251 (MQ=255)                                                                       ccTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAACTCGTCATCCAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGACTATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTAcc                                                                                                                                                                                                            >  1:251489/1‑215 (MQ=255)                                                                             aGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCGCGACGATGGCACCGGGCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCTGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCgag                                                                                                                                                                  <  1:167449/251‑1 (MQ=255)                                                                                   aTGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAgn                                                                                                                                                            <  1:270867/251‑2 (MQ=255)                                                                                                                                                   ggCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAg                                                                                                                                                             <  1:376164/186‑1 (MQ=255)                                                                                                                                                                                       cttcctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCAtctcn                                                                                                  <  1:130374/209‑2 (MQ=255)                                                                                                                                                                                        ntcctCGNNTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACTAACCGGCCTCGCGATTGAGGATAagtgagttc                                                                                                                                                                                                                 <  2:1818/96‑3 (MQ=255)                                                                                                                                                                                             cGACTGGTTGATCCACGGGCAGCGGAATAGCGACCAGCAATGCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCATCTTGTTCAGCTAACGGGCGCTCGAGAAGAGGCGCAGAATCTACAAACGGTCTAAAATTCATCGAAACACAATGTTGCTGAATCAGCGCGCCCAGCTGTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGcc                                                   <  1:293641/250‑1 (MQ=255)                                                                                                                                                                                                        tattccgctgcccgtggatcaaacAGCAATTCGCCTAAAAAAAAGAACGACCCGGCCTCGCGATTGAGGATAAGAGAGAGCTTACATGTCACGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCTCCCAGATTTTTGAGTCGGTTGCGCAGAAGATTGTGATAGTCACGGACCAGCGCATAACg                                       >  1:182134/24‑251 (MQ=255)                                                                                                                                                                                                                    ggAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATACCGGCTAACATAGcc                           >  1:18640/1‑251 (MQ=255)                                                                                                                                                                                                                             nccAGCANNTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTCCAGCTAACg                                                                                                                                                                            <  2:5439/96‑1 (MQ=255)                                                                                                                                                                                                                                 gCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTGCCTGTCCAGCCGAGCCCAGCTTTTTTAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTg                                                                             >  1:230867/1‑188 (MQ=255)                                                                                                                                                                                                                                          cTAAAAAAAACAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGc                                                                                                                                                                                          <  1:398244/70‑1 (MQ=255)                                                                                                                                                                                                                                            aaaaaaaaGAACGCACTGGCCTCGCGNTTGCGTNNAANTGAGTGCTTACCTGTCCAGNCGNNNCNNNCNTTTTCNGCNNACGGGAGCTNGAGAATAGGCGCCGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTNCTGACTCCTNNCGCCCAGCTTTTTGAGACGGTNGCGNNGNNGNTtgng                                                               >  2:96226/1‑191 (MQ=255)                                                                                                                                                                                                                                             aaaaaaAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTGCTg  <  1:416793/251‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                       GATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTGCTG  >  NC_000913/4393312‑4393798

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