Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,187 G→A 24.8% L364L (CTG→CTA yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,1870GA24.8% 267.0 / 88.1 153L364L (CTG→CTAyibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (50/65);  new base A (20/18);  total (70/83)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 163 total aligned reads displayed.

ATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAC  >  minE/2571122‑2571256
                                                                 |                                                                     
aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct                                                                     >  1:177315/1‑68 (MQ=255)
aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                                  <  1:247959/71‑1 (MQ=255)
 tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                                                  >  1:3844818/1‑70 (MQ=255)
  cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                                >  1:3046096/1‑71 (MQ=255)
  cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                                >  1:5194102/1‑71 (MQ=255)
   tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGTGGAAAGTGCTGCTAAGAAt                                                               >  1:5332413/1‑71 (MQ=255)
   tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                               >  1:2021940/1‑71 (MQ=255)
    aCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                                                                 >  1:5349075/1‑68 (MQ=255)
     cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGGTAAGaa                                                                <  1:10340/68‑1 (MQ=255)
     cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                                <  1:4346162/68‑1 (MQ=255)
     cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                                                                 <  1:3103857/67‑1 (MQ=255)
     cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                             >  1:4704239/1‑71 (MQ=255)
     cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                             >  1:3375289/1‑71 (MQ=255)
      gCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAc                                                            <  1:4215956/71‑1 (MQ=255)
       cccAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                                                           <  1:5510290/71‑1 (MQ=255)
           ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCt                                                       <  1:2763793/71‑1 (MQ=255)
           ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCt                                                       >  1:622532/1‑71 (MQ=255)
           ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCt                                                       <  1:2196252/71‑1 (MQ=255)
            gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg                                                      <  1:3993103/71‑1 (MQ=255)
              ggCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTggg                                                    >  1:1178487/1‑70 (MQ=255)
              ggCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTgg                                                     >  1:4524765/1‑70 (MQ=255)
              ggCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                    >  1:3191083/1‑71 (MQ=255)
                cGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                                               <  1:5425119/58‑1 (MQ=255)
                 ggTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGc                                                 <  1:1867771/71‑1 (MQ=255)
                 ggTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGc                                                 <  1:740455/71‑1 (MQ=255)
                   tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                               <  1:1072992/71‑1 (MQ=255)
                   tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                               >  1:5134461/1‑71 (MQ=255)
                   tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                               >  1:1102715/1‑71 (MQ=255)
                    ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                                                       >  1:978845/1‑62 (MQ=255)
                    ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                                               >  1:2106060/1‑70 (MQ=255)
                     tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTgg                                                     >  1:4999038/1‑63 (MQ=255)
                     tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                                              <  1:3043775/70‑1 (MQ=255)
                      aCTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTac                                            >  1:478737/1‑70 (MQ=255)
                      aCTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTac                                            >  1:4997278/1‑70 (MQ=255)
                      aCTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTAt                                            >  1:682438/1‑71 (MQ=255)
                         ccGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt                                          >  1:5092416/1‑70 (MQ=255)
                          cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGc                                                        >  1:3770922/1‑55 (MQ=255)
                            atcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCACTACCttt                                        <  1:1088495/67‑1 (MQ=255)
                            aCCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                                             >  1:589001/1‑48 (MQ=255)
                             tcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTg                                                  <  1:1535765/57‑1 (MQ=255)
                             tcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtc                                     <  1:5010097/70‑2 (MQ=255)
                             ccATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt                                     >  1:4091093/1‑71 (MQ=255)
                             ccATTTCAGCCATGTTTCGGGGATGCGGAAAGTGCTGCTAAGAATGACCAGc                                                        <  1:2560631/52‑1 (MQ=255)
                              cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGAcc                                                           >  1:1493427/1‑48 (MQ=255)
                              cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtct                                    >  1:2350337/1‑69 (MQ=255)
                               aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                      <  1:258812/52‑1 (MQ=255)
                               aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                   <  1:3148849/71‑1 (MQ=255)
                               aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                   <  1:1718755/71‑1 (MQ=255)
                               aTTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtctc                                   >  1:1381078/1‑68 (MQ=255)
                                tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                                >  1:3620312/1‑41 (MQ=255)
                                tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                                                >  1:4178161/1‑41 (MQ=255)
                                tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                    >  1:2274517/1‑53 (MQ=255)
                                tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                    >  1:1763955/1‑53 (MQ=255)
                                tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                   <  1:4417888/70‑1 (MQ=255)
                                tttCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCa                                                          <  1:284447/47‑1 (MQ=255)
                                  tCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                >  1:472049/1‑71 (MQ=255)
                                  tCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGAt                                >  1:5009621/1‑71 (MQ=255)
                                   cAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                <  1:956006/70‑1 (MQ=255)
                                    aGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCa                              <  1:1321618/71‑1 (MQ=255)
                                      ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt                                     <  1:2496535/62‑1 (MQ=255)
                                      ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                            <  1:1570369/71‑1 (MQ=255)
                                      ccATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtt                            <  1:2614570/71‑1 (MQ=255)
                                        aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                      <  1:5676362/43‑1 (MQ=255)
                                        aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                            >  1:3875052/1‑69 (MQ=255)
                                        aTGTTTCGGTGATGCGGAAAGTCCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                   >  1:1980416/1‑62 (MQ=255)
                                         tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                  <  1:4432622/62‑1 (MQ=255)
                                         tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg                                  <  1:353479/62‑1 (MQ=255)
                                         tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAttaa                          <  1:5176130/70‑3 (MQ=255)
                                          gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                                                           <  1:2647854/36‑1 (MQ=255)
                                          gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                                                           <  1:5259237/36‑1 (MQ=255)
                                            ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                   >  1:924756/1‑58 (MQ=255)
                                            ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg                                   >  1:4929438/1‑58 (MQ=255)
                                            ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                          >  1:3630509/1‑67 (MQ=255)
                                             tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                                                    >  1:2804004/1‑40 (MQ=255)
                                             tCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                          >  1:4006483/1‑66 (MQ=255)
                                             tCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCAtt                                              >  1:875924/1‑46 (MQ=255)
                                             tCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCAtt                                              >  1:1497352/1‑46 (MQ=255)
                                              cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCt                                          <  1:5112401/49‑1 (MQ=255)
                                              cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                     <  1:3658348/70‑1 (MQ=255)
                                              cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                     <  1:3664760/70‑1 (MQ=255)
                                              cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                     <  1:5288122/70‑1 (MQ=255)
                                              cGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggg                      >  1:4349559/1‑69 (MQ=255)
                                               ggTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                                                      <  1:2341009/36‑1 (MQ=255)
                                                gTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgagggggagg                  <  1:5671021/71‑1 (MQ=255)
                                                   aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                <  1:4818097/54‑1 (MQ=255)
                                                   aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                <  1:3269735/54‑1 (MQ=255)
                                                   aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGAt                                <  1:1691149/54‑1 (MQ=255)
                                                   aTGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt                            <  1:2620025/58‑1 (MQ=255)
                                                    tGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAtta                           <  1:4467704/58‑2 (MQ=255)
                                                     gCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAttt                                    <  1:5496116/48‑1 (MQ=255)
                                                       ggAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga                          <  1:4302414/56‑1 (MQ=255)
                                                         aaaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGAt          >  1:616329/1‑70 (MQ=255)
                                                          aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTa        >  1:2522013/1‑71 (MQ=255)
                                                          aaGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTa        >  1:2923462/1‑71 (MQ=255)
                                                             tgctgcTAAGAATGACCAGCTGGATGCATTATCttt                                        <  1:5405463/36‑1 (MQ=255)
                                                              gctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTgaggggg                     >  1:4843058/1‑54 (MQ=255)
                                                               ctgctAAGAATGACCAGCTGGATGCATTATCTTTa                                       >  1:5485272/1‑35 (MQ=255)
                                                               ctgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGa    >  1:4544095/1‑70 (MQ=255)
                                                                tgctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAc  >  1:4026501/1‑71 (MQ=255)
                                                                 gctAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCa                              <  1:724465/42‑1 (MQ=255)
                                                                 |                                                                     
ATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAGGGGGAGGGGGGAGATTATGGAAC  >  minE/2571122‑2571256

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: