Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,331,252 A→C 18.1% intergenic (‑157/+273) sodA ← / ← fdhD superoxide dismutase, Mn/formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,331,2520AC18.1% 133.3 / 43.4 83intergenic (‑157/+273)sodA/fdhDsuperoxide dismutase, Mn/formate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base A (40/28);  new base C (6/9);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 2.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

GATTATCAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATCGACAGCAC  >  minE/2331182‑2331320
                                                                      |                                                                    
gATTATCAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCaa                                                                      >  1:1359930/1‑71 (MQ=255)
 aTTATCAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCaa                                                                      <  1:1146457/70‑1 (MQ=255)
  ttATCAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAAcc                                                                    <  1:1637279/71‑1 (MQ=255)
   tATCAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCa                                                                   >  1:1154541/1‑71 (MQ=255)
    aTCAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCa                                                                   <  1:3797395/70‑1 (MQ=255)
      cAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAAc                                                                     >  1:2515366/1‑66 (MQ=255)
      cAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACt                                                                 >  1:3243860/1‑70 (MQ=255)
       aaTGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCa                                                                   >  1:3350491/1‑67 (MQ=255)
       aaTGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCa                                                                   >  1:3236261/1‑67 (MQ=255)
       aaTGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                                <  1:2836477/70‑1 (MQ=255)
       aaTGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                                <  1:84867/70‑1 (MQ=255)
       aaTGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                                <  1:3062736/70‑1 (MQ=255)
       aaTGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                                <  1:1780546/70‑1 (MQ=255)
         tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCa                                                                   >  1:3608548/1‑65 (MQ=255)
         tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAAt                                                             >  1:3225817/1‑71 (MQ=255)
           ccGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCACCCACTTAATg                                                            <  1:1563312/70‑1 (MQ=255)
           ccGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGc                                                           >  1:2728994/1‑71 (MQ=255)
            cGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGcc                                                          <  1:3620749/71‑1 (MQ=255)
            cGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGcc                                                          >  1:432873/1‑71 (MQ=255)
              tACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                        >  1:1083941/1‑71 (MQ=255)
              tACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                        >  1:1240776/1‑71 (MQ=255)
               aCTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                        <  1:407390/70‑1 (MQ=255)
                cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCACCCACTTAATGCCGTg                                                       <  1:3097943/70‑1 (MQ=255)
                cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                                       <  1:938802/70‑1 (MQ=255)
                cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGa                                                      >  1:2820442/1‑71 (MQ=255)
                 ttttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                                <  1:3049609/60‑1 (MQ=255)
                 ttttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACt                                                                 <  1:475927/59‑1 (MQ=255)
                   ttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCaa                                                                      >  1:3452070/1‑52 (MQ=255)
                   ttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCa                                                                   >  1:3544020/1‑55 (MQ=255)
                   ttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCa                                                                   >  1:563977/1‑55 (MQ=255)
                   ttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGAt                                                     >  1:2529322/1‑69 (MQ=255)
                     cGTAAGGGTATGGTTTTGCAGGAAAATGCCCTAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGc                                                  <  1:3141048/70‑1 (MQ=255)
                      gTAAGGGTATGGTTTTGCATGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTg                                                   <  1:3010853/68‑1 (MQ=255)
                      gTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCAGTGATTg                                                   <  1:3434381/68‑1 (MQ=255)
                        aaGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAg                                               >  1:1109911/1‑70 (MQ=255)
                         aGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGt                                              >  1:1317595/1‑70 (MQ=255)
                           ggTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaaa                                           <  1:1317703/71‑1 (MQ=255)
                            gTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAAt                                          >  1:4009271/1‑71 (MQ=255)
                            gTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCACCAACTTAATGCCGTGATTGCCAGTAAAt                                          >  1:3491721/1‑71 (MQ=255)
                             tATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                         <  1:3321874/71‑1 (MQ=255)
                              aTGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                        <  1:3715981/55‑1 (MQ=255)
                               tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCACCCACtt                                                                <  1:204347/46‑1 (MQ=255)
                                ggTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGa                                       >  1:1164973/1‑70 (MQ=255)
                                 gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCACCCACtt                                                                <  1:2465117/44‑1 (MQ=255)
                                 gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCACCCACTTaa                                                              >  1:2362358/1‑46 (MQ=255)
                                 gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGAtt                                                    >  1:3430511/1‑56 (MQ=255)
                                 gTTTTGCAGGAAAATGCCCGAGAAGTGAAGCAAATCACCCACtt                                                                <  1:741845/44‑1 (MQ=38)
                                    ttGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACt                                                                 <  1:548003/40‑1 (MQ=255)
                                    ttGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACt                                                                 <  1:724235/40‑1 (MQ=255)
                                    ttGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACt                                                                 <  1:1754542/40‑1 (MQ=255)
                                    ttGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTaa                                                              >  1:943259/1‑43 (MQ=255)
                                       cAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaaa                                           <  1:2453899/59‑1 (MQ=255)
                                         ggAAAATGCCCGAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                             <  1:2588241/71‑1 (MQ=255)
                                         ggAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                         >  1:767844/1‑59 (MQ=255)
                                         ggAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAAcggcgg                              >  1:564955/1‑70 (MQ=255)
                                          gAAAATGCTCGAGATGTGAAGCAAATCACCCACTTAATGCCGTg                                                       <  1:1852711/44‑1 (MQ=37)
                                          gAAAATGCCCGAGATGTGAAGCAAATCAACCACTTaa                                                              >  1:3527892/1‑37 (MQ=255)
                                          gAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                        <  1:2185982/43‑1 (MQ=255)
                                          gAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGAt                                                     >  1:1121254/1‑46 (MQ=255)
                                          gAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACaa                                    >  1:3235466/1‑63 (MQ=255)
                                          gAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                             >  1:4086578/1‑70 (MQ=255)
                                           aaaaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCaa                           >  1:2351363/1‑71 (MQ=255)
                                            aaaTGCCCGAGATGTGAAGCAAATCAACCACTTaa                                                              >  1:1074531/1‑35 (MQ=255)
                                            aaaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACaa                                    >  1:525720/1‑61 (MQ=255)
                                             aaTGCCCGAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCaa                           >  1:3126002/1‑69 (MQ=255)
                                             aaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACa                         <  1:598756/71‑1 (MQ=255)
                                               tGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGAt                                                     >  1:3609962/1‑41 (MQ=255)
                                               tGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGAt                                                     >  1:1056999/1‑41 (MQ=255)
                                               tGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAgg                       >  1:851087/1‑71 (MQ=255)
                                               tGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAgg                       >  1:2735033/1‑71 (MQ=255)
                                                          tgAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTaa            >  1:1798860/1‑71 (MQ=255)
                                                          tgAAGCAAATCAACCACTTAATGCCGTGATTGCCa                                                >  1:782517/1‑35 (MQ=255)
                                                            aaGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATc          <  1:976579/71‑1 (MQ=255)
                                                            aaGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATc          <  1:411928/71‑1 (MQ=255)
                                                             aGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCaa                           >  1:1358328/1‑53 (MQ=255)
                                                              gCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGa                                       >  1:811371/1‑40 (MQ=255)
                                                              gCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATCg         >  1:1738474/1‑70 (MQ=255)
                                                              gCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATCGa        <  1:3510469/71‑1 (MQ=255)
                                                              gCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATCGa        >  1:3282303/1‑71 (MQ=255)
                                                               cAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCCGCAACAGGCGAAAGGTTAATCGa        >  1:517104/1‑70 (MQ=255)
                                                                 aaTCACCCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATCGACAg     <  1:156011/71‑1 (MQ=255)
                                                                  aTCACCCACTTAATGCCGTGATTGCCAGTAAATCGAc                                      <  1:4078321/37‑1 (MQ=255)
                                                                   tCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATCGACAGCa   >  1:2156980/1‑71 (MQ=255)
                                                                    cAACCACTTAATGCCGTGATTGCCAGTAAATCGACAAcggc                                <  1:2878447/41‑1 (MQ=255)
                                                                    cAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACa                         <  1:593823/48‑1 (MQ=255)
                                                                    cAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATCGACAGCAc  >  1:1313065/1‑71 (MQ=255)
                                                                     aaCCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATCGACAGCAc  <  1:2697115/70‑1 (MQ=255)
                                                                      |                                                                    
GATTATCAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATCGACAGCAC  >  minE/2331182‑2331320

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: