Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,822,088 C→T 7.4% S465S (TCC→TCT norV → flavorubredoxin oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,822,0880CT 7.4% 204.8 / 6.1 67S465S (TCC→TCTnorVflavorubredoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (21/41);  new base T (3/2);  total (24/43)
Fisher's exact test for biased strand distribution p-value = 3.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

GCAGGACGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGATCATCGGTTC  >  minE/1822022‑1822158
                                                                  |                                                                      
gCAGGACGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCg                                                                    >  1:1037953/1‑71 (MQ=255)
gCAGGACGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCg                                                                    >  1:986216/1‑71 (MQ=255)
gCAGGACGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCg                                                                    >  1:79652/1‑71 (MQ=255)
   ggACGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCa                                                                 <  1:1577599/71‑1 (MQ=255)
   ggACGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCa                                                                 <  1:15510/71‑1 (MQ=255)
   ggACGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCa                                                                 <  1:1378293/71‑1 (MQ=255)
     aCGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCTGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCaaa                                                               <  1:1846075/71‑1 (MQ=255)
     aCGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCaaa                                                               <  1:3443201/71‑1 (MQ=255)
     aCGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCaaa                                                               <  1:3507663/71‑1 (MQ=255)
     aCGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCaaa                                                               <  1:593504/71‑1 (MQ=255)
     aCGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCaaa                                                               <  1:3094061/71‑1 (MQ=255)
         tGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACg                                                           >  1:2045001/1‑71 (MQ=255)
         tGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTACTCTGCCCGGAATGCTCCCTCGGCAAAGACg                                                           >  1:915029/1‑71 (MQ=255)
              cAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCttt                                                      >  1:3175304/1‑71 (MQ=255)
                     gCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAAc                                               <  1:204297/71‑1 (MQ=255)
                       cGTGGAGTGAAGTCCCGGATAACTTCCTGTGCCCGGAATGCTCTCTCGGCAAAGACGTCttt                                                      <  1:3934898/62‑1 (MQ=255)
                         tgtaGTGAAGTCCCGGATAACTTCCTGTGCCCGGAATGCTCTCTCGGCAAAGACGTCTTTga                                                    >  1:4098216/4‑62 (MQ=255)
                         tGGAGTGAAGTCCCGGATAACTTCCTGTGCCCGGAATGCTCTCTCGGCAAAGACGTCTTTga                                                    >  1:1682524/1‑62 (MQ=255)
                         tGGAGTGAAGTCCCGGATAACTTCCTGTGCCCGGAATGCTCTCTCGGCAAAGACGTCTTTga                                                    >  1:401771/1‑62 (MQ=255)
                              tGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCgg                                      <  1:428721/71‑1 (MQ=255)
                               gAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCtt                                                       >  1:3798422/1‑53 (MQ=255)
                               gAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCtt                                                       >  1:1187813/1‑53 (MQ=255)
                                aaGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGa                                     <  1:1948168/70‑1 (MQ=255)
                                aaGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGa                                     <  1:704330/70‑1 (MQ=255)
                                aaGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAg                                    >  1:3343282/1‑71 (MQ=255)
                                aaGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAg                                    >  1:3231385/1‑71 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTc                                                         <  1:2881727/47‑1 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTc                                                         <  1:1798427/47‑1 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTc                                                         <  1:508133/47‑1 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTc                                                         <  1:1522549/47‑1 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTc                                                         <  1:1088762/47‑1 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTg                                                     <  1:996357/51‑1 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTg                                                     <  1:2497866/51‑1 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTg                                                     <  1:2087905/51‑1 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTg                                                     <  1:2295438/51‑1 (MQ=255)
                                   tCCCGGATAACTTCCTCTGCCAGGAATGCTCCCTCGGCAAAGACGTc                                                         <  1:1496566/47‑1 (MQ=255)
                                     ccGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTgaaga                                                 <  1:2764190/53‑1 (MQ=255)
                                     ccGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTgaaga                                                 <  1:1527518/53‑1 (MQ=255)
                                     ccGGATAACTTCCTCTGCCAGGAATGCTCCCTCGGCAAAGACGTCTTTgaaga                                                 <  1:1973356/53‑1 (MQ=255)
                                          tAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAg                          >  1:3135539/1‑71 (MQ=255)
                                           aaCTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCttt                                                      >  1:4023595/1‑42 (MQ=255)
                                           aaCTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAg                          >  1:714312/1‑70 (MQ=255)
                                           aaCTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAg                          >  1:3294313/1‑70 (MQ=255)
                                           aaCTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAg                          >  1:2923017/1‑70 (MQ=255)
                                             cTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGt                         <  1:296617/69‑1 (MQ=255)
                                             cTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGt                         <  1:1553436/69‑1 (MQ=255)
                                               tCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGt                         <  1:612268/67‑1 (MQ=255)
                                               tCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTaa                       >  1:1279720/1‑69 (MQ=255)
                                               tCCTCTGCCCGGAATGCTCCATCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTaa                       >  1:3305312/1‑69 (MQ=255)
                                                ccTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACgg                    >  1:112670/1‑71 (MQ=255)
                                                      cccGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATGGGAGGCAAAATGAGt                         >  1:1807656/1‑60 (MQ=255)
                                                      cccGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCAtt                <  1:3719068/69‑1 (MQ=255)
                                                      cccGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATtgt              >  1:3371230/1‑71 (MQ=255)
                                                      cccGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATtgt              >  1:973985/1‑71 (MQ=255)
                                                       ccGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATtgtg             <  1:652644/71‑1 (MQ=255)
                                                       ccGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATtgtg             <  1:3765357/71‑1 (MQ=255)
                                                       ccGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATtgtg             <  1:3749601/71‑1 (MQ=255)
                                                         ggAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGat           <  1:1594890/71‑1 (MQ=255)
                                                         ggAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGat           <  1:2026110/71‑1 (MQ=255)
                                                         ggAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGat           <  1:2343064/71‑1 (MQ=255)
                                                         ggAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGat           <  1:3124507/71‑1 (MQ=255)
                                                         ggAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGat           <  1:1099633/71‑1 (MQ=255)
                                                         ggAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGat           <  1:1055708/71‑1 (MQ=255)
                                                           aaTGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGatca         <  1:2195509/71‑1 (MQ=255)
                                                            aTGCTCTCTCGGCAAAGACGTCTTTGACGAACTGGCATCGGAGGCAAAATGAGCAACGGCATTGTGatca         <  1:3134739/70‑1 (MQ=255)
                                                                tCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCaa                                <  1:2115872/43‑1 (MQ=255)
                                                                tCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCaa                                <  1:3204561/43‑1 (MQ=255)
                                                                  ccTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCaa                                <  1:3273655/41‑1 (MQ=255)
                                                                  ccTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGATCATCGGTTc  >  1:4029227/1‑71 (MQ=255)
                                                                  |                                                                      
GCAGGACGTTGCGCCAGGAACGCCGTGGAGTGAAGTCCCGGATAACTTCCTCTGCCCGGAATGCTCCCTCGGCAAAGACGTCTTTGAAGAACTGGCATCGGAGGCAAAATGAGTAACGGCATTGTGATCATCGGTTC  >  minE/1822022‑1822158

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