Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,268 T→C 9.8% A242A (GCT→GCC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,2680TC 9.8% 138.7 / 6.9 51A242A (GCT→GCCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (2/3);  total (25/26)
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.62e-01

GTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCT  >  minE/832202‑832337
                                                                  |                                                                     
gTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGAc                                                                    <  1:714021/70‑1 (MQ=255)
gTATTGCGACCTGCCGCCCGCCTCGCCCATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGAc                                                                    <  1:1639320/70‑1 (MQ=255)
gTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGAc                                                                    <  1:3989238/70‑1 (MQ=255)
      cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctgct                                                                       <  1:3382371/61‑1 (MQ=255)
      cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                                                              <  1:892855/70‑1 (MQ=255)
       gACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACaa                                                            <  1:2306942/71‑1 (MQ=255)
       gACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACaa                                                            <  1:1342047/71‑1 (MQ=255)
           tGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGAc                                                                    <  1:3017732/59‑1 (MQ=255)
           tGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc                                                         >  1:1647249/1‑70 (MQ=255)
            gCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc                                                         >  1:140720/1‑69 (MQ=255)
              ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGc                                                     >  1:1382133/1‑71 (MQ=255)
                  tGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAca                                                 >  1:3500645/1‑71 (MQ=255)
                    cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAg                                                      >  1:2639592/1‑64 (MQ=255)
                    cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatc                                               <  1:1309794/71‑1 (MQ=255)
                      tctcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc                                                         <  1:2407774/59‑1 (MQ=255)
                        tcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGa                                                                     >  1:2961882/1‑45 (MQ=255)
                        tcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTAc                                                              <  1:2716643/52‑1 (MQ=255)
                          cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatcat                                             >  1:2371338/1‑67 (MQ=255)
                          cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatcat                                             >  1:2258489/1‑67 (MQ=255)
                          cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCg                                           <  1:993401/69‑1 (MQ=255)
                          cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCg                                           <  1:3768986/69‑1 (MQ=255)
                          cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGcc                                         >  1:360151/1‑71 (MQ=255)
                            cAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTg                                       <  1:3718717/71‑1 (MQ=255)
                            atgaCCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATCGCCTg                                       >  1:3698067/4‑71 (MQ=255)
                              aaCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCg                                           <  1:1831548/65‑1 (MQ=255)
                              aaCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCg                                           <  1:1072517/65‑1 (MQ=255)
                               aCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTgggg                                    <  1:428872/71‑1 (MQ=255)
                                 cTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTaaa                                                             >  1:3352082/1‑42 (MQ=255)
                                  tGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGt                                  <  1:2706459/70‑1 (MQ=255)
                                    gggCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCa                                <  1:3266878/70‑1 (MQ=255)
                                       cGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACg                            <  1:2623916/71‑1 (MQ=255)
                                        gAGCAAACCGACGTACCGGAATCTGCCGACTGGTATAACTCCAGCTACATCATAGcc                                         <  1:4011224/57‑1 (MQ=38)
                                           cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCaaa                              <  1:2199315/65‑1 (MQ=255)
                                           cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTg                          <  1:135780/69‑1 (MQ=255)
                                           cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGcc                        >  1:1066559/1‑71 (MQ=255)
                                           cAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCACCTACATCATCGCCTgggg                                    >  1:1687978/1‑59 (MQ=255)
                                             aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGc                                          >  1:1672224/1‑51 (MQ=255)
                                             aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGc                                          >  1:4071052/1‑51 (MQ=255)
                                             aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGc                      <  1:575698/71‑1 (MQ=255)
                                               cTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCg                       >  1:2777381/1‑68 (MQ=255)
                                               cTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCg                       >  1:2705013/1‑68 (MQ=255)
                                                tGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGc                                          >  1:3235057/1‑48 (MQ=255)
                                                tGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCa                     <  1:3609393/69‑1 (MQ=255)
                                                   cGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGAcgc                <  1:1765693/71‑1 (MQ=255)
                                                       ccGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCa                     >  1:2883012/1‑62 (MQ=255)
                                                       ccGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGa                   >  1:3465248/1‑64 (MQ=255)
                                                       ccGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAc            >  1:3120513/1‑71 (MQ=255)
                                                       ccGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAc            >  1:1066170/1‑71 (MQ=255)
                                                       ccGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAc            >  1:482385/1‑71 (MQ=255)
                                                       ccGGAATCTGCCGACTGGTATAACTCCAGCTAcat                                                >  1:1107212/1‑35 (MQ=25)
                                                        cGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTAcc           <  1:3302991/71‑1 (MQ=255)
                                                         ggAATCTGCTGACTGGTACAACTCCAGCTAcatcat                                             >  1:1775448/1‑36 (MQ=255)
                                                               tgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATg    >  1:2155553/1‑71 (MQ=255)
                                                               tgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATg    >  1:2524471/1‑71 (MQ=255)
                                                               tgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATg    >  1:4142192/1‑71 (MQ=255)
                                                               tgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATg    >  1:2769933/1‑71 (MQ=255)
                                                               tgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATg    >  1:3351299/1‑71 (MQ=255)
                                                               tgctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATg    >  1:2936683/1‑71 (MQ=255)
                                                                gctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGc   >  1:3603283/1‑71 (MQ=255)
                                                                 ctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCt  <  1:223210/71‑1 (MQ=255)
                                                                 ctgACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCt  >  1:2700456/1‑71 (MQ=255)
                                                                  |                                                                     
GTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCAGACGCGTACCCCGGATGCT  >  minE/832202‑832337

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: