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%	146.4 / NA	35	intergenic (+19/‑96)	thrT/tufB	tRNA‑Thr/translation elongation factor Tu 2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/19);  total (16/19)

TGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTT  >  NC_000913/4175770‑4175919                                                                               |                                                                        tGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGttt                                                                       >  2:849393/1‑81 (MQ=255)    tAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttttc                                                                    <  2:849337/81‑2 (MQ=255)       agCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGtttttcctt                                                                 <  2:849338/81‑5 (MQ=255)        gcgcACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTc                                                               >  2:849383/1‑81 (MQ=255)         cgcACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCt                                                              >  2:849387/1‑81 (MQ=255)          gcACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTg                                                             >  2:849375/1‑81 (MQ=255)           cACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGt                                                            <  2:849339/81‑1 (MQ=255)            aCCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGtt                                                           <  1:849340/81‑1 (MQ=255)               cTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTAt                                                        <  2:849341/81‑1 (MQ=255)                ttGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTAtt                                                       >  2:849380/1‑81 (MQ=255)                  ggTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTAtt                                                       <  2:849342/79‑1 (MQ=255)                   gTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCa                                                    >  2:849388/1‑81 (MQ=255)                    tAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCAt                                                   <  1:849343/81‑1 (MQ=255)                     aaGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCAtt                                                  >  1:849379/1‑81 (MQ=255)                      aGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTc                                                 <  2:849344/81‑1 (MQ=255)                          tGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTcaac                                              >  1:849395/1‑80 (MQ=255)                                cGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCgg                                       >  2:849390/1‑81 (MQ=255)                                 ggCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCggg                                      >  2:849405/1‑81 (MQ=255)                                 ggCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCggg                                      <  2:849345/81‑1 (MQ=255)                                  gCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGc                                     <  1:849346/81‑1 (MQ=255)                                     gTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATg                                  >  2:849399/1‑81 (MQ=255)                                       tCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGtt                                >  2:849411/1‑81 (MQ=255)                                        cGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGtt                                <  1:849347/80‑1 (MQ=255)                                         gAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGc                              <  2:849349/81‑1 (MQ=255)                                         gAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGc                              <  1:849348/81‑1 (MQ=255)                                           aTCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCt                            >  1:849385/1‑81 (MQ=255)                                           aTCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCt                            <  1:849350/81‑1 (MQ=255)                                                ccTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGt                         <  2:849351/79‑1 (MQ=255)                                                ccTATCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTg                       <  2:849352/81‑1 (MQ=255)                                                    tCAGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGAtgt                   <  2:849353/81‑1 (MQ=255)                                                      aGCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCTACAAGTCGGGCATGTTGCCTGGTTGATgtgg                 >  2:849408/1‑81 (MQ=255)                                                       gCACCACTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATgtggt                <  2:849354/81‑1 (MQ=255)                                                         accacTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGa              <  1:849355/81‑1 (MQ=255)                                                          ccacTTCTTTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGa              >  1:849397/1‑80 (MQ=255)                                                             cttcttTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGatat           >  1:849396/1‑80 (MQ=255)                                                                cttTTCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCAcc       >  2:849410/1‑81 (MQ=255)                                                                   ttCTCCTCCCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCAcc       <  1:849356/78‑1 (MQ=255)                                                                        ctccCTG‑TTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGAttt  <  1:849357/78‑1 (MQ=255)                                                                               |                                                                        TGGTAGAGCGCACCCTTGGTAAGGGTGAGGTCGGCAGTTCGAATCTGCCTATCAGCACCACTTCTTTTCTCCTCCCTGTTTTTTCCTTCTGTTTATTGCATTCAACAAGTCGGGCATGTTGCCTGGTTGATGTGGTGATATCACCGATTT  >  NC_000913/4175770‑4175919

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