Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,855,981	T→C	G109G (GGA→GGG)	yidH ←	conserved inner membrane protein YidH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,855,981	0	T	C	100.0%	160.2 / NA	45	G109G (GGA→GGG)	yidH	conserved inner membrane protein YidH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/25);  total (20/25)

AGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTT  >  NC_000913/3855906‑3856057                                                                            |                                                                             aGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCAt                                                                         <  1:694749/81‑1 (MQ=255) aGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCAt                                                                         <  2:694747/81‑1 (MQ=255) aGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCAt                                                                         <  1:694747/81‑1 (MQ=255)  gTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCAt                                                                         >  1:694801/1‑80 (MQ=255)       gtgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGAt                                                                   >  1:694815/1‑81 (MQ=255)        tgCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATc                                                                  <  2:694749/81‑1 (MQ=255)           gcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCAcc                                                               >  1:694813/1‑81 (MQ=255)           gcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCAcc                                                               <  1:694751/81‑1 (MQ=255)           gcTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCAcc                                                               <  2:694751/81‑1 (MQ=255)              ccGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGc                                                             >  1:694828/1‑80 (MQ=255)               cGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGcc                                                            <  1:694753/80‑1 (MQ=255)               cGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCa                                                           >  1:694802/1‑81 (MQ=255)               cGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCa                                                           >  1:694826/1‑81 (MQ=255)                 gCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCAc                                                          >  2:694822/1‑80 (MQ=255)                 gCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACg                                                         <  2:694753/81‑1 (MQ=255)                  cTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGa                                                        <  2:694754/81‑1 (MQ=255)                    gCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGAc                                                       >  1:694800/1‑80 (MQ=255)                     cAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGAc                                                       <  2:694755/79‑1 (MQ=255)                      aGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCAt                                                    >  1:694820/1‑81 (MQ=255)                       gTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATc                                                   <  2:694756/81‑1 (MQ=255)                         ccTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCaa                                                 >  1:694814/1‑81 (MQ=255)                          cTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCaa                                                 <  1:694758/80‑1 (MQ=255)                          cTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAat                                                <  2:694758/81‑1 (MQ=255)                          cTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAat                                                <  2:694759/81‑1 (MQ=255)                             ggTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAatat                                              <  1:694761/80‑1 (MQ=255)                             ggTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATAtt                                             <  2:694761/81‑1 (MQ=255)                             ggTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATAtt                                             <  1:694763/81‑1 (MQ=255)                                  gCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCt                                        <  2:694763/81‑1 (MQ=255)                                   cAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTa                                       <  1:694765/81‑1 (MQ=255)                                        cGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGAt                                  >  1:694811/1‑81 (MQ=255)                                         gTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGAt                                  >  1:694812/1‑80 (MQ=255)                                         gTCTGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGAt                                  >  2:694828/1‑80 (MQ=255)                                            tGGCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaa                               >  2:694808/1‑80 (MQ=255)                                             ggCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaat                              >  2:694818/1‑80 (MQ=255)                                             ggCTTTGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATtaat                              >  1:694825/1‑80 (MQ=255)                                                tttGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGc                          >  1:694816/1‑81 (MQ=255)                                                 ttGCGGCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCt                         >  1:694817/1‑81 (MQ=255)                                                      gCTATCCGGCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTgg                    <  2:694765/81‑1 (MQ=255)                                                            cGGCATACAACACCAGCCCCATAACGATCGCCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAg               <  1:694767/80‑1 (MQ=255)                                                             ggCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAggg             >  1:694824/1‑81 (MQ=255)                                                              gCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAggg             >  2:694832/1‑80 (MQ=255)                                                              gCATACAACACCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGc            <  2:694767/81‑1 (MQ=255)                                                                   caacaCCAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTc       <  1:694769/81‑1 (MQ=255)                                                                        ccAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCtt  <  2:694769/81‑1 (MQ=255)                                                                        ccAGCCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCtt  <  2:694770/81‑1 (MQ=255)                                                                            |                                                                             AGTGATGTGCGCTCCGGCTGCAGTCCTGGGTCGGCAATGCGTCTGGCTTTGCGGCTATCCGGCATACAACACCAGTCCCATAACGATCACCGCCACGACCATCAATATTAAGCTAATGATTAATAAGCTGTTGGTGTAGGGCAGGTCCTCTT  >  NC_000913/3855906‑3856057

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: ‑