Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,124 A→G 6.9% A194A (GCA→GCG narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1240AG 6.9% 162.0 / 3.8 58A194A (GCA→GCGnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base G (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

CCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCT  >  minE/832054‑832194
                                                                      |                                                                      
ccgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCa                                                                        <  1:3948961/71‑1 (MQ=255)
     tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa                                                                       >  1:1465581/1‑67 (MQ=255)
     tCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc                                                                   >  1:483044/1‑71 (MQ=255)
      cTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTc                                                                   <  1:1805004/70‑1 (MQ=255)
       tAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGa                                                                 >  1:2173058/1‑71 (MQ=255)
       tAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGa                                                                 >  1:1755868/1‑71 (MQ=255)
         aCGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATg                                                               >  1:479996/1‑71 (MQ=255)
          cGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATgg                                                              <  1:2944115/71‑1 (MQ=255)
            tttACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt                                                            >  1:3686150/1‑71 (MQ=255)
              tACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGttt                                                           <  1:130029/70‑1 (MQ=255)
                cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATg                                                                     >  1:2676481/1‑58 (MQ=255)
                cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCt                                                         <  1:1578860/70‑1 (MQ=255)
                cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCt                                                         <  1:640374/70‑1 (MQ=255)
                  ccATCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTAc                                                      <  1:1314850/71‑1 (MQ=255)
                    aTCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGc                                                    <  1:1493212/71‑1 (MQ=255)
                     tCAAAACCTACGGCCCGGATCGCGTTGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGt                                                             <  1:2299396/61‑1 (MQ=255)
                      cAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCAt                                                  >  1:1310227/1‑71 (MQ=255)
                       aaaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCaa                                                                       >  1:3475153/1‑49 (MQ=255)
                       aaaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATc                                                 >  1:1724244/1‑71 (MQ=255)
                       aaaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATc                                                 >  1:694688/1‑71 (MQ=255)
                        aaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCg                                                <  1:2911752/71‑1 (MQ=255)
                           aCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATgg                                                              >  1:4104735/1‑54 (MQ=255)
                            cTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTg                                            <  1:1575124/71‑1 (MQ=255)
                                 gCCCGTACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGc                                       >  1:2774026/1‑71 (MQ=255)
                                  cccGGACCGTGTTGCTGGTTTCCCCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGt                                             >  1:1361467/1‑63 (MQ=255)
                                    cGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAt                                    >  1:221494/1‑71 (MQ=255)
                                       aCCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctc                                 <  1:537180/71‑1 (MQ=255)
                                       aCCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct                                  >  1:1812175/1‑70 (MQ=255)
                                        ccGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctct                                >  1:2979674/1‑71 (MQ=255)
                                         cGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACg                                                     <  1:2179863/49‑1 (MQ=255)
                                           tgtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGa                                                                 >  1:969154/1‑35 (MQ=255)
                                           tgtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGc                             >  1:3946193/1‑71 (MQ=255)
                                            gtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt                                                            <  1:84032/39‑1 (MQ=255)
                                            gtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGtt                                                            <  1:4013773/39‑1 (MQ=255)
                                             ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTa                                                       >  1:1173899/1‑43 (MQ=255)
                                             ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                            >  1:3406514/1‑70 (MQ=255)
                                             ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                            <  1:1586102/70‑1 (MQ=255)
                                             ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                            >  1:459835/1‑70 (MQ=255)
                                             ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                            >  1:3746947/1‑70 (MQ=255)
                                                 tGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                            <  1:2733639/66‑1 (MQ=255)
                                                  ggTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCt                                                         <  1:2113566/36‑1 (MQ=255)
                                                   gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCt                                      <  1:1784274/54‑1 (MQ=255)
                                                    tttCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACTCTATCTCTCGCTGATTGGCg                    <  1:4012040/71‑1 (MQ=255)
                                                     ttCTCGCCAATTCCGGCAATGTCGATGGTTTCTTa                                                       >  1:577285/1‑35 (MQ=255)
                                                      tctcGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCgg                   >  1:1282069/1‑70 (MQ=255)
                                                          gCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTgct                         >  1:1212902/1‑58 (MQ=255)
                                                          gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                          >  1:590717/1‑59 (MQ=255)
                                                          gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                          >  1:570453/1‑59 (MQ=255)
                                                          gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTa                 <  1:3483918/68‑1 (MQ=255)
                                                          gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACt               <  1:471673/70‑1 (MQ=255)
                                                              aTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt          <  1:3176962/71‑1 (MQ=255)
                                                               ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTg                                            <  1:4108291/36‑1 (MQ=255)
                                                               ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                            <  1:1805796/52‑1 (MQ=255)
                                                                tCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCgg                   >  1:3544688/1‑60 (MQ=255)
                                                                 ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa        <  1:1000774/70‑1 (MQ=255)
                                                                 ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa        <  1:1498431/70‑1 (MQ=255)
                                                                 ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa        <  1:4095222/70‑1 (MQ=255)
                                                                 ccGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa        <  1:407164/70‑1 (MQ=255)
                                                                   ggCAATGTCGATGGTTTCTTACGCATCGGGTTCACGCt                                      <  1:3525044/38‑1 (MQ=255)
                                                                     cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTa                 >  1:3461843/1‑57 (MQ=255)
                                                                     cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc   <  1:3395428/71‑1 (MQ=255)
                                                                      aaTGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCt  <  1:1805607/71‑1 (MQ=255)
                                                                      |                                                                      
CCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCT  >  minE/832054‑832194

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