Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 835,414 G→C 9.6% V44V (GTG→GTC 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,4140GC 9.6% 142.4 / 5.8 52V44V (GTG→GTCnarHnitrate reductase 1, beta (Fe‑S) subunit
Reads supporting (aligned to +/- strand):  ref base G (24/23);  new base C (3/2);  total (27/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.84e-01

CTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTG  >  minE/835345‑835479
                                                                     |                                                                 
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)
   ttCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGa                                                                 <  1:3215110/69‑1 (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)
        ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAAcc                                                             <  1:2143466/68‑1 (MQ=255)
        ttACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAg                                                          <  1:942218/71‑1 (MQ=255)
          aCCTGTAAAAACGTCTGGACCGGGCGCGAAGGCATGGAGTACGCATGGTTTAACAACGTCGaaaccaaacc                                                        >  1:690504/1‑71 (MQ=255)
          aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCaa                                                           <  1:1861070/68‑1 (MQ=255)
          aCCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGcc                                                        >  1:3361058/1‑71 (MQ=255)
           ccTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGcc                                                        <  1:2274587/70‑1 (MQ=255)
              ttAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCggg                                                     >  1:3332404/2‑70 (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)
                    aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGc                                              >  1:1035147/1‑71 (MQ=255)
                    aCGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGc                                              >  1:2327977/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:3544092/1‑71 (MQ=255)
                                 ccGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgact                                 <  1:1795505/71‑1 (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)
                                                              acaacGTGGAAACCAAGCCGGGCCAGGGCTTCCCgactgac                                  >  1:1671379/1‑41 (MQ=255)
                                                              acaacGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAggcggc    >  1:1564678/1‑71 (MQ=255)
                                                                aacGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTg  >  1:3576067/1‑71 (MQ=255)
                                                                     |                                                                 
CTGTTCAGTTACCTGTAAAAACGTCTGGACCAGCCGTGAAGGCGTGGAATATGCGTGGTTCAACAACGTGGAAACCAAGCCGGGCCAGGGCTTCCCGACTGACTGGGAAAACCAGGAAAAATACAAAGGCGGCTG  >  minE/835345‑835479

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: