Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,405 C→T 18.3% F41F (TTC→TTT narH → nitrate reductase 1, beta (Fe‑S) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE835,4050CT18.3% 132.8 / 22.2 60F41F (TTC→TTTnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base C (23/26);  new base T (5/6);  total (28/32)
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.71e-01

GCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACA  >  minE/835335‑835469
                                                                      |                                                                
gCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTc                                                                  <  1:2074764/71‑1 (MQ=255)
  tGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaa                                                                >  1:2051205/1‑71 (MQ=255)
  tGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaa                                                                >  1:2398714/1‑71 (MQ=255)
  tGCCACACCTGCTCGATCACCTGCAAGAACGTCTGGACCAGCCGCGAAGGCATGGAGTACGCCTGGTTcaa                                                                >  1:709105/1‑71 (MQ=255)
    ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                               >  1:2892066/1‑70 (MQ=255)
    ccATACCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTtaac                                                               >  1:1854520/1‑70 (MQ=255)
    ccACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaa                                                                <  1:3337043/69‑1 (MQ=255)
       tacCTGTTCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACg                                                           <  1:2538245/70‑1 (MQ=255)
       cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                                           <  1:2316976/71‑1 (MQ=255)
       cacCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACg                                                           <  1:2439287/71‑1 (MQ=255)
        acCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGt                                                          >  1:2166895/1‑71 (MQ=255)
          cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaaca                                                              >  1:1893268/1‑65 (MQ=255)
          cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTgg                                                        <  1:1327911/71‑1 (MQ=255)
          cTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTg                                                         >  1:3401816/1‑70 (MQ=255)
             ttgaGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGaa                                                      <  1:2134511/67‑1 (MQ=255)
             ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaa                                                     >  1:2118523/1‑71 (MQ=255)
             ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaa                                                      <  1:243508/70‑1 (MQ=255)
             ttCGGTGACCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGa                                                       <  1:2919966/69‑1 (MQ=255)
             ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                                                 <  1:3307356/59‑1 (MQ=255)
             ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                       <  1:3215110/69‑1 (MQ=255)
              tCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                                               >  1:699739/1‑60 (MQ=255)
              tCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                                               >  1:1884875/1‑60 (MQ=255)
               cAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAc                                                    <  1:3762977/70‑1 (MQ=255)
                aGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCa                                                  <  1:1729499/71‑1 (MQ=255)
                 gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAcc                                                   >  1:2873644/1‑69 (MQ=255)
                 gTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                                                 <  1:3256987/71‑1 (MQ=255)
                  ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTca                                                                 <  1:1815241/54‑1 (MQ=255)
                  ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAcc                                                   <  1:2143466/68‑1 (MQ=255)
                  ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                                <  1:942218/71‑1 (MQ=255)
                  tgaCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACg                                                           <  1:3412602/58‑1 (MQ=255)
                    aCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                                         >  1:2752404/1‑59 (MQ=255)
                    aCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaaccaaacc                                              >  1:690504/1‑71 (MQ=255)
                    aCCTGTAAAAACGTCTGGACCGGGCGAGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                                         >  1:3463912/1‑59 (MQ=255)
                    aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                                                 <  1:1861070/68‑1 (MQ=255)
                    aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGcc                                              >  1:3361058/1‑71 (MQ=255)
                     ccTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGcc                                              <  1:2274587/70‑1 (MQ=255)
                      cTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGtc                                                         <  1:2574569/58‑2 (MQ=255)
                        ttAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                                           >  1:3332404/2‑70 (MQ=255)
                         tAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTaacaac                                                            <  1:1343900/52‑1 (MQ=255)
                         tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaacaa                                                             >  1:2718493/1‑51 (MQ=255)
                         tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGcc                                         <  1:3938006/71‑1 (MQ=255)
                         tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGcc                                         <  1:188282/71‑1 (MQ=255)
                         tAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                                          >  1:1180552/1‑70 (MQ=255)
                           aaaaCGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGAAACCAAACCGGGcatt                                       >  1:1561972/1‑68 (MQ=255)
                           aaaaCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCa                                        <  1:3446346/70‑1 (MQ=255)
                              aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                                               >  1:1052481/1‑44 (MQ=255)
                              aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTcaac                                                               >  1:1236055/1‑44 (MQ=255)
                              aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGc                                    >  1:2327977/1‑71 (MQ=255)
                              aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGc                                    >  1:1035147/1‑71 (MQ=255)
                              aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAAGGc                                    <  1:1573352/71‑1 (MQ=255)
                               cGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCt                                   <  1:1203280/71‑1 (MQ=255)
                                  cTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCa                                        >  1:2080787/1‑63 (MQ=255)
                                     gACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                               <  1:2819519/69‑1 (MQ=255)
                                          gCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                        >  1:3545396/1‑71 (MQ=255)
                                           ccGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCa                                                  <  1:2673428/44‑1 (MQ=255)
                                           ccGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgact                       <  1:1795505/71‑1 (MQ=255)
                                           ccGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgact                       >  1:3544092/1‑71 (MQ=255)
                                            cGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGactgactg                      >  1:546169/1‑71 (MQ=255)
                                            cGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGactgactg                      >  1:3840697/1‑71 (MQ=255)
                                              tGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCtt                                  >  1:3829742/1‑57 (MQ=255)
                                               gAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                                           >  1:299584/1‑47 (MQ=255)
                                               gAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                                           >  1:295014/1‑47 (MQ=255)
                                               gAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTccc                               >  1:2759982/1‑59 (MQ=255)
                                                   gCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAAc               <  1:1808116/71‑1 (MQ=255)
                                                       ggAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGc                                          >  1:1609513/1‑40 (MQ=255)
                                                        gAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGa                   <  1:979206/62‑1 (MQ=255)
                                                          atatGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaa         <  1:3411636/70‑1 (MQ=255)
                                                             tGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGaaaaa      >  1:729834/1‑70 (MQ=255)
                                                             tGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAAt     <  1:2451514/71‑1 (MQ=255)
                                                              gCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAg            >  1:1110957/1‑63 (MQ=255)
                                                               cGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATa    >  1:733953/1‑70 (MQ=255)
                                                                gTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATAc   >  1:3680573/1‑70 (MQ=255)
                                                                  ggTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACa  >  1:3216736/1‑69 (MQ=255)
                                                                   gTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAgg           <  1:2884836/59‑1 (MQ=255)
                                                                    ttCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGactgactg                      <  1:2581109/47‑1 (MQ=255)
                                                                      |                                                                
GCTGCCACACCTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACA  >  minE/835335‑835469

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: