Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,286,699 C→A 100% 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
*REL6061,286,6990C→A100.0% 69.3 20intergenic (+23/‑160)narI/ECB_01206nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (5/15):  ref base (0/0):  total (5/15)

CTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGT  >  REL606/1286666‑1286731
                                 |                                
cTCGTCACTAAGCGAATTTTAGTTCACATAGACACt                                <  1:5852543/36‑1 (MQ=255)
atCGTCAATAAGCGAATTTTAGTTCACATAGACACt                                <  1:1814977/35‑1 (MQ=37)
  cGTCACTAAGCGAATTTTAGTTCACATAGACACTGc                              <  1:3986453/36‑1 (MQ=255)
    tCACTAAGCGAATTTTAGTTCACATAGACACTGCtt                            >  1:704798/1‑36 (MQ=255)
    tCACTAAGCGAATTTTAGTTCACATAGACACTGCtt                            >  1:353188/1‑36 (MQ=255)
    tCACTAAGCGAATTTTAGTTCACATAGACACTGCtt                            <  1:236757/36‑1 (MQ=255)
    tCACTAAGCGAATTTTAGTTCACATAGACACTGCtt                            <  1:1289053/36‑1 (MQ=255)
     cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                           <  1:1626208/36‑1 (MQ=255)
         aaGCGAATTTTAGTTCACATAGACACTGCTTcggcg                       <  1:2479826/36‑1 (MQ=255)
            cGAATTTTAGTTCACATAGACACTGCTTCGGCgggg                    <  1:3594632/36‑1 (MQ=255)
            cGAATTTTAGTTCACATAGACACTGCTTCGGCgggg                    <  1:4748724/36‑1 (MQ=255)
             gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                   <  1:3241774/36‑1 (MQ=255)
               aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                 >  1:4563735/1‑36 (MQ=255)
               aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                 >  1:1265691/1‑36 (MQ=255)
                 tttAGTTCACATAGACACTGCTTCGGCGGGGttttt               >  1:38932/1‑36 (MQ=255)
                    aGTTCACATAGACACTGCTTCGGCGGGGTTTTTTTa            <  1:2731140/36‑1 (MQ=255)
                      ttCACATAGACACTGCTTCGGCGGGGTTTTTTTATg          <  1:3405869/36‑1 (MQ=255)
                           aTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCAc     <  1:4521072/36‑1 (MQ=255)
                              gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:2603158/36‑1 (MQ=255)
                              gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:5298025/36‑1 (MQ=255)
                                 |                                
CTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGT  >  REL606/1286666‑1286731

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG