Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	1,117,181	0	G→T	30.8%	12.4	40	F195L (TTC→TTA)	csgD	DNA‑binding transcriptional activator in two‑component regulatory system
Reads supporting (aligned to +/- strand):  new base (12/0):  ref base (11/16):  total (23/17)

ACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTC  >  REL606/1117146‑1117216                                    |                                    aCCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTg                                     <  1:4917637/36‑1 (MQ=255) aCCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTg                                     <  1:6893406/36‑1 (MQ=255)    gCTTGTGTCCGGTTTTTTACGGCTATCTTCTTTaaa                                  >  1:6195221/1‑36 (MQ=255)      ttGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAg                                <  1:4653585/36‑1 (MQ=255)      ttGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAg                                <  1:5246257/36‑1 (MQ=255)        gtgtCCGGTTTTTTACGGCTATATTTTTTAAAAAAt                              >  1:1192830/1‑36 (MQ=18)         tgtCCGGTTTTTTACGGCTATCTTCTTGAAAACAtt                             >  1:3525681/1‑36 (MQ=255)          gtCCGGTTTTTTACGGCTATCTTCTTTAAAAAATta                            >  1:834312/1‑36 (MQ=25)          gtCCGGTTTTTTACGGCTATCTTCTTGAAAAGATta                            <  1:6457589/36‑1 (MQ=255)          gtCCGGTTTTTTACGGCTAACTTCTTTAAAAAAtaa                            >  1:2633951/1‑34 (MQ=18)           tCCGGTTTTTTACGGCTATCTTCTTGAAAAGATtat                           <  1:6645093/36‑1 (MQ=255)           tCCGGTTTTTTACGGCTAAATTTTTTTAAAGAtttt                           >  1:3145033/1‑34 (MQ=2)            ccGGTTTTTTACGGCTATCTTCTTGAAAAGATtata                          <  1:5011064/36‑1 (MQ=255)            ccGGTTTTTTACGGCTATATTCTTTTAAAGATtata                          >  1:2982809/1‑36 (MQ=21)             cGGTTTTTTACGGCTATCTTCTTTTAAAAATTTTaa                         >  1:6663229/1‑36 (MQ=18)             cGGTTTTTTACGGCTATCTTCTTGAAAAGATTATaa                         <  1:4065611/36‑1 (MQ=255)              ggTTTTTTACGGCTATCTTCTTTAAAAGATAAAaaa                        >  1:4642692/1‑36 (MQ=21)               gTTTTTTACGGCTATTTTCTTTAAAAGATTTTAAAg                       >  1:3580962/1‑36 (MQ=21)               gTTTTTTACGGCTAACTTTTTTTAAAGATTATaaaa                       >  1:1404639/1‑35 (MQ=17)               gTTTTTTACGGCTAACTTCTTTAAAAGATTATaaaa                       >  1:2993144/1‑35 (MQ=25)                 tttttACGGCTATCTTCTTTAAAAGATTATAAAAAt                     >  1:6620611/1‑36 (MQ=25)                  ttttACGGCTATCTTCTTGAAAAGATTATAAAGATg                    <  1:3009678/36‑1 (MQ=255)                  ttttACGGCTATCTTCTTGAAAAGATTATAAAGATg                    <  1:1600610/36‑1 (MQ=255)                  tttaACGACTATCTTCTTGAAAAGATTATAAAGATg                    <  1:3255305/36‑1 (MQ=25)                     tACGGCTATCTTCTTGAAAAGATTATAAAGATGCGt                 >  1:4547907/1‑36 (MQ=255)                       cGGCTATCTTCTTGAAAAGATTATAAAGATGCGttt               <  1:7014456/36‑1 (MQ=255)                        ggCTATCTTCTTGAAAAGATTATAAAGATGCGtttt              >  1:323520/1‑36 (MQ=255)                         gCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTa             <  1:6307575/36‑1 (MQ=255)                          cTATCTTCTTGAAAAGATTATAAAGATGCGTTTTaa            >  1:1526438/1‑36 (MQ=255)                           tATCTTCTTGAAAAGATTATAAAGATGCGTTTTAAc           <  1:3634026/36‑1 (MQ=255)                            aTCTTCTTGAAAAGATTATAAAGATGCGTTTTAAcc          >  1:4558092/1‑36 (MQ=255)                            aTCTTCTTGAAAAGATTATAAAGATGCGTTTTAAcc          <  1:7294596/36‑1 (MQ=255)                            aTCTTCGTGAAAAGATTATAAAGATGCGTTTTAAcc          <  1:7204453/36‑1 (MQ=255)                              cttcttGAAAAGATTATAAAGATGCGTTTTAACCGt        <  1:3124643/36‑1 (MQ=255)                                tcttGAAAAGATTATAAAGATGCGTTTTAACTGttt      >  1:1150970/1‑34 (MQ=25)                                tcttAAAAAGATTATAAAGATGCGTTTTAACCGTAt      <  1:779668/36‑1 (MQ=255)                                 cttGAAAAGATTATAAAGATGCGTTTTAACCGTAtt     >  1:6563874/1‑36 (MQ=255)                                 cttGAAAAGATTATAAAGATGCGTTTTAACCGTAtt     >  1:4861279/1‑36 (MQ=255)                                  ttGTAAAGATTATAAAGATGCGTTTTAACCGTAttt    <  1:5989438/36‑1 (MQ=255)                                  ttGAAAAGATTATAAAGATGCGTTTTAACCGTAttt    >  1:903028/1‑36 (MQ=255)                                   tGAAAAGATTATAAAGATGCGTTTTAACCGTAtttt   >  1:2484558/1‑36 (MQ=255)                                   tGAAAAGATTATAAAGATGCGTTTTAACCGTAtttt   >  1:1872664/1‑36 (MQ=255)                                    gAAAAGATTATAAAGATGCGTTTTAACCGTATTTTc  <  1:6416477/36‑1 (MQ=255)                                    |                                    ACCGCTTGTGTCCGGTTTTTTACGGCTATCTTCTTGAAAAGATTATAAAGATGCGTTTTAACCGTATTTTC  >  REL606/1117146‑1117216 Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 16 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG