Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,286,699 C→A 93.1% 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→A93.1% 3.6 ‑2.0 ‑0.249intergenic (+23/‑160)narI/ECB_01206nitrate reductase 1, gamma (cytochrome b(NR)) subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  new base (15/30):  ref base (3/1):  total (18/31)
Fisher's exact test strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 6.32e-01

CTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGT  >  REL606/1286666‑1286731
                                 |                                
ctagTCACTAAGCGAATTTTAGTTCACATAGACACt                                <  1:2089961/33‑1 (MQ=38)
cTCGTCACTAAGCGAATTTTAGTTCACATAGACACt                                >  1:3848927/1‑36 (MQ=255)
cTCGTCACTAAGCGAATTTTAGTTCACATAGACACt                                <  1:1746794/36‑1 (MQ=255)
   gTCACTAAGCGAATTTTAGTTCACATAGACACTGCt                             <  1:5040305/36‑1 (MQ=255)
   gTCACTAAGCGAATTTTAGTTCACATAGACACTGCt                             <  1:837865/36‑1 (MQ=255)
   gTCACTAAGCGAATTTTAGTTCACATAGACACTGCt                             <  1:2971519/36‑1 (MQ=255)
     cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                           <  1:5892558/36‑1 (MQ=255)
     cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                           <  1:6743314/36‑1 (MQ=255)
     cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                           <  1:3136583/36‑1 (MQ=255)
     cACTAAGCGAATTTTAGTTCACATAGACACTGCTTc                           <  1:1693917/36‑1 (MQ=255)
      aCTAAGTGAATTTTAGTTCACATAGACACTGCTTcg                          <  1:6869881/36‑1 (MQ=38)
       cTAAGCGAATTTTAGTTCACATAGACACTGCTTcgg                         <  1:3689718/36‑1 (MQ=255)
        tAAGCGAATTTTAGTTCACATAGACACTGCTTcggc                        >  1:4259560/1‑36 (MQ=255)
        tAAGCGAATTTTAGTTCACATAGACACTGCTTcggc                        <  1:7922416/36‑1 (MQ=255)
         aaGCGAATTTTAGTTCACATAGACACTGCTTcggcg                       <  1:1817074/36‑1 (MQ=255)
         aaGCGAATTTTAGTTCACATAGACACTGCTTcggcg                       >  1:6208847/1‑36 (MQ=255)
          aGCGAATTTTAGTTCACATAGACCCTGCTTcggcgg                      >  1:4259022/1‑36 (MQ=255)
           gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                     >  1:1468870/1‑36 (MQ=255)
           gCGAATTTTAGTTCACATAGACACTGCTTCGGCggg                     <  1:6530717/36‑1 (MQ=255)
            cGAATTTTAGTTCACATAGACCCTGCTTCGGCgggg                    <  1:4898156/36‑1 (MQ=255)
             gAATTTTAGTTCACATAGACACTGCTTCGGCGGGGt                   <  1:2517595/36‑1 (MQ=255)
              aaTTTTAGTTCACATAGACACTGCTTCGGCGGGGtt                  <  1:7665839/36‑1 (MQ=255)
              aaTTTTAGTTCACATAGACACTGCTTCGGCGGGGtt                  <  1:4755305/36‑1 (MQ=255)
               aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                 <  1:2952910/36‑1 (MQ=255)
               aTTTTAGTTCACATAGACACTGCTTCGGCGGGGttt                 >  1:1082475/1‑36 (MQ=255)
                 tttAGTTCACATAGACCCTGCTTCGGCGGGGttttt               >  1:3178623/1‑36 (MQ=255)
                 tttAGTTCACATAGACACTGCTTCGGCGGGGttttt               <  1:7409515/36‑1 (MQ=255)
                  ttAGTTCACATAGACCCTGCTTCGGCGGGGtttttt              >  1:3950724/1‑36 (MQ=255)
                  ttAGTTCACATAGACACTGCTTCGGCGGGGtttttt              >  1:7448567/1‑36 (MQ=255)
                   tAGTTCACATAGACACTGCTTCGGCGGGGttttttt             <  1:8157667/36‑1 (MQ=255)
                   tAGTTCACATAGACACTGCTTCGGCGGGGttttttt             <  1:1068933/36‑1 (MQ=255)
                      ttCACATAGACACTGCTTCGGCGGGGTTTTTTTAtt          >  1:1378409/1‑35 (MQ=255)
                      ttCACATAGACACTGCTTCGGCGGGGTTTTTTTAtt          >  1:4642612/1‑35 (MQ=255)
                      ttCACATAGACACTGCTGCGGCGGGGTTTTTTTATg          <  1:4299734/36‑1 (MQ=37)
                       tCACATAGACACTGCTTCGGCGGGGTTTTTTTATgg         <  1:2489385/36‑1 (MQ=255)
                        cacaTAGACACTGCTTCGGCGGGGTTTTTTTATggg        <  1:960377/36‑1 (MQ=255)
                         acaTAGACACTGCTTCGGCGGGGTTTTTTTATTGGc       >  1:1421669/1‑36 (MQ=25)
                         acaTAGACACTGCTTCGGCGGGGTTTTTTTATGGGc       <  1:4966621/36‑1 (MQ=255)
                          caTAGACACTGCTTCGGCGGGGTTTTTTTATTGGCa      >  1:7308201/1‑36 (MQ=25)
                          caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa      >  1:481629/1‑36 (MQ=255)
                          caTAGACACTGCTTCGGCGGGGTTTTTTTATGGGCa      <  1:7870263/36‑1 (MQ=255)
                           aTAGACACTGCTTCGGCGGGGTTTTTTTATTTGCAc     >  1:8169907/1‑36 (MQ=18)
                           aTAGACACTGCTTCGGCGGGGTTTTTTTATTGGAAc     >  1:2613685/1‑36 (MQ=18)
                            tAGACACTGCTTCGGCGGGGTTTTTTTTTGGGCagg    >  1:2482278/1‑34 (MQ=21)
                            tAGACACTGCTTCGGCGGGGTTTTTTTATTGGCACg    >  1:7961363/1‑36 (MQ=25)
                              gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:350876/36‑1 (MQ=255)
                              gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:6212519/36‑1 (MQ=255)
                              gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:4831082/36‑1 (MQ=255)
                              gACACTGCTTCGGCGGGGTTTTTTTATGGGCACGGt  <  1:5604073/36‑1 (MQ=255)
                                 |                                
CTCGTCACTAAGCGAATTTTAGTTCACATAGACCCTGCTTCGGCGGGGTTTTTTTATGGGCACGGT  >  REL606/1286666‑1286731

Base quality scores:  ATCG  < 3 ≤  ATCG  < 17 ≤  ATCG  < 26 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG