Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,331,736 C→T 6.9% S208N (AGC→AAC)  fdhD ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,331,7360CT 6.9% 373.8 / 13.0 115S208N (AGC→AAC) fdhDformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base C (66/41);  new base T (4/4);  total (70/45)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

Only 100 of 118 total aligned reads displayed.

CACATCGCCGACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCTTCATGCCC  >  minE/2331667‑2331806
                                                                     |                                                                      
cacaTCGCCGACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCt                                                                       >  1:3451955/1‑71 (MQ=255)
 acaTCGCCGACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCtt                                                                      >  1:3574946/1‑71 (MQ=255)
    tCGCCGACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTcc                                                                   >  1:1772590/1‑71 (MQ=255)
     cGCCGACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTccc                                                                  >  1:3813449/1‑71 (MQ=255)
     cGCCGACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTccc                                                                  >  1:1915779/1‑71 (MQ=255)
        cGACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCctt                                                               >  1:361760/1‑71 (MQ=255)
        cGACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCctt                                                               >  1:1835755/1‑71 (MQ=255)
         gACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCcttc                                                              <  1:3148712/71‑1 (MQ=255)
           cTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCcttctt                                                            <  1:276145/71‑1 (MQ=255)
           cTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCcttctt                                                            >  1:145558/1‑71 (MQ=255)
            ttttGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTcc                                                                   >  1:449595/1‑63 (MQ=255)
            ttttGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTtg                                                           >  1:542750/1‑71 (MQ=255)
               tGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTtgtg                                                         <  1:906628/70‑1 (MQ=255)
               tGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGa                                                        >  1:1290177/1‑71 (MQ=255)
               tGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGa                                                        >  1:435656/1‑71 (MQ=255)
                gCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCctt                                                               >  1:2664594/1‑63 (MQ=255)
                   ccATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCcttc                                                              <  1:762105/61‑1 (MQ=255)
                    cATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCttt                                                                     >  1:799525/1‑53 (MQ=255)
                      tttCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGACAGCTATCCCCTTCTTGTGAc                                                       >  1:1638786/1‑65 (MQ=255)
                       ttCATAACTGTCACGACTGGAAACCAGTACCGCACCTTGCTGCCATCTTTCCCCTTCTTGTGACcggcgg                                                 <  1:2436186/70‑1 (MQ=255)
                       ttCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACcgg                                                    <  1:370275/67‑1 (MQ=255)
                        tCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGcc                                               >  1:277272/1‑71 (MQ=255)
                         cATAACTGGCACGACTGGAAACCAGTACTGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGAcc                                                      >  1:1048777/1‑63 (MQ=255)
                            atcTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTa                                             <  1:2036577/67‑1 (MQ=255)
                            aaCTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTcccc                                                                 >  1:3614803/1‑49 (MQ=255)
                                ggCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCcttc                                                              <  1:908145/48‑1 (MQ=255)
                                 gCACGACTGGAAACCAGTACCTCACCTTGCTGCCAGCTTTCCCCTTCTTGTGAc                                                       >  1:3313243/1‑54 (MQ=255)
                                 gCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGa                                                        >  1:2332240/1‑53 (MQ=255)
                                 gCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGAc                                                       >  1:1841499/1‑54 (MQ=255)
                                 gCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGAc                                                       >  1:89859/1‑54 (MQ=255)
                                 gCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGttt                                      <  1:1457689/71‑1 (MQ=255)
                                 gCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGttt                                      >  1:1849901/1‑71 (MQ=255)
                                 gCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGttt                                      <  1:774381/71‑1 (MQ=255)
                                 gCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGtt                                       <  1:3905759/70‑1 (MQ=255)
                                  cACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTg                                     <  1:2174037/71‑1 (MQ=255)
                                  aaCGACTGGAAACCAGTACCGCACCTTGCTGCCAGTTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTg                                     <  1:2022655/70‑1 (MQ=255)
                                    cGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCcttc                                                              <  1:2468088/44‑1 (MQ=255)
                                    cGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCcttc                                                              <  1:1550522/44‑1 (MQ=255)
                                    cGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTtgtg                                                         <  1:1105628/49‑1 (MQ=255)
                                      aCTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTTTGa                                                        >  1:1977239/1‑48 (MQ=255)
                                      aCTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGAc                                                       >  1:600320/1‑49 (MQ=255)
                                      aCTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGAc                                                       >  1:1117043/1‑49 (MQ=255)
                                      aCTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGAc                                                       >  1:1773643/1‑49 (MQ=255)
                                      aCTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGt                                        <  1:973454/64‑1 (MQ=255)
                                      aCTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCa                                 <  1:422822/71‑1 (MQ=255)
                                        tGGAAACCAGTACCGCACCTTGCTGCCAGTTTTCCCcttct                                                             >  1:1047021/1‑41 (MQ=255)
                                        tGGAAACCAGTACCGCACCTTGCTGCCAGTTTTCCCcttct                                                             >  1:2400467/1‑41 (MQ=255)
                                        tGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAgc                               >  1:2333392/1‑71 (MQ=255)
                                         ggAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAgcg                              >  1:3807887/1‑71 (MQ=255)
                                          gAAACCATTACCGCACCTTGCTGCCAGCTTTCCCcttc                                                              <  1:1927279/38‑1 (MQ=255)
                                          gAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCcttc                                                              <  1:737948/38‑1 (MQ=255)
                                          gAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAgc                               <  1:3560176/69‑1 (MQ=255)
                                           aaaCCAGTACCGCACCTTGCTGCCAGCTTTCCCcttct                                                             <  1:2734029/38‑1 (MQ=255)
                                           aaaCCAGTACCGCACCTTGCTGCCAGCTTTCCCcttct                                                             <  1:1258787/38‑1 (MQ=255)
                                           aaaCCAGTACCGCACCTTGCTGCCAGCTTTCCCcttct                                                             <  1:2916118/38‑1 (MQ=255)
                                           aaaCCAGTACCGCACCTTGCTGCCAGCTTTCCCcttct                                                             <  1:518563/38‑1 (MQ=255)
                                             aCCAGTACCGCACCTTGCTGCCAGTTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGACCAGCGCca                           >  1:1000365/1‑70 (MQ=255)
                                             aCCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCca                           >  1:989709/1‑70 (MQ=255)
                                             aCCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCca                           >  1:1186005/1‑70 (MQ=255)
                                             aCCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCca                           >  1:825658/1‑70 (MQ=255)
                                               cAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGtt                                       >  1:4043930/1‑56 (MQ=255)
                                               cAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACAt                        <  1:363817/71‑1 (MQ=255)
                                               cAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACAt                        <  1:2201815/71‑1 (MQ=255)
                                                aGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCggcggc                                                >  1:3178384/1‑46 (MQ=255)
                                                aGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTaa                                            >  1:3008798/1‑50 (MQ=255)
                                                 gTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATgg                      >  1:3447728/1‑71 (MQ=255)
                                                 gTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATgg                      >  1:1540043/1‑71 (MQ=255)
                                                 gTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATg                       >  1:2412415/1‑70 (MQ=255)
                                                  tACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGtt                                       >  1:444952/1‑53 (MQ=255)
                                                  tACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGt                                        <  1:761509/52‑1 (MQ=255)
                                                  tACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCa                                 <  1:4017574/59‑1 (MQ=255)
                                                   aCCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGc                     >  1:283724/1‑70 (MQ=255)
                                                    ccGCACGTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGttt                                      >  1:3884762/1‑52 (MQ=255)
                                                    ccGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGttt                                      >  1:1405225/1‑52 (MQ=255)
                                                    ccGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCcaca                         >  1:1927488/1‑65 (MQ=255)
                                                    ccGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCg                    <  1:3353191/70‑1 (MQ=255)
                                                    ccGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGa                   <  1:2013326/71‑1 (MQ=255)
                                                    ccGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCAGCCTAACAGTTTGTCCAGCGCCACATGGCg                    <  1:2637819/70‑1 (MQ=255)
                                                     cGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGAc                  >  1:2672942/1‑71 (MQ=255)
                                                      gCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGAc                  >  1:2320481/1‑70 (MQ=255)
                                                      gCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGAc                  >  1:3233027/1‑70 (MQ=255)
                                                      gCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGAc                  >  1:3057376/1‑70 (MQ=255)
                                                       cACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGttt                                      <  1:1960179/49‑1 (MQ=255)
                                                       cACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACAt                        >  1:2347344/1‑63 (MQ=255)
                                                          cTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAg                                         <  1:3858869/43‑1 (MQ=255)
                                                            tgctgcCTGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCca                           <  1:3441692/55‑1 (MQ=255)
                                                            tgctgcCAGTTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTc           >  1:2163415/1‑71 (MQ=255)
                                                            tgctgcCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGtt                                       >  1:2662725/1‑43 (MQ=255)
                                                              ctgcCAGTTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCtt         <  1:3246759/71‑1 (MQ=255)
                                                              ctgcCAGTTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCtt         <  1:3685717/71‑1 (MQ=255)
                                                                gcCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGt                                        <  1:1634702/38‑1 (MQ=255)
                                                                 ccAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCTTCa       <  1:2487443/70‑1 (MQ=255)
                                                                 ccAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCTTCa       <  1:3621604/70‑1 (MQ=255)
                                                                 ccAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCTTCAt      >  1:535382/1‑71 (MQ=255)
                                                                 ccAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCTTCAt      >  1:2474336/1‑71 (MQ=255)
                                                                   aGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTc                                   >  1:1545011/1‑40 (MQ=255)
                                                                   aGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTc                                   >  1:1434210/1‑40 (MQ=255)
                                                                     cTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCTTCATGccc  >  1:1811461/1‑71 (MQ=255)
                                                                     cTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCTTCATGcc   >  1:472672/1‑70 (MQ=255)
                                                                     cTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCTTCATGc    <  1:2007397/69‑1 (MQ=255)
                                                                     |                                                                      
CACATCGCCGACTTTTGCACCATTTCATAACTGGCACGACTGGAAACCAGTACCGCACCTTGCTGCCAGCTTTCCCCTTCTTGTGACCGGCGGCCTAACAGTTTGTCCAGCGCCACATGGCGACCCACGTCTTCATGCCC  >  minE/2331667‑2331806

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: