Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 831,492 T→C 15.1% intergenic (+465/‑51) narK → / → narG nitrate/nitrite transporter/nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE831,4920TC15.1% 75.3 / 4.6 46intergenic (+465/‑51)narK/narGnitrate/nitrite transporter/nitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (26/13);  new base C (4/3);  total (30/16)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATT  >  minE/831427‑831553
                                                                 |                                                             
ccccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGc                                                            <  1:1105379/69‑1 (MQ=255)
ccccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGc                                                            <  1:839195/69‑1 (MQ=255)
ccccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGc                                                            <  1:765831/69‑1 (MQ=255)
ccccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGc                                                            <  1:2621705/69‑1 (MQ=255)
ccccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGc                                                            <  1:1320065/69‑1 (MQ=255)
ccccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGc                                                            <  1:1157931/69‑1 (MQ=255)
ccccATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCATAAACATTAGc                                                            <  1:770254/69‑1 (MQ=255)
      cACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCaa                                                          >  1:669979/1‑65 (MQ=255)
      cACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCaa                                                          >  1:2562072/1‑65 (MQ=255)
      cACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCaa                                                          >  1:2242704/1‑65 (MQ=255)
      cACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCaa                                                          >  1:2418872/1‑65 (MQ=255)
      cACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCaa                                                          >  1:3339636/1‑65 (MQ=255)
       aCTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCg                                                     >  1:3360360/1‑69 (MQ=255)
                     tCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGa                                                    >  1:1793583/1‑56 (MQ=255)
                     tCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAggg                                        >  1:2767337/1‑68 (MQ=255)
                     tCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAggg                                        >  1:322852/1‑68 (MQ=255)
                     tCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAggg                                        >  1:1405530/1‑68 (MQ=255)
                     tCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAggg                                        >  1:671993/1‑68 (MQ=255)
                     tCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAggg                                        >  1:2373698/1‑68 (MQ=255)
                      cAATTCCCACGCTGTTTCAGCGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:1283298/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:3402179/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:331870/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:569375/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:1100638/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:3153206/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:2947076/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:1335786/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:2717475/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:2626410/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:1250225/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:2474365/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:2275491/1‑58 (MQ=255)
                      cAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGAttt                                                 >  1:2077029/1‑58 (MQ=255)
                           cccACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATCAGCAATGTCGATTTATCAGAGGGCCGACAg                                 >  1:1338329/1‑69 (MQ=255)
                           cccACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATCAGCAATGTCGATTTATCAGAGGGCCGACAg                                 >  1:1514420/1‑69 (MQ=255)
                           cccACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATCAGCAATGTCGATTTATCAGAGGGCCGACAg                                 >  1:2261098/1‑69 (MQ=255)
                           cccACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATCAGCAATGTCGATTTATCAGAGGGCCGACAg                                 >  1:749773/1‑69 (MQ=255)
                                                       ccTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTa      <  1:893024/68‑1 (MQ=255)
                                                       ccTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTa      <  1:1545139/68‑1 (MQ=255)
                                                         ttAAACATCAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAAt   <  1:2250242/69‑1 (MQ=255)
                                                         ttAAACATCAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAAt   <  1:3061013/69‑1 (MQ=255)
                                                         tcaaaCATCAGCAATGTCGATTTATCAGAGGCCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAAt   <  1:3230210/67‑1 (MQ=255)
                                                          tAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAAtt  <  1:194405/69‑1 (MQ=255)
                                                          tAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAAtt  <  1:1586028/69‑1 (MQ=255)
                                                          tAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAAtt  <  1:3100885/69‑1 (MQ=255)
                                                          tAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAAtt  <  1:298327/69‑1 (MQ=255)
                                                                 |                                                             
CCCCATCACTCTTGATCGTTATCAATTCCCACGCTGTTTCAGAGCGTTACCTTGCCCTTAAACATTAGCAATGTCGATTTATCAGAGGGCCGACAGGCTCCCACAGGAGAAAACCGATGAGTAAATT  >  minE/831427‑831553

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: