Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,154 A→G 10.7% A204A (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,1540AG10.7% 246.7 / 17.9 93A204A (GCA→GCG) ‡narGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (46/37);  new base G (4/6);  total (50/43)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

CGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCG  >  minE/832085‑832223
                                                                     |                                                                     
cGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCa                                                                       <  1:1010444/70‑1 (MQ=255)
    ccGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTa                                                                  >  1:436276/1‑71 (MQ=255)
      ggACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtc                                                                >  1:616363/1‑71 (MQ=255)
      ggACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtc                                                                >  1:3637623/1‑71 (MQ=255)
       gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct                                                               <  1:1066680/71‑1 (MQ=255)
       gACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtct                                                               <  1:1695461/71‑1 (MQ=255)
        aCCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctc                                                              <  1:14805/71‑1 (MQ=255)
          cGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAtctctc                                                            <  1:1882145/71‑1 (MQ=255)
             gtTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAt                                                                 >  1:2504553/1‑63 (MQ=255)
              ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                         >  1:1338445/1‑70 (MQ=255)
              ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                         >  1:1370585/1‑70 (MQ=255)
              ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCt                                                         >  1:3640714/1‑70 (MQ=255)
              ttGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCAAGCTATCTCTCGCt                                                         >  1:3011493/1‑70 (MQ=255)
               tGCCGGTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTgc                                                       >  1:1516964/1‑70 (MQ=255)
                gCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt                                                      >  1:3524794/1‑71 (MQ=255)
                 ccggTTTCTCGCCGATCCCGGCTATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCtt                                                     <  1:1692044/69‑1 (MQ=255)
                 cTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAtt                                                     <  1:276185/71‑1 (MQ=255)
                 cTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt                                                      <  1:3205353/70‑1 (MQ=255)
                  cggTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTg                                                    <  1:1763539/70‑1 (MQ=255)
                    gTTTCTCGCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGc                                                  <  1:1921210/71‑1 (MQ=255)
                    gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                       >  1:688888/1‑66 (MQ=255)
                    gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                       >  1:1134028/1‑66 (MQ=255)
                    gTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGc                                                  <  1:3392997/71‑1 (MQ=255)
                       tctcGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGATATCTCTCGCTGATTGGCgg                                                >  1:503240/1‑70 (MQ=255)
                        ctcGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTa                                              >  1:34490/1‑71 (MQ=255)
                           gCCGATCCCGGCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCAcct                                           <  1:2778830/71‑3 (MQ=255)
                           gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTa                                                                  >  1:3001070/1‑48 (MQ=255)
                           gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                       >  1:754515/1‑59 (MQ=255)
                           gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTAc                                             >  1:3864520/1‑69 (MQ=255)
                           gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTAc                                             >  1:1179512/1‑69 (MQ=255)
                           gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt                                           >  1:2991987/1‑71 (MQ=255)
                           gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt                                           <  1:604534/71‑1 (MQ=255)
                           gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt                                           <  1:741763/71‑1 (MQ=255)
                           gCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt                                           <  1:3831561/71‑1 (MQ=255)
                            ccAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTg                                          >  1:2483212/1‑71 (MQ=255)
                             cAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGc                                         <  1:3323066/71‑1 (MQ=255)
                             cAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGc                                         >  1:3388541/1‑71 (MQ=255)
                              aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTgg                                                   >  1:3641737/1‑60 (MQ=255)
                              aaTTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCt                                        >  1:1457850/1‑71 (MQ=255)
                                ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCt                                        >  1:3323626/1‑69 (MQ=255)
                                ttCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGAGGTACTTGCtt                                       >  1:3808029/1‑70 (MQ=255)
                                 tCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAt                                                                 >  1:1182559/1‑43 (MQ=255)
                                 tCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCtt                                       >  1:2817410/1‑69 (MQ=255)
                                   cGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTAc                                             <  1:3705275/61‑1 (MQ=255)
                                   cGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAg                                    <  1:2844012/70‑1 (MQ=255)
                                    ggCGATGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAgtt                                  <  1:2549326/71‑3 (MQ=255)
                                    ggCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTAt                                                                 <  1:878847/40‑1 (MQ=255)
                                     gcgaTGTCGATGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGttt                                 >  1:3752547/4‑71 (MQ=255)
                                     gCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt                                           <  1:2404714/61‑1 (MQ=255)
                                     gCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACtt                                           <  1:3367791/61‑1 (MQ=255)
                                     gCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCtt                                 >  1:2286859/1‑71 (MQ=255)
                                      cAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCtt                                 >  1:3806575/1‑70 (MQ=255)
                                       aaTGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTgg                                         >  1:816553/1‑60 (MQ=255)
                                         tGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGa                                                       <  1:2568209/45‑1 (MQ=255)
                                          gTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTaa                                     <  1:3395008/62‑1 (MQ=255)
                                           tCGATGGTTTCTTACGCATCGGGTGCACGCTAtctc                                                              <  1:2996612/36‑1 (MQ=255)
                                           tCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                           <  1:1302671/71‑1 (MQ=255)
                                              aTGGTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTAAGTTTCTACGACTg                        <  1:1520559/71‑1 (MQ=255)
                                              aTGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGAc                          <  1:2846263/69‑1 (MQ=255)
                                               tGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt                                                      >  1:2253290/1‑40 (MQ=255)
                                               tGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGAt                                                      >  1:2485685/1‑40 (MQ=255)
                                                ggTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTaa                                     >  1:3128065/1‑56 (MQ=255)
                                                ggTTTCTTACGCCGCCGGAACGCGTTATCTGTCGCTGCTTGGCGGCACCTGTTTaa                                     >  1:83013/1‑56 (MQ=255)
                                                ggTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGc                                         >  1:1223452/1‑52 (MQ=255)
                                                    tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACt                                            >  1:1491344/1‑45 (MQ=255)
                                                    tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACt                   <  1:1480589/70‑1 (MQ=255)
                                                    tCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTg                  >  1:2595552/1‑71 (MQ=255)
                                                      ttACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTc                                <  1:585165/55‑1 (MQ=255)
                                                        aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGc                 <  1:1296926/68‑1 (MQ=255)
                                                        aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGa               <  1:1987235/70‑1 (MQ=255)
                                                        aCGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGa               <  1:2233866/70‑1 (MQ=255)
                                                         cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc              <  1:3524280/70‑1 (MQ=255)
                                                         cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt             <  1:1830429/71‑1 (MQ=255)
                                                         cGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt             <  1:2150949/71‑1 (MQ=255)
                                                          gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt            >  1:1258687/1‑71 (MQ=255)
                                                          gCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt            >  1:917353/1‑71 (MQ=255)
                                                          gCACCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTAc                             >  1:3469226/1‑54 (MQ=255)
                                                           cATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTa                     <  1:1785476/61‑1 (MQ=255)
                                                            aTCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGc          >  1:396847/1‑71 (MQ=255)
                                                              cGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcc         >  1:246326/1‑70 (MQ=255)
                                                               gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc       >  1:2472942/1‑71 (MQ=255)
                                                               gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc       <  1:2526776/71‑1 (MQ=255)
                                                               gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcctc       <  1:3345869/71‑1 (MQ=255)
                                                               gggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGcccc       <  1:1947401/71‑3 (MQ=255)
                                                                ggTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACCTGcctcc      >  1:2531308/1‑71 (MQ=255)
                                                                 gTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGTTTCTACGACTGGTACTGCGACTTGcctcct     <  1:2330743/71‑1 (MQ=255)
                                                                  tGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGAc              <  1:3529190/61‑1 (MQ=255)
                                                                  tGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACt             <  1:3922766/62‑1 (MQ=255)
                                                                    cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGa                           >  1:3525635/1‑46 (MQ=255)
                                                                    cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACtt            >  1:1317014/1‑61 (MQ=255)
                                                                    cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg  >  1:2843832/1‑71 (MQ=255)
                                                                    cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg  >  1:1408826/1‑71 (MQ=255)
                                                                    cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg  >  1:2483430/1‑71 (MQ=255)
                                                                    cACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCg  <  1:2010019/71‑1 (MQ=255)
                                                                     |                                                                     
CGGCCCGGACCGTGTTGCTGGTTTCTCGCCAATTCCGGCAATGTCGATGGTTTCTTACGCATCGGGTGCACGCTATCTCTCGCTGATTGGCGGTACTTGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCG  >  minE/832085‑832223

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: