Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_017635	1,677,854	C→A	intergenic (‑110/+21)	digH ← / ← gadC	glycoside hydrolase family 10 protein/acid resistance gamma‑aminobutyrate antiporter GadC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_017635	1,677,854	0	C	A	100.0%	153.7 / NA	42	intergenic (‑110/+21)	digH/gadC	glycoside hydrolase family 10 protein/acid resistance gamma‑aminobutyrate antiporter GadC
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/22);  total (20/22)

CGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCAC  >  NC_017635/1677808‑1677898                                               |                                             cGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATta                                           <  2:13167/50‑1 (MQ=255)    cAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATTATaa                                        <  2:375285/50‑1 (MQ=255)    cAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATTATaa                                        <  1:1367257/50‑1 (MQ=255)      gcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATTATAAcc                                      <  1:452893/50‑1 (MQ=255)        aagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATTATAACCAt                                    <  2:790554/50‑1 (MQ=255)         agcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATTATAACCAtt                                   <  2:809844/50‑1 (MQ=255)         agcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATTATAACCAtt                                   >  2:1476968/1‑50 (MQ=255)          gcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATTATAACCATTc                                  <  1:274011/50‑1 (MQ=255)             cGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATTATAACCATTCAGa                               >  2:555910/1‑50 (MQ=255)              gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTATTATAACCATTCAGa                               >  2:1828700/1‑49 (MQ=255)               aaaaaGGGAGCGATGAATTATCGCTCCCTTGTATTATAACCATTCAGACa                             >  2:1124016/1‑50 (MQ=255)                    gggAGCGATGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGtt                        <  1:1275713/50‑1 (MQ=255)                    gggAGCGATGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGtt                        <  2:1674871/50‑1 (MQ=255)                     ggAGCGATGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGt                         >  2:596964/1‑48 (MQ=255)                       aGCGATGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAgt                     >  1:1416/1‑50 (MQ=255)                        gCGATGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAgtg                    >  2:1783005/1‑50 (MQ=255)                        gCGATGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAgtg                    <  1:1066771/50‑1 (MQ=255)                         cGATGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAgtgt                   >  2:696964/1‑50 (MQ=255)                          gATGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGtt                  >  1:902374/1‑50 (MQ=255)                          gATGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGtt                  >  1:1593851/1‑50 (MQ=255)                            tGAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTc                >  1:1098237/1‑50 (MQ=255)                             gAATTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCt               >  1:1145299/1‑50 (MQ=255)                              aaTTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCtt              >  1:1261315/1‑50 (MQ=255)                              aaTTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCtt              <  1:1411128/50‑1 (MQ=255)                               aTTATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTg             <  1:2056020/50‑1 (MQ=255)                                ttATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGt            >  1:656825/1‑50 (MQ=255)                                ttATCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGt            <  1:1777968/50‑1 (MQ=255)                                  aTCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtca          >  1:1426083/1‑50 (MQ=255)                                  aTCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtca          >  1:140777/1‑50 (MQ=255)                                  aTCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtc           <  2:2037367/49‑1 (MQ=255)                                   tCGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtca          <  2:1790229/49‑1 (MQ=255)                                    cGCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtcatt        >  1:1190414/1‑50 (MQ=255)                                    cGCTCCCTTGTATTATAACCATTCAGACATGGTAAGTGTTTCTTTccatt        >  2:561680/1‑50 (MQ=255)                                     gCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattc       <  1:2009909/50‑1 (MQ=255)                                     gCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtcatt        <  2:1595387/49‑1 (MQ=255)                                     gCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtcatt        <  2:438394/49‑1 (MQ=255)                                     gCTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtcatt        >  1:1184236/1‑49 (MQ=255)                                      cTCCCTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattca      >  2:782146/1‑50 (MQ=255)                                         ccTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattcat     <  1:984874/48‑1 (MQ=255)                                         ccTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatc    <  2:1067299/49‑1 (MQ=255)                                          cTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATcac  <  1:309742/50‑1 (MQ=255)                                          cTTGTATTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATcac  <  1:609623/50‑1 (MQ=255)                                               |                                             CGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCAC  >  NC_017635/1677808‑1677898

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