Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_012967	3,733,506	A→C	intergenic (‑189/‑25)	kbl ← / → rfaD	glycine C‑acetyltransferase/ADP‑glyceromanno‑heptose 6‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_012967	3,733,506	0	A	C	100.0%	273.6 / NA	66	intergenic (‑189/‑25)	kbl/rfaD	glycine C‑acetyltransferase/ADP‑glyceromanno‑heptose 6‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (35/31);  total (35/31)

CTCACTTTTGGCGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGG  >  NC_012967/3733430‑3733582                                                                             |                                                                             ctcACTTTTGGCGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCa                                                                         <  1:1220702/82‑1 (MQ=255)  tcACTTTTGGCGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCAt                                                                        >  1:970194/1‑82 (MQ=255)    aCTTTTGGCGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATac                                                                      >  1:2657531/1‑81 (MQ=255)     cTTTTGGCGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAAt                                                                     <  1:4441586/82‑1 (MQ=255)      ttttGGCGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAAtt                                                                    >  1:1913886/1‑82 (MQ=255)          ggCGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGaa                                                                >  1:532573/1‑82 (MQ=255)          ggCGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGaa                                                                <  1:1847719/82‑1 (MQ=255)            cGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAgg                                                              >  1:2595731/1‑82 (MQ=255)            cGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAga                                                              >  1:1055967/1‑81 (MQ=255)             gAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGt                                                             >  1:1186811/1‑82 (MQ=255)             gAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGt                                                             <  1:569553/82‑1 (MQ=255)              aaTGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGtt                                                            >  1:3376255/1‑82 (MQ=255)                tGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTAc                                                          <  1:2425000/82‑1 (MQ=255)                tGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTAc                                                          <  1:2425001/82‑1 (MQ=255)                     gAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTa                                                     >  1:4171885/1‑82 (MQ=255)                     gAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTa                                                     <  1:3135157/82‑1 (MQ=255)                      aaGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTAt                                                    >  1:153952/1‑82 (MQ=255)                       aGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATg                                                   >  1:1547602/1‑82 (MQ=255)                       aGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATg                                                   <  1:2214829/82‑1 (MQ=255)                        gCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGa                                                  <  1:1318309/82‑1 (MQ=255)                         ccAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGat                                                 <  1:26131/82‑1 (MQ=255)                          cAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGatc                                                <  1:522461/82‑1 (MQ=255)                             cGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGatcatc                                             >  1:2000187/1‑82 (MQ=255)                                   gtAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTAcc                                       >  1:3807506/1‑82 (MQ=255)                                   gtAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTAcc                                       >  1:3385785/1‑82 (MQ=255)                                    tAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcg                                      >  1:2905637/1‑82 (MQ=255)                                      gTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcggc                                    >  1:1414957/1‑82 (MQ=255)                                      gTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcggc                                    >  1:3437874/1‑82 (MQ=255)                                        gatatgataGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcggcgg                                  <  1:2052511/82‑1 (MQ=255)                                        gatatgataGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACcggcgg                                  >  1:1136768/1‑82 (MQ=255)                                         atatgataGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCggcggc                                 >  1:1101512/1‑82 (MQ=255)                                         atatgataGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCggcggc                                 <  1:1267117/82‑1 (MQ=255)                                         atatgataGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCggcggc                                 <  1:3442910/82‑1 (MQ=255)                                          tatgataGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGgcggcg                                >  1:1108157/1‑82 (MQ=255)                                            tgataGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGcgcg                              <  1:3801895/82‑1 (MQ=255)                                             gataGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCgg                             <  1:1699387/82‑1 (MQ=255)                                              ataGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCggg                            >  1:2973118/1‑82 (MQ=255)                                              ataGAAGTATTCGTGTCTGAGATTGCCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCggg                            <  1:720169/82‑1 (MQ=255)                                                aGAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCt                          >  1:2185057/1‑82 (MQ=255)                                                 gAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCtt                         >  1:2653636/1‑82 (MQ=255)                                                 gAAGTATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCtt                         >  1:244088/1‑82 (MQ=255)                                                     tATTCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATc                     <  1:4432138/82‑1 (MQ=255)                                                       ttCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTCTGATCATCGTTACCAGAGGCGCGGGCTTTATCgg                   >  1:3506894/1‑82 (MQ=255)                                                       ttCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCgg                   >  1:3506895/1‑82 (MQ=255)                                                        tCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGgc                  >  1:2999501/1‑82 (MQ=255)                                                        tCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGgc                  <  1:491486/82‑1 (MQ=255)                                                        tCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGgc                  <  1:1266687/82‑1 (MQ=255)                                                        tCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGgc                  <  1:1974194/82‑1 (MQ=255)                                                        tCGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGgc                  <  1:159667/82‑1 (MQ=255)                                                         cGTGTCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGgca                 <  1:4381058/82‑1 (MQ=255)                                                             tCTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCaa             >  1:1908472/1‑82 (MQ=255)                                                             tCTGAGATTGTCTCTGCCCCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCaa             >  1:1074682/1‑82 (MQ=255)                                                              ttGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAAc            <  1:4102713/81‑1 (MQ=255)                                                              cTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAAc            <  1:3845916/82‑1 (MQ=255)                                                              cTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAAc            <  1:3287858/82‑1 (MQ=255)                                                              cTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAAc            >  1:2838994/1‑82 (MQ=255)                                                              cTGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAAc            <  1:1502093/82‑1 (MQ=255)                                                               tGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACa           >  1:1087807/1‑82 (MQ=255)                                                               tGAGATTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTATATCGGCAGCAACa           >  1:410887/1‑82 (MQ=255)                                                                gagaTTGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACAt          <  1:553907/82‑1 (MQ=255)                                                                   gttGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGt       <  1:119385/81‑1 (MQ=255)                                                                     tGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTa     >  1:3859444/1‑82 (MQ=255)                                                                     tGTCTCTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTa     >  1:3742811/1‑82 (MQ=255)                                                                       tctcTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAg   <  1:3029118/82‑1 (MQ=255)                                                                       tctcTGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAg   >  1:685015/1‑82 (MQ=255)                                                                        ctctGCCTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAgg  <  1:4202553/82‑1 (MQ=255)                                                                             |                                                                             CTCACTTTTGGCGAATGGCTGAAGCCAACGCCGTGTAGTGATATGATAGAAGTATTCGTGTCTGAGATTGTCTCTGACTCCATAATTCGAAGGTTACAGTTATGATCATCGTTACCGGCGGCGCGGGCTTTATCGGCAGCAACATCGTTAAGG  >  NC_012967/3733430‑3733582

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