Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,118 T→C 7.0% I192I (ATT→ATC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,1180TC 7.0% 160.4 / 3.9 57I192I (ATT→ATCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (30/23);  new base C (2/2);  total (32/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

ATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTT  >  minE/832050‑832186
                                                                    |                                                                    
aTCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTcc                                                                    >  1:3972932/1‑71 (MQ=255)
   gccgcATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGc                                                                 <  1:3349131/71‑1 (MQ=255)
    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)
                    cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCgg                                                                  >  1:636605/1‑53 (MQ=255)
                    cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATg                                                             >  1:2676481/1‑58 (MQ=255)
                    cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCt                                                 <  1:640374/70‑1 (MQ=255)
                    cacCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCt                                                 <  1:1578860/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)
                            aaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGc                                                                 <  1:3257603/46‑1 (MQ=255)
                            aaaaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCg                                        <  1:2911752/71‑1 (MQ=255)
                              aaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCg                                                                   >  1:4105704/1‑42 (MQ=255)
                              aaCTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCg                                                                   >  1:689351/1‑42 (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:3746947/1‑70 (MQ=255)
                                                 ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                    >  1:459835/1‑70 (MQ=255)
                                                 ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                    <  1:1586102/70‑1 (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:570453/1‑59 (MQ=255)
                                                              gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                  >  1:590717/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)
                                                                    |                                                                    
ATCGCCGCATCTAACGTTTACACCATCAAAAACTACGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTT  >  minE/832050‑832186

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: