Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,571,166 T→C 30.8% V357V (GTT→GTC yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,571,1660TC30.8% 248.5 / 156.0 198V357V (GTT→GTCyibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (52/85);  new base C (37/24);  total (89/109)
Fisher's exact test for biased strand distribution p-value = 3.42e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 206 total aligned reads displayed.

TGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAG  >  minE/2571097‑2571233
                                                                     |                                                                   
tgatATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGttt                                                                    >  1:4241635/1‑71 (MQ=255)
      cGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGa                                                              >  1:222200/1‑71 (MQ=255)
       gATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGAt                                                             >  1:2979381/1‑71 (MQ=255)
         tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGc                                                           <  1:2141511/71‑1 (MQ=255)
         tGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGc                                                           >  1:4482994/1‑71 (MQ=255)
            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:1555635/71‑1 (MQ=255)
            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:3007160/71‑1 (MQ=255)
            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:2004031/71‑1 (MQ=255)
            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:3091865/71‑1 (MQ=255)
            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:1480555/71‑1 (MQ=255)
            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:5572070/71‑1 (MQ=255)
            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:109937/71‑1 (MQ=255)
            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:1087588/71‑1 (MQ=255)
            cTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:4756545/71‑1 (MQ=255)
             ttACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                                                        >  1:3570946/1‑68 (MQ=255)
              tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                                                        >  1:2958551/1‑67 (MQ=255)
              tACCAGACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCgca                                                        >  1:2442065/1‑67 (MQ=255)
                ccagACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGAt                                                             >  1:2823070/4‑62 (MQ=255)
                 ccGACGGCATCTACGCCCAGGTGTCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                <  1:3640913/58‑1 (MQ=255)
                 ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                <  1:3453689/58‑1 (MQ=255)
                 ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                <  1:4489800/58‑1 (MQ=255)
                 ccGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                <  1:101163/58‑1 (MQ=255)
                  cGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGtgc                                                  >  1:3225277/1‑71 (MQ=255)
                  agACGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc                                                  <  1:1571938/70‑1 (MQ=255)
                   gACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                <  1:1387534/56‑1 (MQ=255)
                     cGGTATTTACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                                               >  1:5093293/1‑69 (MQ=255)
                         aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTgctgct                                              >  1:1785017/1‑68 (MQ=255)
                         aTCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                           <  1:247959/71‑1 (MQ=255)
                          tCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAg                                           >  1:3844818/1‑70 (MQ=255)
                           cTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                         >  1:5194102/1‑71 (MQ=255)
                            tACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                                               >  1:335830/1‑62 (MQ=255)
                            tACGCCCAGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgctat                                               >  1:2071736/1‑62 (MQ=255)
                            tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGa                                                        <  1:3293001/55‑1 (MQ=255)
                            tACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAAt                                        >  1:1936429/1‑71 (MQ=255)
                            aaCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCg                                                                  >  1:3970743/2‑45 (MQ=255)
                             aCGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                                          >  1:5349075/1‑68 (MQ=255)
                              cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                         <  1:5279253/68‑1 (MQ=255)
                              cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                         <  1:4346162/68‑1 (MQ=255)
                              cGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGa                                          <  1:3103857/67‑1 (MQ=255)
                                cccAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                                    <  1:5510290/71‑1 (MQ=255)
                                   aGGTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGt                                                                <  1:3777817/40‑1 (MQ=37)
                                    ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGa                                                              >  1:261139/1‑41 (MQ=37)
                                    ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCt                                >  1:622532/1‑71 (MQ=255)
                                    ggTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCt                                <  1:2763793/71‑1 (MQ=255)
                                    ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCgg                                                         <  1:3166395/46‑1 (MQ=255)
                                    ggTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGt                                                    >  1:3931399/1‑51 (MQ=255)
                                     gTGGCGGTTTACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTg                               <  1:4966286/71‑1 (MQ=255)
                                     gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                <  1:2550300/38‑1 (MQ=255)
                                     gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                <  1:2095646/38‑1 (MQ=255)
                                     gTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGt                                                                <  1:5104158/38‑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)
                                           gTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGaa                                                       <  1:4706375/41‑1 (MQ=255)
                                            tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                        <  1:1072992/71‑1 (MQ=255)
                                            tttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAt                        >  1:5134461/1‑71 (MQ=255)
                                             ttACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCt                                >  1:4916584/1‑62 (MQ=255)
                                              tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTgg                              >  1:4999038/1‑63 (MQ=255)
                                              tACTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTa                      >  1:1495673/1‑71 (MQ=255)
                                              tACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCAtt                       <  1:3043775/70‑1 (MQ=255)
                                               aCTCTGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTac                     >  1:478737/1‑70 (MQ=255)
                                                ctctGATCATTTCAGCCATGTCTCGGTGATGCGCAAAGtgc                                                  >  1:2373363/1‑41 (MQ=255)
                                                 tCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                               <  1:696623/59‑1 (MQ=255)
                                                   cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGc                                 >  1:3770922/1‑55 (MQ=255)
                                                   cGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCtt                  <  1:1151396/70‑1 (MQ=255)
                                                     aCCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGa                                      >  1:5036744/1‑48 (MQ=255)
                                                      tcATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTAtc              <  1:3609129/70‑2 (MQ=255)
                                                      ccATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt              >  1:4091093/1‑71 (MQ=255)
                                                      ccATTTCAGCCATGTTTCGGGGATGCGGAAAGTGCTGCTAAGAATGACCAGc                                 <  1:2560631/52‑1 (MQ=255)
                                                       cATTTCAGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCa                                   >  1:887306/1‑49 (MQ=255)
                                                        aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                               <  1:258812/52‑1 (MQ=255)
                                                        aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg            <  1:1718755/71‑1 (MQ=255)
                                                        aTTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg            <  1:3148849/71‑1 (MQ=255)
                                                         tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGaa                                         >  1:4178161/1‑41 (MQ=255)
                                                         tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                             >  1:1763955/1‑53 (MQ=255)
                                                         tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGa                             >  1:2274517/1‑53 (MQ=255)
                                                         tttCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTg            <  1:4417888/70‑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:2616552/71‑1 (MQ=255)
                                                             aGCCATGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCa       <  1:1321618/71‑1 (MQ=255)
                                                               ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTAtt              <  1:2496535/62‑1 (MQ=255)
                                                               ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt     >  1:4362758/1‑71 (MQ=255)
                                                               ccATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAt      >  1:3999890/1‑70 (MQ=255)
                                                                cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                               <  1:3886935/44‑1 (MQ=255)
                                                                cATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTg    <  1:5332099/71‑1 (MQ=255)
                                                                 aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTg                               <  1:5676362/43‑1 (MQ=255)
                                                                 aTGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCAtt     >  1:603299/1‑69 (MQ=255)
                                                                  tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg           <  1:4604329/62‑1 (MQ=255)
                                                                  tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg           <  1:353479/62‑1 (MQ=255)
                                                                  tGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTgg           <  1:4432622/62‑1 (MQ=255)
                                                                  tGTCTCGGTGATGCGCAAAGTGCTATTACGCATGACCAGCTGGGTGCATTACCTTTATCTCGATCAttaa   <  1:5176130/70‑3 (MQ=255)
                                                                   gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                                    <  1:783149/36‑1 (MQ=255)
                                                                   gTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga   >  1:2894390/1‑69 (MQ=255)
                                                                   gTTTAGGTGATGCGGAAAGTGCTGCTAAGAATGAcc                                    <  1:2705178/36‑1 (MQ=255)
                                                                     ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga   >  1:665949/1‑67 (MQ=255)
                                                                     ttCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTga   >  1:3630509/1‑67 (MQ=255)
                                                                     |                                                                   
TGATATCGATGCCTTACCCGACGGCATCTACGCCCAGGTGGCGGTTTACTCCGACCATTTCAGCCATGTTTCGGTGATGCGGAAAGTGCTGCTAAGAATGACCAGCTGGATGCATTATCTTTATTTGGATCATTGAG  >  minE/2571097‑2571233

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: