Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 982,617 C→T 24.7% D21D (GAC→GAT ydfA → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE982,6170CT24.7% 284.4 / 93.6 162D21D (GAC→GATydfAhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (54/68);  new base T (15/25);  total (69/93)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 168 total aligned reads displayed.

AGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGG  >  minE/982547‑982686
                                                                      |                                                                     
aGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAc                                                                       >  1:560113/1‑71 (MQ=255)
 ggTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACg                                                                      <  1:254718/71‑1 (MQ=255)
     aGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAc                                                                  >  1:4293902/1‑71 (MQ=255)
      gCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACg                                                                 <  1:1306891/71‑1 (MQ=255)
       cATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGt                                                                >  1:1875911/1‑71 (MQ=255)
         tGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTc                                                              <  1:3036222/71‑1 (MQ=255)
           gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                                                            >  1:5263158/1‑71 (MQ=255)
           gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                                                            <  1:2339448/71‑1 (MQ=255)
           gATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAc                                                            <  1:2359212/71‑1 (MQ=255)
            aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGAc                                                                       >  1:767283/1‑59 (MQ=255)
            aTACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCAcc                                                           >  1:93754/1‑71 (MQ=255)
             tACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCg                                                          <  1:1743708/71‑1 (MQ=255)
               cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                         <  1:165066/70‑1 (MQ=255)
               cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                         <  1:3261743/70‑1 (MQ=255)
               cTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                         <  1:4439219/70‑1 (MQ=255)
                tATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCa                                                             >  1:5086101/1‑65 (MQ=255)
                aaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGTAACGTTCACTGcc                                                        >  1:2488450/2‑69 (MQ=255)
                aaTCGACCTCGGCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGcca                                                       >  1:5656746/2‑69 (MQ=255)
                 aTCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                      >  1:4801204/1‑71 (MQ=255)
                  tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGtt                                                               >  1:103855/1‑61 (MQ=255)
                  tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                      >  1:2985504/1‑70 (MQ=255)
                  tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                      >  1:2598837/1‑70 (MQ=255)
                  tCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAAt                                                      >  1:5238324/1‑70 (MQ=255)
                   cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATg                                                     >  1:2933886/1‑70 (MQ=255)
                   cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGa                                                    >  1:2927463/1‑71 (MQ=255)
                   cGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGa                                                    >  1:3132958/1‑71 (MQ=255)
                    gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                   >  1:1627117/1‑71 (MQ=255)
                    gATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGAc                                                   >  1:486061/1‑71 (MQ=255)
                       ctcggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACg                                                  <  1:3665795/66‑1 (MQ=255)
                       cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                >  1:4769378/1‑71 (MQ=255)
                       cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                >  1:4194718/1‑71 (MQ=255)
                       cTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                >  1:2658389/1‑71 (MQ=255)
                         cggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGa                                                    >  1:2261816/2‑65 (MQ=255)
                         cggCAACAACGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGa                                                    >  1:1535482/2‑65 (MQ=255)
                          ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCAGGTTCACCGCAATGACGTATAc                                             >  1:5165757/1‑71 (MQ=255)
                          ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                  <  1:5432214/66‑1 (MQ=255)
                          ggCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                  <  1:1106693/66‑1 (MQ=255)
                           gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                            <  1:4958940/71‑1 (MQ=255)
                           gCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                            <  1:5574020/71‑1 (MQ=255)
                            caacaaCGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGc                                                         <  1:342520/57‑1 (MQ=255)
                            caacaaCGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATAcc                                            >  1:4368369/1‑70 (MQ=255)
                            caacaaCGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTc                                                              >  1:1368259/1‑52 (MQ=255)
                             aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                          <  1:4636744/71‑1 (MQ=255)
                             aacaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaa                                          <  1:587474/71‑1 (MQ=255)
                             aacaacGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCaa                                          <  1:3829569/71‑1 (MQ=255)
                             aacaacGAATCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCaa                                          <  1:5612030/71‑1 (MQ=255)
                              acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaaa                                         <  1:1357286/71‑1 (MQ=255)
                              acaacGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCaaa                                         <  1:3543836/71‑1 (MQ=255)
                                   gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                    <  1:348569/71‑1 (MQ=255)
                                   gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                    <  1:411584/71‑1 (MQ=255)
                                   gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                    <  1:3691746/71‑1 (MQ=255)
                                   gAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                    <  1:147057/71‑1 (MQ=255)
                                    aaTCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaa                                    >  1:5408693/1‑70 (MQ=255)
                                     aTCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCaaaa                                  <  1:4875973/71‑1 (MQ=255)
                                     aTCCCTGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCaaaa                                  <  1:5534185/71‑1 (MQ=255)
                                      tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGc                                                                    <  1:672586/36‑1 (MQ=255)
                                      tctcTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGc                                      <  1:4885688/66‑1 (MQ=255)
                                      tccctGGTGTGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCaa                                    <  1:3567720/65‑1 (MQ=255)
                                       ctctGGTGTACGTCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                       >  1:693582/1‑48 (MQ=255)
                                       ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                                   >  1:2412398/1‑68 (MQ=255)
                                       ctctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaa                                   >  1:1740954/1‑68 (MQ=255)
                                        tctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGtata                                              >  1:1169894/1‑56 (MQ=255)
                                        tctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGtata                                              >  1:1999018/1‑56 (MQ=255)
                                         ctGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTAAACCAAAAGCAAAACGtt                              >  1:476508/1‑71 (MQ=255)
                                            gtgtGCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                           <  1:5040189/71‑1 (MQ=255)
                                            gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                           <  1:5059234/71‑1 (MQ=255)
                                            gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                           <  1:3045327/71‑1 (MQ=255)
                                            gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                           <  1:3155242/71‑1 (MQ=255)
                                            gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaa                           <  1:4165288/71‑1 (MQ=255)
                                            gtgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTa                            <  1:5164881/70‑1 (MQ=255)
                                             tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaa                          <  1:3289837/71‑1 (MQ=255)
                                             tgtgCGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaa                           >  1:1653029/1‑70 (MQ=255)
                                             tgtACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTaaa                          <  1:5295608/71‑1 (MQ=255)
                                               tACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                  <  1:3186065/45‑1 (MQ=255)
                                                 cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCaa                                                       >  1:1703812/1‑38 (MQ=255)
                                                 cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCg                      <  1:2057759/71‑1 (MQ=255)
                                                 cGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCg                      <  1:2442740/71‑1 (MQ=255)
                                                 cGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCg                      <  1:2175992/71‑1 (MQ=255)
                                                 cGGCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCg                      <  1:1420245/71‑1 (MQ=255)
                                                  ggCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                  <  1:3771331/42‑1 (MQ=255)
                                                  ggCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACg                                                  <  1:2668792/42‑1 (MQ=255)
                                                  ggCGTGTTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGa                     <  1:2237662/71‑1 (MQ=255)
                                                    cGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaa                   <  1:1957919/71‑1 (MQ=255)
                                                      tgtTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGt                               <  1:3731539/57‑1 (MQ=255)
                                                      tgtTCCCCAACCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTaaaa                         >  1:2095174/1‑63 (MQ=255)
                                                        tttCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGttt                             >  1:4496144/1‑57 (MQ=255)
                                                          tCCAAACCAGGACGGCACGTTCACCGCAATGACGta                                                >  1:1068631/1‑36 (MQ=255)
                                                            cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCaaaa                                  >  1:1941308/1‑48 (MQ=255)
                                                            cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc           >  1:4718072/1‑71 (MQ=255)
                                                            cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGcc           >  1:578358/1‑71 (MQ=255)
                                                            cAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGc            <  1:1170512/70‑1 (MQ=255)
                                                             caaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACAttt                             >  1:3571703/2‑52 (MQ=255)
                                                              aaCCAGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGcgcg         <  1:4637559/71‑1 (MQ=255)
                                                              aaCCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACg                                <  1:2747890/48‑1 (MQ=255)
                                                                  aGGATGGAACGTTCACTGCCATGACGTATACCAAAAGCAAAACATTTAAAACCGAAACTGGTGCGCGCCga     <  1:3936469/71‑2 (MQ=255)
                                                                  aGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                    <  1:4531709/56‑1 (MQ=255)
                                                                  aGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                    <  1:2312492/56‑1 (MQ=255)
                                                                    gACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaa                    <  1:667521/54‑1 (MQ=255)
                                                                     aCGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCgctgg  <  1:3627447/71‑1 (MQ=255)
                                                                      cGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGaaaa                  >  1:4481924/1‑54 (MQ=255)
                                                                      |                                                                     
AGGTCAGCATGGATACTATCGATCTTGGCAACAACGAATCTCTGGTGTACGGCGTGTTTCCAAACCAGGACGGCACGTTCACCGCAATGACGTATACCAAAAGCAAAACGTTTAAAACCGAAAATGGTGCCCGTCGCTGG  >  minE/982547‑982686

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: