Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,331,252 A→C 28.4% 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,2520AC28.4% 162.1 / 90.3 141intergenic (‑157/+273)sodA/fdhDsuperoxide dismutase, Mn/formate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base A (22/79);  new base C (31/9);  total (53/88)
Fisher's exact test for biased strand distribution p-value = 1.23e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.04e-01

Only 100 of 141 total aligned reads displayed.

TATCAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATC  >  minE/2331185‑2331312
                                                                   |                                                            
    aaTGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCAc                                                          >  1:2311238/1‑68 (MQ=255)
      tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                        <  1:2812575/68‑1 (MQ=255)
      tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                        <  1:234/68‑1 (MQ=255)
      tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                        <  1:637722/68‑1 (MQ=255)
      tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                        <  1:151071/68‑1 (MQ=255)
      tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                        <  1:2025150/68‑1 (MQ=255)
      tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                        <  1:952468/68‑1 (MQ=255)
      tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                        <  1:1456733/68‑1 (MQ=255)
      tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                        <  1:375570/68‑1 (MQ=255)
      tGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACtt                                                        <  1:848861/68‑1 (MQ=255)
       gCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTaa                                                      >  1:3067767/1‑69 (MQ=255)
             cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:176998/69‑1 (MQ=255)
             cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:1565447/69‑1 (MQ=255)
             cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:3202349/69‑1 (MQ=255)
             cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:30908/69‑1 (MQ=255)
             cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:1152967/69‑1 (MQ=255)
             cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:2317227/69‑1 (MQ=255)
             cTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:2939990/69‑1 (MQ=255)
             cTTTTCGTAAGGGTATGGTTTCGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:1357343/69‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:2029063/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:568305/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:1839011/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:642152/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:1803176/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:100743/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:69310/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:1556872/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:758520/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:89737/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:1229087/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:933232/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAAACACTTAATGCCGTg                                               <  1:631347/68‑1 (MQ=255)
               tttCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGAAGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:2695935/68‑1 (MQ=255)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATg                                                    >  1:489358/1‑53 (MQ=38)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:2853072/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:1011868/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:2454609/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:3030314/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:3040938/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:3432328/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:1966751/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:1488380/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:2061365/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTaa                                    >  1:929521/1‑69 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTa                                     >  1:2721625/1‑68 (MQ=39)
                         gTATGGTTTTGCAGGAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTa                                     >  1:1136613/1‑68 (MQ=39)
                           aTGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaaa                                   >  1:1111667/1‑68 (MQ=255)
                           aTGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaaa                                   >  1:1003128/1‑68 (MQ=255)
                            tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                 <  1:3078142/69‑1 (MQ=255)
                            tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                 <  1:2250589/69‑1 (MQ=255)
                            tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                 <  1:1053948/69‑1 (MQ=255)
                            tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                 <  1:2240782/69‑1 (MQ=255)
                            tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                 <  1:2024632/69‑1 (MQ=255)
                            tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                 <  1:2101830/69‑1 (MQ=255)
                            tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                 <  1:2551534/69‑1 (MQ=255)
                            tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                 <  1:406141/69‑1 (MQ=255)
                            tGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATc                                 <  1:737526/69‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaa                                    <  1:960750/64‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaa                                    <  1:3389071/64‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaa                                    <  1:729361/64‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaa                                    <  1:2619383/64‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaa                                    <  1:1307793/64‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaa                                    <  1:2450278/64‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaa                                    <  1:1845998/64‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaa                                    <  1:1955347/64‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaa                                    <  1:2147769/64‑1 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCg                                >  1:896158/1‑68 (MQ=255)
                              gTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCg                                >  1:2052217/1‑68 (MQ=255)
                                 ttGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTg                                               <  1:1728069/50‑1 (MQ=255)
                                       gAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAgcggcg                       >  1:2110040/1‑68 (MQ=38)
                                       gAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAgcggcg                       >  1:1524862/1‑68 (MQ=38)
                                       gAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAgcggcg                       >  1:1967790/1‑68 (MQ=38)
                                       gAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAgcggcg                       >  1:2133122/1‑68 (MQ=38)
                                       gAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAgcggcg                       >  1:330616/1‑68 (MQ=38)
                                       gAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAgcggcg                       >  1:3302624/1‑68 (MQ=38)
                                       gAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAgcggcg                       >  1:3109226/1‑68 (MQ=38)
                                       gAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAgcggcg                       >  1:3195261/1‑68 (MQ=38)
                                       gAAAATGCCCAAGATGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAGcggc                        >  1:139477/1‑67 (MQ=38)
                                        aaaaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTaaa                                   >  1:1462987/1‑55 (MQ=255)
                                        aaaaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                     >  1:3323141/1‑69 (MQ=255)
                                        aaaaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                     >  1:834821/1‑69 (MQ=255)
                                        aaaaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                     <  1:2211711/69‑1 (MQ=255)
                                        aaaaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                     <  1:2198574/69‑1 (MQ=255)
                                        aaaaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                     <  1:2757265/69‑1 (MQ=255)
                                        aaaaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                     <  1:3370406/69‑1 (MQ=255)
                                          aaTGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:2418787/40‑1 (MQ=255)
                                             gCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGt                                                <  1:2333504/37‑1 (MQ=255)
                                                   gaTGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAGCGGCGGCAACa                 <  1:3302953/62‑1 (MQ=255)
                                                   gaTGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAGCGGCGGCAACa                 <  1:1707818/62‑1 (MQ=255)
                                                   gaTGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAGCGGCGGCAACa                 <  1:2378630/62‑1 (MQ=255)
                                                   gaTGTGAAGCAAATCACCCACTTAATGCCGTGATTGCCAGTAAATCGGCAGCGGCGGCAACa                 <  1:333034/62‑1 (MQ=255)
                                                       tgAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGtt      >  1:1263073/1‑69 (MQ=255)
                                                       tgAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGtt      >  1:3015227/1‑69 (MQ=255)
                                                       tgAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGtt      >  1:1122582/1‑69 (MQ=255)
                                                         aaGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTaa    >  1:216525/1‑69 (MQ=255)
                                                           gCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATc  <  1:907481/69‑1 (MQ=255)
                                                           gCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATc  <  1:1233972/69‑1 (MQ=255)
                                                            cAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATc  <  1:667084/68‑1 (MQ=255)
                                                                 cAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                     >  1:1391333/1‑44 (MQ=255)
                                                                 cAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACggcggc                     >  1:1807393/1‑44 (MQ=255)
                                                                   |                                                            
TATCAATGCCGTACTTTTCGTAAGGGTATGGTTTTGCAGGAAAATGCCCGAGATGTGAAGCAAATCAACCACTTAATGCCGTGATTGCCAGTAAATCGACAACGGCGGCAACAGGCGAAAGGTTAATC  >  minE/2331185‑2331312

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: