Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	1,286,699	C→A	intergenic (+23/‑160)	narI/ECB_01206	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein

Read alignment evidence...
	seq id	position		change	freq	score	reads	annotation	genes	product
*	REL606	1,286,699	0	C→A	100.0%	152.7	52	intergenic (+23/‑160)	narI/ECB_01206	nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (18/32):  ref base (0/0):  total (18/34)

CTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGT  >  REL606/1286666‑1286731                                  |                                 cTCGTCACTAAGCGAATTTTAGTTCACATAGACACt                                <  1:1338398/36‑1 (MQ=255) cTCGTCACTAAGCGAATTTTAGTTCACATAGACACt                                <  1:432456/36‑1 (MQ=255) cTCGTCACTAAGCGAATTTTAGTTCACATAGACACt                                <  1:3984244/36‑1 (MQ=255)   cGTCACTAAGCGAATTTTAGTTCACATAGACACTGc                              <  1:7172464/36‑1 (MQ=255)    gTCACTAAGCGAATTTTAGTTCACATAGACACTGCt                             >  1:3568780/1‑36 (MQ=255)     tCACTAAGCGAATTTTAGTTCACATAGACACTGCtt                            <  1:7321552/36‑1 (MQ=255)     tCACTAAGCGAATTTTAGTTCACATAGACACTGCtt                            >  1:5577654/1‑36 (MQ=255)      cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                           <  1:6993915/36‑1 (MQ=255)      cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                           >  1:4844848/1‑36 (MQ=255)      cACTAAGAGAATTTTAGTTCACATAGACACTGCTTc                           >  1:6134185/1‑36 (MQ=25)       aCTAAGCGAATTTTAGTTCACATAGACACTGCTTcg                          >  1:4313312/1‑36 (MQ=255)       aCTAAGCGAATTTTAGTTCACATAGACACTGCTTcg                          >  1:2179331/1‑36 (MQ=255)       aCTAAGCGAATTTTAGTTCACATAGACACTGCTTcg                          <  1:100917/36‑1 (MQ=255)         tAAGCGAATTTTAGTTCTCATAGACACTGCTTcggc                        <  1:1630250/36‑1 (MQ=25)          aaGCGAATTTTAGTTCACATAGACACTGCTTcggcg                       <  1:4501687/36‑1 (MQ=255)          aaGCGAATTTTAGTTCACATAGACACTGCTAcagcg                       >  1:3326879/1‑36 (MQ=21)           aGCGAATTTTAGTTCACATAGACACTGCCTcggcgg                      >  1:7079706/1‑36 (MQ=25)            gCGTATTTTAGTTCACATAGACACTGCTTCGGCggg                     <  1:6120965/36‑1 (MQ=255)            gCGTATTTTAGTTCACATAGACACTGCTTCGGCggg                     <  1:4165360/36‑1 (MQ=255)            gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                     >  1:2703177/1‑36 (MQ=255)            gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                     <  1:2576795/36‑1 (MQ=255)            gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                     <  1:6295019/36‑1 (MQ=255)             cGAATTTTAGTTCACATAGACACTGCTTCGGCgggg                    <  1:7318502/36‑1 (MQ=255)              gtaTTTTAGTTCACATAGACACTGCTTCGGCGGGGt                   <  1:454/34‑1 (MQ=37)              gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                   >  1:7313874/1‑36 (MQ=255)              gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                   <  1:4536138/36‑1 (MQ=255)              gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                   <  1:2646949/36‑1 (MQ=255)               aaTTTTAGTTCGCATAGACACTGCTTCGGCGGGGtt                  <  1:4166011/36‑1 (MQ=25)               aaTTTTAGTTCACATAGACACTGCTTCGGCGGGGtt                  <  1:1383104/36‑1 (MQ=255)               aaTTTTAGTTCACATAGACACTGCTTCGCCGGGGtt                  >  1:6672761/1‑36 (MQ=25)                aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                 >  1:690117/1‑36 (MQ=255)                 ttttAGTTCACATAGACGCTGCTTCGGCGGGGtttt                <  1:6382735/36‑1 (MQ=255)                 ttttAGTTCACATAGACACTGCTTCGGCGGGGtttt                <  1:974897/36‑1 (MQ=255)                 ttttAGTTCACATAGACACTGCTTCGGCGGGGtttt                <  1:187094/36‑1 (MQ=255)                  tttAGTTCACATAGACGCTGCTTCGGCGGGGttttt               <  1:3639647/36‑1 (MQ=255)                   ttAGTTCACATAGACACTGCTTCGGCGGGGtttttt              <  1:7060926/36‑1 (MQ=255)                     aGTTCACATAGACACTGCTTCGGCGGGGTTTTTTTa            <  1:4948525/36‑1 (MQ=255)                      gTTCGCATAGACACTGCTTCGGCGGGGTTTTTTTAt           <  1:6294349/36‑1 (MQ=21)                       ttCACATAAACACTGCTTCGGCTGGGTTTTTTTATg          >  1:3145471/1‑36 (MQ=18)                        tCACATAGACACTGCTTCGGCGGGGTTTTTTTATgg         <  1:647771/36‑1 (MQ=255)                         cacaTAGACACTGCTTCGGCGGGGTTTTTTTATggg        <  1:6474932/36‑1 (MQ=255)                           caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa      <  1:6326730/36‑1 (MQ=255)                           caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa      <  1:5845514/36‑1 (MQ=255)                           caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa      >  1:3739890/1‑36 (MQ=255)                              aGACACTGCTTCGGCGGGGTTTTTTTATTGGCTCgg   >  1:5897368/1‑36 (MQ=18)                              aGACACTGCTTCGGCGGGGTTTTTTTATGGGCTCgg   >  1:4887121/1‑36 (MQ=21)                              aGACACTGCTTCGGCGGGGTTTTTTTATGGGCGCgg   >  1:7133379/1‑36 (MQ=21)                              aGACACTGCTTCGGCGGGGTTTTTTTATGGGCACgg   >  1:5220990/1‑36 (MQ=255)                               gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:5405476/36‑1 (MQ=255)                               gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:4821635/36‑1 (MQ=255)                               gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:620856/36‑1 (MQ=255)                               gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:2190891/36‑1 (MQ=255)                                  |                                 CTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGT  >  REL606/1286666‑1286731 Base quality scores:  ATCG  < 3 ≤  ATCG  < 9 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG