Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 833,210 T→C 9.4% R556R (CGT→CGC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE833,2100TC 9.4% 204.0 / 11.7 74R556R (CGT→CGCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (29/38);  new base C (2/5);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATG  >  minE/833142‑833277
                                                                    |                                                                   
cACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTgg                                                                   <  1:1031779/71‑1 (MQ=255)
  cGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                                                                 >  1:3584905/1‑71 (MQ=255)
    gTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTg                                                               <  1:1479669/71‑1 (MQ=255)
      cGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGAt                                                             >  1:3966255/1‑71 (MQ=255)
       gctCGATGATTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATc                                                            <  1:1547010/69‑1 (MQ=255)
        ctCGATGATTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCa                                                           <  1:3507793/70‑1 (MQ=255)
          cGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                                                          >  1:743291/1‑70 (MQ=255)
          cGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                                                          <  1:838513/70‑1 (MQ=255)
          cGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                                                          >  1:1957518/1‑70 (MQ=255)
           gatgatTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCaa                                                          <  1:3581984/69‑1 (MQ=255)
              gatTATCCTCGGCGCAGGTGTTAACCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGc                                                     >  1:957155/1‑71 (MQ=255)
                ttATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGa                                                              >  1:4110667/1‑60 (MQ=255)
                  aTCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                                                                 >  1:3046272/1‑55 (MQ=255)
                  aTCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAt                                                 >  1:2868683/1‑71 (MQ=255)
                     gtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTgg                                                                   <  1:2633458/50‑1 (MQ=255)
                     gtcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                                                                  <  1:3842512/51‑1 (MQ=255)
                      tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTg                                                                    <  1:1881647/48‑1 (MQ=255)
                      tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTg                                                                    <  1:1120413/48‑1 (MQ=255)
                      tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTg                                                                    <  1:2505111/48‑1 (MQ=255)
                      tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                                                                  <  1:382336/50‑1 (MQ=255)
                      tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                                                                  <  1:3373450/50‑1 (MQ=255)
                      tcgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                             >  1:3716220/1‑71 (MQ=255)
                       cgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTc                                                                 >  1:598088/1‑50 (MQ=255)
                       cgGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCt                                                                >  1:2940499/1‑51 (MQ=255)
                        ggTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                                                                  >  1:2095725/1‑48 (MQ=255)
                        ggTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCt                                                                >  1:103783/1‑50 (MQ=255)
                           gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCaa                                                          <  1:2675433/53‑1 (MQ=255)
                           gCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCg                                         >  1:3018486/1‑70 (MQ=255)
                            cGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTg                                                               <  1:4102981/47‑1 (MQ=255)
                            cGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTg                                                               <  1:1163247/47‑1 (MQ=255)
                                gCTGAACCACTGGTATCACCTCGATATGAACTATCGTggg                                                                  <  1:3460443/40‑2 (MQ=255)
                                gCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGt                                                                  <  1:539666/40‑1 (MQ=255)
                                  ttaaCCACTGGTATCACATGGACATGAACTACCGCGGGATGATCAATATGCTGATCTTCTg                                           <  1:2124134/59‑1 (MQ=255)
                                   taaCCACTGGTATCACATGGACATGAACTACCGCGGGATGATc                                                            <  1:750418/42‑1 (MQ=255)
                                    aaCCACTGGTATCACCTCGATATGAACTATCGTGGt                                                                  <  1:4112180/36‑1 (MQ=255)
                                    aaCCACTGGTATCACCTCGATATGAACTATCGTGGTCt                                                                >  1:3871801/1‑38 (MQ=255)
                                    aaCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCt                                      <  1:2514496/64‑1 (MQ=255)
                                    aaCCACTGGTATCACATGGACATGAACTACCGCGGGATGATCaata                                                        >  1:2231859/1‑44 (MQ=12)
                                     aCCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                              <  1:860394/71‑1 (MQ=255)
                                         cTGGTATCACCTCGATTTGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTc                             <  1:1862111/68‑1 (MQ=255)
                                         cTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCaga                          <  1:945602/71‑1 (MQ=255)
                                          tGGTATCACCTCGATATGAACTATCGTGGTCTGATCa                                                           <  1:648578/37‑1 (MQ=255)
                                          tGGTATCACCTCGATATGAACTATCGTGGTCTGATCa                                                           <  1:977748/37‑1 (MQ=255)
                                          tGGTATCACCTAGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCagag                         >  1:4014660/1‑71 (MQ=255)
                                            gTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGAtttt                                              >  1:2034287/1‑48 (MQ=255)
                                            gTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCagag                         <  1:3434237/69‑1 (MQ=255)
                                            gTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCg                       <  1:509547/71‑1 (MQ=255)
                                            gTATCACCTCGATATGAACTATCGTGGTATGATCAACATGCTGAtttt                                              >  1:2822667/1‑48 (MQ=255)
                                             tATCACCTCGATATGAACTATCGTGGTCTGATCaa                                                          <  1:3568320/35‑1 (MQ=255)
                                                 aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTgtgt                                  >  1:1628008/1‑55 (MQ=255)
                                                 aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGggc                  <  1:315065/71‑1 (MQ=255)
                                                 aCCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGggc                  >  1:684835/1‑71 (MQ=255)
                                                   cTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTc                                             >  1:1870835/1‑42 (MQ=255)
                                                   cTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGggcgg                >  1:4045942/1‑71 (MQ=255)
                                                     cGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCt              >  1:3053443/1‑71 (MQ=255)
                                                        tatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCg                       <  1:348373/59‑1 (MQ=255)
                                                        tatGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTggg           >  1:3839502/1‑71 (MQ=255)
                                                          tGAACTATCGTGGTCTGATCAACATGCTGATTTTCt                                            <  1:980404/36‑1 (MQ=255)
                                                          tGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGgcg         <  1:3991403/71‑1 (MQ=255)
                                                          tGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGgcg         <  1:1466397/71‑1 (MQ=255)
                                                              cTATCGTGGTCTGATCAACATGCTGATTTTCTGCgg                                        <  1:3909963/36‑1 (MQ=255)
                                                              cTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGt                              <  1:3546475/46‑1 (MQ=255)
                                                              cTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCag                           <  1:3397444/49‑1 (MQ=255)
                                                              cTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTg             >  1:2107738/1‑63 (MQ=255)
                                                              cTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCAc      >  1:2678059/1‑70 (MQ=255)
                                                               tATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTggg           <  1:1948615/64‑1 (MQ=255)
                                                                 tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCtgtg                                   <  1:3384509/38‑1 (MQ=255)
                                                                 tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggg                     >  1:2073244/1‑52 (MQ=255)
                                                                 tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg  <  1:2088908/71‑1 (MQ=255)
                                                                 tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg  <  1:801615/71‑1 (MQ=255)
                                                                 tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg  <  1:843544/71‑1 (MQ=255)
                                                                 tCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg  >  1:1648506/1‑71 (MQ=255)
                                                                   gTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCggg                     >  1:1695867/1‑50 (MQ=255)
                                                                   gTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATg  >  1:1082347/1‑69 (MQ=255)
                                                                    |                                                                   
CACGGTCGTTCGATGATTATCGTCGGTGCGGGGCTGAACCACTGGTATCACCTCGATATGAACTATCGTGGTCTGATCAACATGCTGATTTTCTGCGGCTGTGTCGGTCAGAGCGGGGGCGGCTGGGCGCACTATG  >  minE/833142‑833277

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