Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	638,399	A→G	Y59H (TAC→CAC)	dsbG ←	thiol:disulfide interchange protein, periplasmic

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	638,399	0	A	G	100.0%	133.3 / NA	36	Y59H (TAC→CAC)	dsbG	thiol:disulfide interchange protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/18);  total (18/18)

GTTCCATCCGTTGCCACATTTCGCGTCCGGCTGGTGCGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTAGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGG  >  NC_000913/638267‑638513                                                                                                                                     |                                                                                                                   gTTCCATCCGTTGCCACATTTCGCGTCCGGCTGGTGCGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGt                                                                                                              <  1:671193/139‑1 (MQ=255)  ttCCATCCGTTGCCACATTTCGCGTCCGGCTGGTGCGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTg                                                                                                             >  1:225072/1‑139 (MQ=255)            tGCCACATTTCGCGTCCGGCTGGTGCGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCatat                                                                                                   >  1:655344/1‑139 (MQ=255)                 caTTTCGCGTCCGGCTGGTGCGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGa                                                                                              <  2:114233/139‑1 (MQ=255)                        cgTCCGGCTGGTGCGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGa                                                                                                                  >  1:744254/1‑112 (MQ=255)                        cgTCCGGCTGGTGCGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTc                                                                                       <  2:477840/139‑1 (MQ=255)                              gCTGGTGCGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGAt                                                                                 <  2:779526/139‑1 (MQ=255)                                     cGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCttt                                                                          >  2:281871/1‑139 (MQ=255)                                                         aGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGa                                                      <  1:281871/139‑1 (MQ=255)                                                             tgtTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACg                                                                                                          >  1:15820/1‑83 (MQ=255)                                                             tgtTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACg                                                                                                          <  2:15820/83‑1 (MQ=255)                                                                tACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGtttt                                               <  2:648002/139‑1 (MQ=255)                                                                     aGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTtgat                                            <  1:674736/137‑1 (MQ=255)                                                                     aGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTtgat                                            >  2:674736/1‑137 (MQ=255)                                                                          ttCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAAt                                                    <  2:361554/124‑1 (MQ=255)                                                                          ttCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAAt                                                    >  1:361554/1‑124 (MQ=255)                                                                               cTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGcc                                <  2:655344/139‑1 (MQ=255)                                                                                tttCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGccc                               <  2:225072/139‑1 (MQ=255)                                                                                            cATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATc                   >  1:340404/1‑139 (MQ=255)                                                                                                  aCCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGtt                                                 >  1:700479/1‑103 (MQ=255)                                                                                                  aCCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGtt                                                 <  2:700479/103‑1 (MQ=255)                                                                                                  aCCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCATCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATCGCtttt             >  1:124938/1‑139 (MQ=255)                                                                                                       agaTAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGa                                                                                                            >  1:770034/1‑39 (MQ=255)                                                                                                       agaTAGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGa                                                                                                            <  2:770034/39‑1 (MQ=255)                                                                                                           aGCGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCa    >  2:304432/1‑139 (MQ=255)                                                                                                             cGTGCTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAgg  >  2:582410/1‑139 (MQ=255)                                                                                                                 cTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCgaga                                                                                  <  1:777974/55‑1 (MQ=255)                                                                                                                 cTTACCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCgaga                                                                                  >  2:777974/1‑55 (MQ=255)                                                                                                                    aCCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATCGCTTTTACTGGa       >  2:257176/1‑127 (MQ=255)                                                                                                                    aCCATCTGGAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATCGCTTTTACTGGa       <  1:257176/127‑1 (MQ=255)                                                                                                                            gAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCAt                                                                       <  2:686501/55‑1 (MQ=255)                                                                                                                            gAGTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCAt                                                                       >  1:686501/1‑55 (MQ=255)                                                                                                                              gTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCgaga                                                                                  >  1:87329/1‑42 (MQ=255)                                                                                                                              gTCAGGTGGATGGTGACGCCCATATCCTGATACTTTCCgaga                                                                                  <  2:87329/42‑1 (MQ=255)                                                                                                                                cAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAg      <  1:432414/116‑1 (MQ=255)                                                                                                                                cAGGTGGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAg      >  2:432414/1‑116 (MQ=255)                                                                                                                                     |                                                                                                                   GTTCCATCCGTTGCCACATTTCGCGTCCGGCTGGTGCGTAAATTTCTTTTTCGATAAGTGTGTTACTCAGGTTTTCACCTTTCTCGTTGTACATGTAACCAGAGATAGCGTGCTTACCATCTGGAGTCAGGTAGATGGTGACGCCCATATCCTGATACTTTCCGAGATAACCTTTCATTCCTCCGGGGGCATCGAATGTTTTGATGATTGTAATGCCCTGTTTTTCAATCGCTTTTACTGGAGCAGG  >  NC_000913/638267‑638513

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