Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,175,853	(T)6→5	intergenic (+24/‑91)	thrT → / → tufB	tRNA‑Thr/translation elongation factor Tu 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,175,848	0	T	.	100.0%	134.7 / NA	32	intergenic (+19/‑96)	thrT/tufB	tRNA‑Thr/translation elongation factor Tu 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (25/7);  total (25/7)

GTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATC  >  NC_000913/4175768‑4175922                                                                                 |                                                                           gTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGt                                                                            <  1:781543/81‑1 (MQ=255)    ggTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttt                                                                         >  2:781579/1‑81 (MQ=255)     gTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGttttt                                                                        <  1:781544/81‑1 (MQ=255)         agCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttttcctt                                                                    >  2:781575/1‑77 (MQ=255)            gcACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTg                                                                <  2:781544/81‑1 (MQ=255)                 cTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTAt                                                           >  2:781581/1‑81 (MQ=255)                    ggTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTg                                                         >  1:781585/1‑80 (MQ=255)                        aGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCAtt                                                     >  2:781574/1‑80 (MQ=255)                        aGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTc                                                    <  2:781545/81‑1 (MQ=255)                             gAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTcaac                                                 >  1:781600/1‑79 (MQ=255)                             gAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTcaac                                                 >  1:781589/1‑79 (MQ=255)                              aGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAg                                              >  2:781576/1‑81 (MQ=255)                               ggTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGt                                             <  2:784698/81‑1 (MQ=255)                                gTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTc                                            >  2:781602/1‑81 (MQ=255)                                gTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTc                                            >  1:781605/1‑81 (MQ=255)                                  cGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCgg                                          >  1:781567/1‑81 (MQ=255)                                         tCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGtt                                   >  2:781606/1‑81 (MQ=255)                                          cGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGtt                                   <  2:781546/80‑1 (MQ=255)                                          cGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTg                                  >  2:781577/1‑81 (MQ=255)                                           gAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGc                                 >  2:781587/1‑81 (MQ=255)                                           gAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGc                                 >  2:781586/1‑81 (MQ=255)                                           gAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGc                                 >  2:781582/1‑81 (MQ=255)                                            aaTCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGcc                                >  1:781590/1‑81 (MQ=255)                                               cTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTgg                             >  1:781609/1‑81 (MQ=255)                                                 gCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGtt                           >  2:781580/1‑81 (MQ=255)                                                    tATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAt                        >  2:781590/1‑81 (MQ=255)                                                     aTCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAtg                       >  2:781596/1‑81 (MQ=255)                                                      tCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAt                        <  2:781547/79‑1 (MQ=255)                                                           accacTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGa                 >  2:781593/1‑81 (MQ=255)                                                                ttcttTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCa            >  1:781596/1‑81 (MQ=255)                                                                  cttTTCTCCTCCCTG‑TTTTTCCTTCTGTTTGTTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCAc           <  2:781548/80‑1 (MQ=255)                                                                  cttTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCAc           <  2:781549/80‑1 (MQ=255)                                                                    tttCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGa        >  1:781593/1‑81 (MQ=255)                                                                      tctcCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAtt      >  1:781584/1‑81 (MQ=255)                                                                      tctcCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAtt      >  2:781594/1‑81 (MQ=255)                                                                          ctccCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATc  >  2:781607/1‑81 (MQ=255)                                                                                 |                                                                           GTTGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTTATC  >  NC_000913/4175768‑4175922

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