Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 832,250 T→C 7.1% T236T (ACT→ACC narG → nitrate reductase 1, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE832,2500TC 7.1% 200.0 / 6.2 70T236T (ACT→ACCnarGnitrate reductase 1, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (29/36);  new base C (2/3);  total (31/39)
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.24e-01

TGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCA  >  minE/832182‑832318
                                                                    |                                                                    
tGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGa                                                                    >  1:3774770/1‑71 (MQ=255)
 gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTg                                                                     <  1:3651997/69‑1 (MQ=255)
 gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGa                                                                    <  1:703097/70‑1 (MQ=255)
 gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGa                                                                    <  1:1584377/70‑1 (MQ=255)
 gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGAc                                                                   >  1:4146358/1‑71 (MQ=255)
 gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGAc                                                                   >  1:2519174/1‑71 (MQ=255)
 gCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGAc                                                                   >  1:2070949/1‑71 (MQ=255)
  cTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                                                                  <  1:3013797/71‑1 (MQ=255)
    tAAGTTTCTACGACTGGTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTa                                                                >  1:3256040/1‑71 (MQ=255)
    tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACt                                                                      >  1:4072358/1‑65 (MQ=255)
    tAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACt                                                                      >  1:3910294/1‑65 (MQ=255)
     aaGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTAc                                                               >  1:505476/1‑71 (MQ=255)
      aGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTAcc                                                              <  1:1728745/71‑1 (MQ=255)
          tCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                                           <  1:3987790/70‑1 (MQ=255)
           cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                                           <  1:2528593/69‑1 (MQ=255)
           cTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGa                                                           <  1:3622156/69‑1 (MQ=255)
              cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACg                                                                  >  1:1656659/1‑59 (MQ=255)
              cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctg                                                      <  1:2087521/71‑1 (MQ=255)
              cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctg                                                      <  1:3218110/71‑1 (MQ=255)
              cGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctg                                                      <  1:2378867/71‑1 (MQ=255)
                    gTATTGCGACCTGCCGCCCGCCTCGCCGATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGAc                                                 <  1:714021/70‑1 (MQ=255)
                    gTATTGCGACCTGCCGCCCGCCTCGCCCATGACCTGGGGCGAGCAAACCGACGTACCGGAATCTGCCGAc                                                 <  1:1639320/70‑1 (MQ=255)
                    gTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGAc                                                 <  1:3989238/70‑1 (MQ=255)
                        tGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGa                                                                    <  1:1191767/47‑1 (MQ=255)
                         gCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                                         >  1:2831502/1‑57 (MQ=255)
                          cGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTg                                                                     >  1:1475862/1‑44 (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)
                                  ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctg                                                      >  1:3121036/1‑51 (MQ=255)
                                  ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATctg                                                      >  1:404510/1‑51 (MQ=255)
                                  ctcctGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGc                                  >  1:1382133/1‑71 (MQ=255)
                                      tGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAca                              >  1:3500645/1‑71 (MQ=255)
                                       gCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                                         >  1:1109659/1‑43 (MQ=255)
                                        cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAg                                   >  1:2639592/1‑64 (MQ=255)
                                        cGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTAcatc                            <  1:1309794/71‑1 (MQ=255)
                                          tctcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTc                                      <  1:2407774/59‑1 (MQ=255)
                                            tcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                                         >  1:1001760/1‑38 (MQ=255)
                                            tcCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAAt                                                         >  1:1706253/1‑38 (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:3768986/69‑1 (MQ=255)
                                              cGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCg                        <  1:993401/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)
                                                 aaaCCTGGGGCGAGCAAACTGACGTACCGGAATctg                                                      <  1:650133/36‑1 (MQ=255)
                                                 aaaCCTGGGGCGAGCAAACTGACGTACCGGAATctg                                                      <  1:1644265/36‑1 (MQ=255)
                                                 aaaCCTGGGGCGAGCAAACTGACGTACCGGAATctg                                                      <  1:1122517/36‑1 (MQ=255)
                                                  aaCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCg                        <  1:1072517/65‑1 (MQ=255)
                                                  aaCCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCg                        <  1:1831548/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:4071052/1‑51 (MQ=255)
                                                                 aaCTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGc                       >  1:1672224/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)
                                                                    |                                                                    
TGCTTAAGCTTCTACGACTGGTACTGCGACTTGCCTCCTGCGTCTCCGCAAACCTGGGGCGAGCAAACTGACGTACCGGAATCTGCTGACTGGTACAACTCCAGCTACATCATCGCCTGGGGGTCAAACGTGCCGCA  >  minE/832182‑832318

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: