Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,353,113 A→C 17.8% F434L (TTT→TTG yihF ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,353,1130AC17.8% 181.6 / 34.6 101F434L (TTT→TTGyihFconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/56);  new base C (8/10);  total (35/66)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 117 total aligned reads displayed.

AAACAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTGCGCTTTTTCCGCA  >  minE/2353046‑2353179
                                                                   |                                                                  
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:2639475/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:956085/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:874444/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:769031/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:3322241/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:2836765/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:2801141/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:24572/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:1671696/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:1584856/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:1395419/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGa                                                                    <  1:1039083/68‑1 (MQ=255)
aaaCAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGATGTCAATCGTGATTAACTGa                                                                    <  1:3435225/68‑1 (MQ=255)
  acaacTTTACCCGGTGTATAACGCAGTTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGCAAc                                                                 >  1:1331800/1‑69 (MQ=255)
  acaacTTTACCCGGTGTATAACGCAGTTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGCAAc                                                                 >  1:2623888/1‑69 (MQ=255)
  acaacTTTACCCGGTGTATAACGCAGTTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGCAAc                                                                 >  1:3319074/1‑69 (MQ=255)
  acaacTTTACCCGGTGTATAACGCAGTTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGCAAc                                                                 >  1:2539628/1‑69 (MQ=255)
  acaacTTTACCCGGTGTATAACGCAGTTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGCAAc                                                                 >  1:353050/1‑69 (MQ=255)
  acaacTTTACCCGGTGTATAACGCAGTTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGCAAc                                                                 >  1:1087545/1‑69 (MQ=255)
  acaacTTTACCCGGTGTATAACGCAGTTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGCAAc                                                                 >  1:586501/1‑69 (MQ=255)
  acaacTTTACCCGGTGTAAAACGCAGTTCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGCAACa                                                                >  1:3082889/1‑69 (MQ=255)
       tttACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTg                                                            <  1:2342946/69‑1 (MQ=255)
            ccGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCt                                                        <  1:3033139/68‑1 (MQ=255)
            ccGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCt                                                        <  1:1651095/68‑1 (MQ=255)
            ccGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCt                                                        <  1:543363/68‑1 (MQ=255)
              ggTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAAt                                                     <  1:2286237/69‑1 (MQ=255)
              ggTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAAt                                                     <  1:2131906/69‑1 (MQ=255)
                tgtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTg                                                            <  1:1468599/60‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:1810677/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:1958004/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:1779831/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:2898168/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:115316/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:281840/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:2223495/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:383659/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:2475134/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTaa                                                      <  1:1604517/65‑1 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtc                                                  >  1:2984165/1‑69 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtc                                                  >  1:2381403/1‑69 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtc                                                  >  1:2361634/1‑69 (MQ=255)
                 gtATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGt                                                   >  1:1646045/1‑68 (MQ=255)
                    taACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtcat                                                <  1:232473/68‑1 (MQ=255)
                    taACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtcat                                                <  1:2181200/68‑1 (MQ=255)
                    taACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtcat                                                <  1:1347987/68‑1 (MQ=255)
                    taACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtcat                                                <  1:2733752/68‑1 (MQ=255)
                      aCGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtcat                                                >  1:572329/1‑66 (MQ=255)
                      aCGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtcat                                                >  1:3089474/1‑66 (MQ=255)
                      aCGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGtcat                                                >  1:910955/1‑66 (MQ=255)
                      aCGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTcatcat                                             >  1:1008086/1‑69 (MQ=255)
                      aCGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTcatcat                                             >  1:1067174/1‑69 (MQ=255)
                                  cGAGGCGGTATTGTTGTCAATCGTGATTAAGTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTg                                 <  1:523856/69‑1 (MQ=255)
                                  cGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTg                                 <  1:2721234/69‑1 (MQ=255)
                                  cGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTg                                 <  1:1275035/69‑1 (MQ=255)
                                  cGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTg                                 <  1:802862/69‑1 (MQ=255)
                                  cGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTg                                 <  1:2087723/69‑1 (MQ=255)
                                  cGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTg                                 <  1:2911029/69‑1 (MQ=255)
                                    aGGCTGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:252966/69‑1 (MQ=255)
                                    aGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:1190224/69‑1 (MQ=255)
                                    aGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:3069673/69‑1 (MQ=255)
                                    aGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:901054/69‑1 (MQ=255)
                                     ggCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGttt                              <  1:661100/69‑1 (MQ=255)
                                     ggCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGttt                              <  1:1739094/69‑1 (MQ=255)
                                     ggCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGttt                              <  1:1493328/69‑1 (MQ=255)
                                           aTTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGc                        >  1:1221455/1‑69 (MQ=255)
                                           aTTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGc                        >  1:2603352/1‑69 (MQ=255)
                                           aTTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGc                        >  1:1992609/1‑69 (MQ=255)
                                              gttgtCTATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTg                     >  1:441756/1‑69 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTaa                                                      >  1:2713190/1‑36 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTaa                                                      >  1:747480/1‑36 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:2980799/59‑1 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:541541/59‑1 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:328692/59‑1 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:1267005/59‑1 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:2759567/59‑1 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:3159886/59‑1 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTg                     >  1:1116874/1‑69 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTg                     >  1:3396769/1‑69 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTg                     >  1:3135151/1‑69 (MQ=255)
                                              gttgtCAATCGTGATTAACTGAAACATCTGACCTAATGTAATCATCCCGCTGATTTGtt                               <  1:1204933/59‑1 (MQ=255)
                                              gttgtCAATCCTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGtt                               <  1:88691/59‑1 (MQ=255)
                                                       cGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTgcgc            >  1:1607441/1‑69 (MQ=255)
                                                       cGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTgcgc            >  1:2627386/1‑69 (MQ=255)
                                                       cGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTgcgc            >  1:309248/1‑69 (MQ=255)
                                                       cGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTgcgc            >  1:2565101/1‑69 (MQ=255)
                                                              aaCTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTGCGCTTTTTc      >  1:3285482/1‑68 (MQ=255)
                                                                cTGCAACATCTGACCTAATGTCATCATCCCGCTGATTTTTTTATCAGCCCGCTTTTGCGCttttt       <  1:3026796/65‑1 (MQ=255)
                                                                cTGCAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCCGCTTTTGCGCttttt       <  1:1265111/65‑1 (MQ=255)
                                                                cTGCAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCCGCTTTTGCGCttttt       <  1:731085/65‑1 (MQ=255)
                                                                cTGCAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCCGCTTTTGCGCttttt       <  1:1292313/65‑1 (MQ=255)
                                                                cTGCAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCCGCTTTTGCGCttttt       <  1:494336/65‑1 (MQ=255)
                                                                cTGCAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCCGCTTTTGCGCttttt       <  1:3424626/65‑1 (MQ=255)
                                                                cTGCAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCCGCTTTTGCGCttttt       <  1:2602919/65‑1 (MQ=255)
                                                                cTGCAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCCGCTTTTGCGCttttt       <  1:318541/65‑1 (MQ=255)
                                                                cTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTGCGCTTTTTCCg    <  1:166311/68‑1 (MQ=255)
                                                                cTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTGCGCTTTTTCCg    <  1:2240645/68‑1 (MQ=255)
                                                                cTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTGCGCTTTTTCCg    <  1:808768/68‑1 (MQ=255)
                                                                cTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTGCGCTTTTTCCg    <  1:340806/68‑1 (MQ=255)
                                                                cTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTGCGCTTTTTCCg    <  1:3012229/68‑1 (MQ=255)
                                                                  gAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTGCGCTTTTTCCGCa  >  1:1549852/1‑68 (MQ=255)
                                                                   |                                                                  
AAACAACTTTACCCGGTGTATAACGCAGCTGCAGCGAGGCGGTATTGTTGTCAATCGTGATTAACTGAAACATCTGACCTAATGTCATCATCCCGCTGATTTGTTTATCAGCCTGCTTTTGCGCTTTTTCCGCA  >  minE/2353046‑2353179

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: