Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,411,465 C→T 25.9% S245S (AGC→AGT hemY → predicted protoheme IX synthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,411,4650CT25.9% 156.0 / 47.5 128S245S (AGC→AGThemYpredicted protoheme IX synthesis protein
Reads supporting (aligned to +/- strand):  ref base C (49/45);  new base T (20/13);  total (70/58)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

Only 100 of 128 total aligned reads displayed.

CTGGAACAACAGGCATGGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGTG  >  minE/2411409‑2411531
                                                        |                                                                  
cTGGAACAACAGGCATGGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGt                                                        <  1:2108678/69‑1 (MQ=255)
             cATTGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGaaaa                                           <  1:366624/69‑1 (MQ=255)
             cATGGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGaaaa                                           <  1:2784192/69‑1 (MQ=255)
             cATGGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGaaaa                                           <  1:2823960/69‑1 (MQ=255)
             cATGGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGaaaa                                           <  1:2777185/69‑1 (MQ=255)
             cATGGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGaaaa                                           <  1:3173085/69‑1 (MQ=255)
             cATGGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGaaaa                                           <  1:113314/69‑1 (MQ=255)
             cATGGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGaaaa                                           <  1:1879876/69‑1 (MQ=255)
             cATGGATTGCCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGaaaa                                           <  1:1777360/69‑1 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGa                                                                  >  1:2929137/1‑41 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:2681819/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:2693148/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:1728841/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:2729228/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:2733275/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:2511336/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:2297374/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:2068896/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:2925855/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:728277/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:899361/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:98576/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCaa                                      >  1:990673/1‑69 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCa                                       >  1:2932164/1‑68 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCa                                       >  1:683269/1‑68 (MQ=255)
                  aTTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGAGaa                                                                 >  1:55033/1‑42 (MQ=255)
                   ttGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCaaa                                     <  1:787078/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:3376984/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:678038/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:268615/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:2140783/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:997223/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:1053123/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:1282067/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:1283686/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:190249/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:1748947/69‑1 (MQ=255)
                    tGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGTGAAGGTTTGCGTAACTGGTGGAAAAACCAAAg                                    <  1:1795538/69‑1 (MQ=255)
                             ggATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAAc                           >  1:2008241/1‑69 (MQ=255)
                             ggATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAAc                           >  1:1629711/1‑69 (MQ=255)
                              gATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTTCGTAACTGGTGGAAAAACCAAGCCGGAAAAcgc                         <  1:1259908/69‑1 (MQ=255)
                                tCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTTGTGGAAAAACCAAAGCCGGAAAAcgc                         >  1:3426587/1‑68 (MQ=255)
                                tCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAAcgc                         >  1:2992883/1‑68 (MQ=255)
                                tCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAAcgc                         >  1:1006389/1‑68 (MQ=255)
                                tCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAAcgc                         >  1:1358541/1‑68 (MQ=255)
                                tCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAAcgc                         >  1:3242146/1‑68 (MQ=255)
                                tCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAAcgc                         >  1:3052726/1‑68 (MQ=255)
                                tCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAAcgc                         <  1:1577247/68‑1 (MQ=255)
                                tCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAAcgc                         >  1:2089806/1‑68 (MQ=255)
                                     cgcgTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGtcatc                   <  1:2175564/69‑1 (MQ=255)
                                     cgcgTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGtcatc                   <  1:3356868/69‑1 (MQ=255)
                                     cgcgTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGtcatc                   >  1:2018083/1‑69 (MQ=255)
                                     cgcgTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGtcatc                   <  1:1437322/69‑1 (MQ=255)
                                           ccGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAg             >  1:3307426/1‑69 (MQ=255)
                                             gATAACGGTAGCGAATGTTTGCGTAACTGGTGGAAAAACCaaa                                     <  1:3409652/43‑1 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:433591/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:2184925/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:1367936/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:3150558/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:3130802/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:666058/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:3085538/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:693218/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:2981919/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:2915217/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:1499863/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:1831645/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:876911/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:1891713/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:1431235/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:2450002/1‑67 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTg        <  1:1014518/69‑1 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTg        <  1:2874620/69‑1 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTg        <  1:303168/69‑1 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTg        <  1:196789/69‑1 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTg        <  1:223121/69‑1 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTg        <  1:2704880/69‑1 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTg        <  1:3311901/69‑1 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTg        <  1:3319661/69‑1 (MQ=255)
                                                aaCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAAACAAAGCCGGAAAACGCGTCATCAGGTAGCGt          >  1:2055257/1‑67 (MQ=255)
                                                 aCGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTg        <  1:853884/68‑1 (MQ=255)
                                                   ggTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCa      <  1:1031164/68‑1 (MQ=255)
                                                   ggTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCa      <  1:2651319/68‑1 (MQ=255)
                                                   ggTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCa      <  1:769316/68‑1 (MQ=255)
                                                   ggTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCa      <  1:1220412/68‑1 (MQ=255)
                                                   ggTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCa      <  1:1574881/68‑1 (MQ=255)
                                                   ggTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCa      <  1:2257573/68‑1 (MQ=255)
                                                   ggTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCa      <  1:2568727/68‑1 (MQ=255)
                                                      aGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGt   >  1:1925303/1‑68 (MQ=255)
                                                      aGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGt   >  1:3120733/1‑68 (MQ=255)
                                                      aGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGt   >  1:2242513/1‑68 (MQ=255)
                                                      aGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGt   >  1:1241044/1‑68 (MQ=255)
                                                      aGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGt   >  1:3026851/1‑68 (MQ=255)
                                                      aGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGt   >  1:1479141/1‑68 (MQ=255)
                                                      aGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGt   >  1:2863584/1‑68 (MQ=255)
                                                      aGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGt   >  1:1543833/1‑68 (MQ=255)
                                                      aGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGt   >  1:950315/1‑68 (MQ=255)
                                                       gCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGTg  >  1:1080798/1‑68 (MQ=255)
                                                       gCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGTg  >  1:2314429/1‑68 (MQ=255)
                                                       gCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGTg  >  1:205093/1‑68 (MQ=255)
                                                        |                                                                  
CTGGAACAACAGGCATGGATTGGCCTGATGGATCAGGCGCGTGCCGATAACGGTAGCGAAGGTTTGCGTAACTGGTGGAAAAACCAAAGCCGGAAAACGCGTCATCAGGTAGCGTTGCAGGTG  >  minE/2411409‑2411531

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: