Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,191,432	C→T	G72S (GGC→AGC)	yehC ←	putative periplasmic pilin chaperone

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,191,432	0	C	T	100.0%	211.0 / NA	59	G72S (GGC→AGC)	yehC	putative periplasmic pilin chaperone
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/27);  total (32/27)

ATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCG  >  NC_000913/2191335‑2191525                                                                                                  |                                                                                              aTTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtat                                                                                             <  2:1261664/100‑1 (MQ=255)  tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTTGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                            >  2:1491599/1‑100 (MQ=255)  tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                            <  1:581910/100‑1 (MQ=255)  tttttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatc                                                                                            >  2:1491940/1‑100 (MQ=255)    tttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatcat                                                                                          <  2:1050529/100‑1 (MQ=255)    tttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCtatcat                                                                                          <  2:308336/100‑1 (MQ=255)     ttACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTatcatc                                                                                         <  1:18192/100‑1 (MQ=255)         ttGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATc                                                                                     <  1:1330519/100‑1 (MQ=255)          tgctgcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATcc                                                                                    <  2:964995/100‑1 (MQ=255)              gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                >  1:6954/1‑100 (MQ=255)              gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                >  1:196167/1‑100 (MQ=255)              gcCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGc                                                                                >  2:544679/1‑100 (MQ=255)                  ggAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgc                                                                            <  1:54940/100‑1 (MQ=255)                     aTTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                         <  2:678014/100‑1 (MQ=255)                      tttGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTgcagc                                                                         <  1:1075145/99‑1 (MQ=255)                         gCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCg                                                                     >  1:323489/1‑100 (MQ=255)                          ccccTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGATTTTTTTCTGGTGGTAATGACGTAGCGCTATCATCATTCCACGCCTGCAGCAGCga                                                                    <  2:656831/100‑1 (MQ=255)                             cTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCgaaga                                                                 <  2:1016328/100‑1 (MQ=255)                             cTATTTTTGCCACTGGTGGAGTTAACTTGTAAGGAACCTGAATTTTTTCTGGTGCTAATGACGTATCGCTATCATCTATCCACGCCTGCagcagc                                                                      <  2:1342008/95‑1 (MQ=255)                                 ttttGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACgg                                                             >  1:1182641/1‑100 (MQ=255)                                     gCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTc                                                         >  1:184686/1‑100 (MQ=255)                                      ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                        >  1:893222/1‑100 (MQ=255)                                      ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTcc                                                        >  1:1535917/1‑100 (MQ=255)                                      ccACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTATcc                                                        >  1:40508/1‑100 (MQ=255)                                       cACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTccc                                                       >  2:1244489/1‑100 (MQ=255)                                         cTGGTGGCGTTAACATGAAAGGAACCTGAATTTCTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTg                                                     <  2:1405724/100‑1 (MQ=255)                                          tggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTgg                                                    <  1:647269/100‑1 (MQ=255)                                          tggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTg                                                     <  2:313776/99‑1 (MQ=255)                                           ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                   <  1:422273/100‑1 (MQ=255)                                           ggtggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                   >  1:919965/1‑100 (MQ=255)                                             tggCGTTAACATGAAATGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGt                                                   <  2:196166/98‑1 (MQ=255)                                             tggCGTTAACATGAAAGGAACCTGAATTTTTTTTTTTTTTATTGACGTATAGCTATCATCAATTTATTTTTGCAGCAGCGAAGAACGGttttctttttt                                                  >  2:667018/1‑91 (MQ=255)                                              ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtca                                                >  1:201800/1‑100 (MQ=255)                                              ggCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtca                                                >  1:1151666/1‑100 (MQ=255)                                                cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                              <  1:490325/100‑1 (MQ=255)                                                cGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTtcatc                                              >  1:235253/1‑100 (MQ=255)                                                   tAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAAc                                           >  1:306893/1‑100 (MQ=255)                                                    aaCATGAAAGGAATCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACt                                          >  2:974413/1‑100 (MQ=255)                                                      cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                        >  2:670334/1‑100 (MQ=255)                                                      cATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGc                                        >  1:383136/1‑100 (MQ=255)                                                          aaaGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCa                                    <  2:6954/100‑1 (MQ=255)                                                            aGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATc                                  <  2:323487/100‑1 (MQ=255)                                                                      aTTTTTTCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTAt                        <  1:646942/100‑1 (MQ=255)                                                                          tttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                    <  2:1448386/100‑1 (MQ=255)                                                                          tttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                    <  2:742669/100‑1 (MQ=255)                                                                          tttCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTc                    <  1:358209/100‑1 (MQ=255)                                                                            tCTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTg                  <  2:243485/100‑1 (MQ=255)                                                                             cTGGTGGTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGc                 >  1:791387/1‑100 (MQ=255)                                                                              tggtggTAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGcc                >  1:1484841/1‑100 (MQ=255)                                                                               ggtggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCg               >  2:1220318/1‑100 (MQ=255)                                                                                 tggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCggg             >  1:320159/1‑100 (MQ=255)                                                                                 tggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCggg             >  1:128829/1‑100 (MQ=255)                                                                                 tggtAATGACGTATCGCTATCATCAACCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCggg             >  2:399003/1‑100 (MQ=255)                                                                                  ggtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGt            <  2:673040/100‑1 (MQ=255)                                                                                   gtAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTa           >  1:1288120/1‑100 (MQ=255)                                                                                    tAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTaa          <  2:145928/100‑1 (MQ=255)                                                                                    tAATGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTaa          >  1:1270149/1‑100 (MQ=255)                                                                                       tGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaata       >  2:309833/1‑100 (MQ=255)                                                                                       tGACGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTGCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaata       <  2:1078792/100‑1 (MQ=255)                                                                                         aCGTATCGCTATCATCAATCCACGCCTGCAGCAGCGACGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGGAaataat     >  2:844921/1‑100 (MQ=255)                                                                                         aCGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAaataat     >  2:482162/1‑100 (MQ=255)                                                                                          cGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATg    <  1:309835/100‑1 (MQ=255)                                                                                          cGTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATg    >  1:360506/1‑100 (MQ=255)                                                                                           gTATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgc   <  1:1114147/100‑1 (MQ=255)                                                                                            tATCGCTATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATgcg  >  1:658731/1‑100 (MQ=255)                                                                                                  |                                                                                              ATTTTTACTTGCTGCCCGGAATTTGCCCCTATTTTTGCCACTGGTGGCGTTAACATGAAAGGAACCTGAATTTTTTCTGGTGGTAATGACGTATCGCCATCATCAATCCACGCCTGCAGCAGCGAAGAACGGTTTCCCTGGTTCATCAACTGCACCATCACTTCTTTATTTTCTGCCGGGTAAATAATGCG  >  NC_000913/2191335‑2191525

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