Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,277 C→T 24.7% Y245Y (TAC→TAT narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,2770CT24.7% 145.0 / 46.9 81Y245Y (TAC→TATnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (23/38);  new base T (11/9);  total (34/47)
Fisher's exact test for biased strand distribution p-value = 2.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.51e-01

CGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTT  >  minE/832208‑832346
                                                                     |                                                                     
cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                                                                       <  1:3645032/70‑1 (MQ=255)
cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                                                                       <  1:1217644/70‑1 (MQ=255)
    ttGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc                                                                  >  1:2335376/1‑71 (MQ=255)
       cctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa                                                                >  1:1639362/1‑70 (MQ=255)
         tcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa                                                                <  1:2042996/68‑1 (MQ=255)
         tcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCa                                                                <  1:2338139/68‑1 (MQ=255)
           ctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAc                                                           <  1:2468480/71‑1 (MQ=255)
             gCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat                                                         >  1:2493073/1‑71 (MQ=255)
              cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcat                                                         >  1:2353916/1‑70 (MQ=255)
              cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAca                                                          <  1:3653117/69‑1 (MQ=255)
               gTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatca                                                       <  1:784003/71‑1 (MQ=255)
               gTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatc                                                        >  1:3764119/1‑70 (MQ=255)
                   ccGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCt                                                             >  1:2402554/1‑61 (MQ=255)
                   ccGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACatcatc                                                     <  1:2893811/69‑1 (MQ=255)
                   ccGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGc                                                   >  1:3814001/1‑71 (MQ=255)
                     gCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGc                                                   >  1:907856/1‑69 (MQ=255)
                     gCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCt                                                 <  1:700540/71‑1 (MQ=255)
                     gCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCt                                                 <  1:844787/71‑1 (MQ=255)
                      cAAACCTGGGGCGAGCAAACTGACGTACCGTAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg                                                <  1:2627818/71‑1 (MQ=255)
                      cAAACCTGGGGCGAGCAAACTGACGTACCGTAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg                                                <  1:1783589/71‑1 (MQ=255)
                      cAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg                                                >  1:457520/1‑71 (MQ=255)
                      atgaCCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTg                                                <  1:634745/68‑1 (MQ=255)
                       aaaCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTgg                                               >  1:224929/1‑71 (MQ=255)
                       aaaCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTgg                                               <  1:2405470/71‑1 (MQ=255)
                       aaaCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg                                                <  1:1740443/70‑1 (MQ=255)
                        gaCCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTggg                                              >  1:2159452/2‑71 (MQ=255)
                         aCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACaa                                                                     <  1:403347/47‑1 (MQ=255)
                         aCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTgggg                                             >  1:1181728/1‑71 (MQ=255)
                         aCCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTgggg                                             <  1:216741/71‑1 (MQ=255)
                           cTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTgggg                                             >  1:3310098/1‑69 (MQ=255)
                            tGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGt                                           >  1:2265544/1‑70 (MQ=255)
                              gggCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCa                                         <  1:3053453/70‑1 (MQ=255)
                              gggCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCa                                         <  1:3753141/70‑1 (MQ=255)
                                 cGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACg                                     <  1:259794/71‑1 (MQ=255)
                                  gAGCAAACCGACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACg                                     >  1:98572/1‑70 (MQ=255)
                                    gctaaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGc                                  <  1:1355083/68‑1 (MQ=255)
                                    gCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGc                                  <  1:2662689/71‑1 (MQ=255)
                                     cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                                 <  1:618696/71‑1 (MQ=255)
                                     cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                                 <  1:1191578/71‑1 (MQ=255)
                                     cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                                 >  1:1545580/1‑71 (MQ=255)
                                     cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGc                                  >  1:1217330/1‑70 (MQ=255)
                                      aaaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCaa                                        <  1:2915844/63‑1 (MQ=255)
                                       aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGc                               <  1:3461661/71‑1 (MQ=255)
                                       aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCACATCATCGCCTGGGGGTCAAAc                                      >  1:2568647/1‑63 (MQ=255)
                                       aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCACATCATCGCCTGGGGGTCAAAc                                      >  1:2409408/1‑63 (MQ=255)
                                        accgACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCg                                >  1:2536065/4‑69 (MQ=255)
                                        aCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCg                                >  1:2502871/1‑69 (MQ=255)
                                            aCGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGa                            <  1:3652509/69‑1 (MQ=255)
                                            aCGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAc                           <  1:3415441/70‑1 (MQ=255)
                                              gTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTgggg                                             <  1:120462/50‑1 (MQ=38)
                                              gTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAcacg                        <  1:2666827/71‑1 (MQ=255)
                                               tACCGGAATCTGCTGACTGGTACAACTCCAGCTACatcatc                                                     <  1:3269632/41‑1 (MQ=255)
                                               tACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGt                       <  1:225416/71‑1 (MQ=255)
                                               tACCGGAATCTGCCGACTGGTATAACTCCAGCTACTTCATCGCCTGGGGGTCTAACGTACCGCAGAcac                         <  1:228291/69‑3 (MQ=255)
                                               tACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAcacg                        <  1:1049604/70‑1 (MQ=255)
                                                aCCGGAATCTGCCGACTGGTATAACTCCAGCTAGATCATCGCCTGGGGGTCTAACGTACCGCAGACACGt                       <  1:1323966/70‑1 (MQ=255)
                                                aCCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTa                      >  1:2342308/1‑71 (MQ=255)
                                                aCCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTa                      >  1:2193418/1‑71 (MQ=255)
                                                 ccGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGa                            <  1:3919279/64‑1 (MQ=255)
                                                 ccGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAcacgaa                      >  1:1522485/1‑68 (MQ=255)
                                                  cGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGt                                           >  1:161398/1‑48 (MQ=255)
                                                  cGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAcc                    <  1:2679312/71‑1 (MQ=255)
                                                  cGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGAca                          >  1:1701875/1‑64 (MQ=255)
                                                   ggAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAccc                   >  1:3743391/1‑71 (MQ=255)
                                                    gAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTAcgcc                  <  1:390411/71‑1 (MQ=255)
                                                     aaTCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGTTCAAACGTGCCGCAGACGCGTACCCCg                 <  1:485721/71‑1 (MQ=255)
                                                     aaTCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCg                 <  1:3692813/71‑1 (MQ=255)
                                                     aaTCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGa                            >  1:2412859/1‑60 (MQ=255)
                                                      aactgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCa                              <  1:211184/55‑1 (MQ=255)
                                                       tCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTct                                         <  1:188414/45‑2 (MQ=38)
                                                         tgccgACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGacg             >  1:170482/1‑69 (MQ=255)
                                                          gccgACTGGTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGACACGTACGCCGGACGc            >  1:1510936/4‑71 (MQ=255)
                                                            tgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATg             <  1:887414/68‑1 (MQ=255)
                                                               cTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCAc        <  1:971018/70‑1 (MQ=255)
                                                                tGGTACAACTCCAGCTACATCATCGCCTGGGGGTCa                                         <  1:2318042/36‑1 (MQ=255)
                                                                tGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActt      >  1:2942217/1‑71 (MQ=255)
                                                                 ggTACAACTCCAGCTACATCATCGCCTGGGGGTCa                                         <  1:2889815/35‑1 (MQ=255)
                                                                  gTATAACTCCAGCTACATCATCGCCTGGGGGTCTAACGTACCGCAGCCACGTACGCCGGACGCCCActtct    >  1:2614210/1‑71 (MQ=255)
                                                                  gTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtct    >  1:2291785/1‑71 (MQ=255)
                                                                  gTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACACGTACCCCGGATGCTCActtct    <  1:1890624/71‑1 (MQ=255)
                                                                  gTACAACTCCAGCTAAATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtc     <  1:3567884/70‑1 (MQ=255)
                                                                   tACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCActtctt   <  1:1598700/71‑1 (MQ=255)
                                                                    aCAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCAGGATGATCActtctt   >  1:3898303/1‑70 (MQ=255)
                                                                     cAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCttt  <  1:599472/70‑1 (MQ=255)
                                                                     |                                                                     
CGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCTCACTTCTTT  >  minE/832208‑832346

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: