Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,730,976 C→T 5.2% V87V (GTG→GTA yfhB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,730,9760CT 5.2% 361.1 / 6.1 115V87V (GTG→GTAyfhBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (45/64);  new base T (2/4);  total (47/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

Only 100 of 116 total aligned reads displayed.

AACAGGTAGGTGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCCACAGAAGCA  >  minE/1730911‑1731043
                                                                 |                                                                   
aaCAGGTAGGTGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGa                                                                 <  1:195768/70‑1 (MQ=255)
 aCAGGTAGGTGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGaa                                                                <  1:396533/70‑1 (MQ=255)
 aCAGGTAGGTGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGa                                                                 <  1:3354926/69‑1 (MQ=255)
  cAGGTAGGTGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCACGAAACCAACGTACGAAAt                                                              <  1:4559766/71‑1 (MQ=255)
  cAGGTAGGTGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCACGAAACCAACGTACGAAAt                                                              <  1:2632336/71‑1 (MQ=255)
    ggtaggtGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCg                                                            <  1:3968556/71‑1 (MQ=255)
     gtaggtGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCg                                                            <  1:1869840/70‑1 (MQ=255)
       aggtGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCACGAAACCAACGTACGAAAt                                                              <  1:1983355/66‑1 (MQ=255)
       aggtGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGT‑‑CGAAACCAGCGCACGAAATCGGCCTg                                                       >  1:1793219/1‑71 (MQ=255)
         gtggtTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACg                                                                  <  1:3938748/60‑1 (MQ=255)
         gtggtTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACg                                                                  <  1:1758789/60‑1 (MQ=255)
         gtggtTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGaa                                                                <  1:1393675/62‑1 (MQ=255)
         gtggtTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGc                                                          <  1:4449003/68‑1 (MQ=255)
          tggtTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGc                                                      <  1:4404904/71‑1 (MQ=255)
          tggtTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGc                                                      <  1:3823864/71‑1 (MQ=255)
            gtTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCa                                                     <  1:608872/70‑1 (MQ=255)
             ttAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAAc                                                   >  1:2308828/1‑71 (MQ=255)
             ttAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAAc                                                   <  1:4476514/71‑1 (MQ=255)
               aaTCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGAACGTc                                                >  1:4570563/1‑71 (MQ=255)
                aTCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTc                                                >  1:2230014/1‑71 (MQ=255)
                 tCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCt                                                        <  1:1487056/62‑1 (MQ=255)
                 tCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGc                                                      <  1:369365/64‑1 (MQ=255)
                 tCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGc                                                      <  1:3371368/64‑1 (MQ=255)
                  cGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAAc                                                   >  1:579856/1‑66 (MQ=255)
                   gCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAAt                                                              >  1:1592637/1‑54 (MQ=255)
                   gCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTg                                              >  1:2076465/1‑70 (MQ=255)
                      ccTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCa                                                                    >  1:2751103/1‑45 (MQ=255)
                       cTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTaa                                           <  1:3986754/69‑1 (MQ=255)
                         gAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGt                                       >  1:3702402/1‑71 (MQ=255)
                          aaCCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACg                                        <  1:4089535/69‑1 (MQ=255)
                             cAGCGGATAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCtt                                   >  1:525245/1‑71 (MQ=255)
                             cAGCGGAAAGGTGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCtt                                   >  1:3747770/1‑71 (MQ=255)
                             cAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCtt                                   >  1:4205628/1‑71 (MQ=255)
                               gCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAAt                                                              <  1:458470/42‑1 (MQ=255)
                               gCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGc                                     >  1:3109011/1‑67 (MQ=255)
                               gCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCg                                 <  1:1307877/71‑1 (MQ=255)
                               gCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGAACGTCTGTAAACGTGCTTCGc                                <  1:52963/71‑1 (MQ=255)
                                 ggAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCt                               <  1:2648884/71‑1 (MQ=255)
                                 ggAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCt                               <  1:112759/71‑1 (MQ=255)
                                  gAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCaa                                                    >  1:4394517/1‑49 (MQ=255)
                                  gAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCtg                              <  1:457597/71‑1 (MQ=255)
                                  gAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCtg                              <  1:2563332/71‑1 (MQ=255)
                                  gAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCtg                              >  1:597706/1‑71 (MQ=255)
                                   aaaGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGTCCTGCAACGTCTGTa                                            <  1:2282389/56‑1 (MQ=255)
                                   aaaGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTa                                            <  1:4199133/56‑1 (MQ=255)
                                   aaaGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCtgt                             <  1:2522904/71‑1 (MQ=255)
                                   aaaGGCGGTAACATTGTCACGAAACCAACGTACGAAATCGGCCTGCAACGCCTGCAAAc                                         >  1:1365486/1‑59 (MQ=255)
                                   aaaGGCGGTAACATTGTCACGAAACCAACGTACGAAATCGGCCTGCAACGCCTGCAAAc                                         >  1:1662872/1‑59 (MQ=255)
                                    aaGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAAc                                         >  1:4430056/1‑58 (MQ=255)
                                     aGGCGTTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGa                           >  1:2169385/1‑71 (MQ=255)
                                     aGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTaaa                                          <  1:2181783/56‑1 (MQ=255)
                                       gCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGAcc                         <  1:3355850/71‑1 (MQ=255)
                                       gCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGAcc                         <  1:3361439/71‑1 (MQ=255)
                                         ggTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaa                       <  1:2097716/71‑1 (MQ=255)
                                         ggTAACATTGTCGCGAAACCAGCCCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaa                       <  1:419404/71‑1 (MQ=255)
                                          gTAACATTGTCGCGAACCCAGCGCACGAAATCGGCCTGCAACGTCTGtt                                            <  1:1230132/49‑2 (MQ=255)
                                          gTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGt                                                 <  1:3420914/44‑1 (MQ=255)
                                          gTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTctct                               <  1:1348043/62‑1 (MQ=255)
                                          gTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCt                               <  1:2310990/62‑1 (MQ=255)
                                          gTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCt                               <  1:2316440/62‑1 (MQ=255)
                                          gTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaa                       <  1:226426/70‑1 (MQ=255)
                                          gTAACATTGTCGCGAAACCAGCGCACAAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaaa                      <  1:2370295/71‑1 (MQ=255)
                                           tAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaaaa                     <  1:2142580/71‑1 (MQ=255)
                                           tAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaaaa                     <  1:648364/71‑1 (MQ=255)
                                            aaCATTGTCGCGAAACCAGCGCACGAAATCGGCCTGc                                                      <  1:643913/37‑1 (MQ=255)
                                            aaCATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGt                                                 >  1:4223662/1‑42 (MQ=255)
                                            aaCATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTaaa                                          >  1:3270307/1‑49 (MQ=255)
                                            aaCATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTaa                                           >  1:4087585/1‑48 (MQ=255)
                                            aaCATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAg                    >  1:771745/1‑71 (MQ=255)
                                            aaCATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAg                    >  1:2895023/1‑71 (MQ=255)
                                             aCATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAg                    <  1:1433584/70‑1 (MQ=255)
                                             aCATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGt                   >  1:2285043/1‑71 (MQ=255)
                                               aTTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaa                       >  1:4403098/1‑65 (MQ=255)
                                                ttGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTa                                            <  1:3870381/43‑1 (MQ=255)
                                                ttGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCa                >  1:2740419/1‑71 (MQ=255)
                                                  gTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTaa                                           <  1:25514/42‑1 (MQ=255)
                                                  gTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAAc                                         >  1:2213488/1‑44 (MQ=255)
                                                  gTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCtg                              >  1:2040178/1‑55 (MQ=255)
                                                  gTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCAcc              >  1:3162841/1‑71 (MQ=255)
                                                  gTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCAcc              >  1:635963/1‑71 (MQ=255)
                                                  gTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCAcc              >  1:4170451/1‑71 (MQ=255)
                                                   tCGCGAAACCATCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCAccc             <  1:1322238/71‑1 (MQ=255)
                                                     gcgAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGAcc                         >  1:636868/1‑57 (MQ=255)
                                                     gcgAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaaaa                     >  1:4259676/1‑61 (MQ=255)
                                                     gcgAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaaaa                     >  1:1423753/1‑61 (MQ=255)
                                                      cgAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCca          <  1:2107634/71‑1 (MQ=255)
                                                       gAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaaa                      <  1:4503128/58‑1 (MQ=255)
                                                       gAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCca          >  1:1940133/1‑70 (MQ=255)
                                                       gAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCca          >  1:172215/1‑70 (MQ=255)
                                                       gAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCca          >  1:1443659/1‑70 (MQ=255)
                                                       gAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCca          >  1:3862642/1‑70 (MQ=255)
                                                       gAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTACCAAAAGTg                  <  1:159598/61‑1 (MQ=255)
                                                       gAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGGGACCAAAAGTGCACCCCca          >  1:177574/1‑70 (MQ=255)
                                                         aaCCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACg                                        <  1:343484/38‑1 (MQ=255)
                                                         aaCCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCca          >  1:2073484/1‑68 (MQ=255)
                                                         aaCCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCCACAg       <  1:514049/71‑1 (MQ=255)
                                                         aaCCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCCACAg       <  1:4480378/71‑1 (MQ=255)
                                                           ccAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCaaaa                     >  1:570675/1‑55 (MQ=255)
                                                            cAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCCACAGAAg    <  1:774149/71‑1 (MQ=255)
                                                              gcgcACGAAATCGTCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCCACAGAAGCa  <  1:3131968/71‑1 (MQ=255)
                                                                 |                                                                   
AACAGGTAGGTGGTTAATCGCTCCTGAACCAGCGGAAAGGCGGTAACATTGTCGCGAAACCAGCGCACGAAATCGGCCTGCAACGTCTGTAAACGTGCTTCGCTGTGACCAAAAGTGCACCCCCACAGAAGCA  >  minE/1730911‑1731043

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: