Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	411,493	0	T→C	30.0%	10.3	40	M191V (ATT→GTT)	thiJ	hypothetical protein
Reads supporting (aligned to +/- strand):  new base (1/11):  ref base (16/12):  total (17/23)

TTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAT  >  REL606/411458‑411528                                    |                                    ttGCCGCATCCGACACTACTCGTAATAATTATAAAt                                     >  1:6317599/1‑36 (MQ=255)    ccGCATCCGACATTACTCGTAATAATTATAAATccc                                  <  1:523025/36‑1 (MQ=255)    ccGCATCCGACACTACTCGTAATAATTATAAATccc                                  >  1:5888478/1‑36 (MQ=255)     cGCATGCGACACTACTCGTAATAATTATAAATCCCt                                 <  1:4026509/36‑1 (MQ=255)      gaaTCCGACACTACTCGTAATAATTATAAATCCCTg                                <  1:1542195/34‑1 (MQ=255)        aTCCGACACTACTCGTAATAATTATAAATCCCTgcc                              >  1:4127672/1‑36 (MQ=255)           cGACACTACTCGTAATAATTATAAATCCCTgccgcc                           <  1:411202/36‑1 (MQ=255)           cGACACTACTCGTAATAATTATAAACCCCCgccgcc                           >  1:505525/1‑36 (MQ=25)              cacTACTCGTAATAATTATAAATCCCTGCCGCCATc                        <  1:5143978/36‑1 (MQ=255)               actactCGTAATAATTATAAATCCCTGCCGCCATCa                       >  1:6362682/1‑36 (MQ=255)               actactCGTAATAATTATAAATCCCTGCCGCCATCa                       <  1:5145814/36‑1 (MQ=255)                ctactCGTAATAATTATAAATCCCTGCCGCCATCAc                      <  1:2750670/36‑1 (MQ=255)                ctactCGTAATAATTATAAATCCCTGCCGCCATCAc                      >  1:5273358/1‑36 (MQ=255)                 tactCGTAATAATTATAAATCCCTGCCGCCATCAcc                     >  1:1955624/1‑36 (MQ=255)                  tctCGTAATAATTATAAATCCCTGCCGCCATCACCa                    <  1:388444/35‑1 (MQ=255)                   ctCGTAAAAATTAAAAACCCCTGCCGCCATCACCAg                   <  1:6652967/36‑1 (MQ=21)                     cGTAATTATTAAAAATCCCTGCCGCCATCACCAGtt                 <  1:6091211/36‑1 (MQ=25)                       taataaAAATAAACCCCTGCCGCCATCACCAGTtgt               <  1:6802973/36‑1 (MQ=21)                        aataatTATAAATCCCTGCCGCCATCACCAGTtgtg              >  1:1538199/1‑36 (MQ=255)                         ataatTATAAATCCCTCCCGCCATCACCAGTTGTGa             <  1:1868360/36‑1 (MQ=255)                          taatTATAAATCCCTGCCGCCATCACCAGTTGTGAt            >  1:1696398/1‑36 (MQ=255)                          taatTATAAATCCCTGCCGCCATCACCAGTTGTGAt            >  1:6637432/1‑36 (MQ=255)                          taatTATAAACCCCTGCCGCCATCGCCAGTTGTGAt            <  1:5090803/36‑1 (MQ=25)                           aatTATAAATCCCTGCCGCCATCACCANTNTTNATg           >  1:2226424/1‑36 (MQ=25)                           aatTATAAATCCCTGCCGCCATCACCAGTTGTGATg           >  1:3055060/1‑36 (MQ=255)                           aatTAGAAACCCCTGCCGCCATCACCAGTTGTGATg           <  1:3311961/36‑1 (MQ=25)                            ataataAACCCCCGCCGCCATCACCAGTTGTGATGc          <  1:3388362/33‑1 (MQ=21)                            atTATAAATCCCTGCCGCCATCACCAGTTGTGATGc          >  1:2891681/1‑36 (MQ=255)                            atTATAAATCCCTGCCGCCATCACCAGTTGTGATGc          <  1:2623670/36‑1 (MQ=255)                             ttATAAATCCCTGCCGCCATCACCAGTTGTGATGcc         >  1:5987408/1‑36 (MQ=255)                              tataAATCCCTGCCGCCATCACCAGTTGTGATGCCa        >  1:2352816/1‑36 (MQ=255)                              tataAACCCCTGCCGCCATCACCAGTTGTGATGCCa        <  1:6817023/36‑1 (MQ=255)                              tagaaaCCCCTGCCGCCATCACCAGTTGTGATGCCa        <  1:7395077/33‑1 (MQ=25)                              taaaaaCCCCTGCCGCCATCACCAGTTGTGATGCCa        <  1:1246166/33‑1 (MQ=25)                               ataAATCCCTGCCGCCATCACCAGTTGTGATGCCAc       >  1:6987872/1‑36 (MQ=255)                                taAATCCCTGCCGCCATCACCAGTTGTGATGCCACt      <  1:3350787/36‑1 (MQ=255)                                aaaaCCCCTGCCGCCATCACCAGTTGTGATGCCACt      <  1:1879913/35‑1 (MQ=255)                                 aaaTCCCTGCCGCCATCACCAGTTGTGATGCCACtt     >  1:5692191/1‑36 (MQ=255)                                 aaaCCCCTGCCGCCATCACCAGTTGTGATGCCACtt     <  1:6104674/36‑1 (MQ=255)                                 aaaCCCCTGCCGCCATCACCAGTTGTGATGCCACtt     <  1:4027652/36‑1 (MQ=255)                                  aaTCCCTGCCGCCATCACCAGTTGTGATGCCACTTc    >  1:4643068/1‑36 (MQ=255)                                    tCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAt  >  1:7053620/1‑36 (MQ=255)                                    tCCCTGCCGCCATCACCAGTTGTGATGCCACTCCAt  >  1:4301324/1‑36 (MQ=255)                                    tCCCTGCCGCCATCACCAGTTGTGATGCCACTACAt  >  1:4787699/1‑36 (MQ=255)                                    |                                    TTGCCGCATCCGACACTACTCGTAATAATTATAAATCCCTGCCGCCATCACCAGTTGTGATGCCACTTCAT  >  REL606/411458‑411528 Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 18 ≤  ATCG  < 29 ≤  ATCG  < 32 ≤  ATCG