Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 539,433 T→C 7.5% V132V (GTA→GTG dps ← Fe‑binding and storage protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE539,4330TC 7.5% 276.9 / 11.1 93V132V (GTA→GTGdpsFe‑binding and storage protein
Reads supporting (aligned to +/- strand):  ref base T (45/41);  new base C (2/5);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 2.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

CAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTTGTGGATGT  >  minE/539364‑539500
                                                                     |                                                                   
cAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGt                                                                     <  1:3318879/70‑1 (MQ=255)
cAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGt                                                                     <  1:1756777/70‑1 (MQ=255)
 aGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTAc                                                                   <  1:1620344/71‑1 (MQ=255)
   gTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGt                                                                 <  1:2555747/71‑1 (MQ=255)
   gTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGt                                                                 >  1:2034545/1‑71 (MQ=255)
   gTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGt                                                                 <  1:2577041/71‑1 (MQ=255)
     cgcgAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCa                                                               >  1:63409/1‑71 (MQ=255)
     cgcgAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCa                                                               >  1:2101181/1‑71 (MQ=255)
      gcgAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCa                                                               <  1:4476042/70‑1 (MQ=255)
      gcgAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCa                                                               <  1:2196974/70‑1 (MQ=255)
      gcgAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCa                                                               <  1:2044145/70‑1 (MQ=255)
       cgAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCAtt                                                             <  1:4674261/71‑1 (MQ=255)
       cgAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCAtt                                                             <  1:3606796/71‑1 (MQ=255)
        gagaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCAt                                                              >  1:2082601/1‑69 (MQ=255)
        gagaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGCACGTCAtt                                                             <  1:373994/70‑1 (MQ=255)
        gagaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGCACGTCAtt                                                             <  1:1691773/70‑1 (MQ=255)
        gagaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGCACGTCAtt                                                             <  1:1218232/70‑1 (MQ=255)
        gagaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGCACGTCAtt                                                             <  1:4148569/70‑1 (MQ=255)
        gagaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGCACGTCAtt                                                             <  1:551324/70‑1 (MQ=255)
          gaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGc                                                          >  1:1655016/1‑71 (MQ=255)
          gaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGc                                                          >  1:3175187/1‑71 (MQ=255)
          gaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGCACGTCATTAg                                                           >  1:2520313/1‑70 (MQ=255)
          gaCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTAGCCAATCGCTTTGCGCACGTCATTAg                                                           >  1:4373151/1‑70 (MQ=255)
           aCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGc                                                          >  1:3532098/1‑70 (MQ=255)
           aCGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGc                                                          >  1:270191/1‑70 (MQ=255)
             gcggcgGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGcgac                                                       >  1:3613412/1‑71 (MQ=255)
              cggcggTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGcgacg                                                      <  1:699837/71‑1 (MQ=255)
              cggcggTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGcgacg                                                      <  1:223168/71‑1 (MQ=255)
                   gTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGcgacg                                                      <  1:1012496/66‑1 (MQ=255)
                    tCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCg                                                >  1:3373484/1‑71 (MQ=255)
                     cAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTa                                                            >  1:1109928/1‑58 (MQ=255)
                       ggATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTaa                                             <  1:625700/71‑1 (MQ=255)
                        gATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAAc                                            >  1:3421991/1‑71 (MQ=255)
                        gATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAAc                                            >  1:2541540/1‑71 (MQ=255)
                          tatCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACgg                                          >  1:1012183/1‑71 (MQ=255)
                              tGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGcgac                                                       >  1:3114703/1‑54 (MQ=255)
                              tGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGcgac                                                       >  1:959799/1‑54 (MQ=255)
                                cGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGc                                     >  1:1853102/1‑70 (MQ=255)
                                     tcgtcATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCCCCGATTGCGt                                               >  1:3799222/1‑55 (MQ=255)
                                       gtcATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTc                               <  1:1221106/69‑1 (MQ=255)
                                       gtcATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTc                               <  1:3872477/69‑1 (MQ=255)
                                       gtcATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTc                               <  1:269894/69‑1 (MQ=255)
                                       gtcATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCtt                             >  1:235059/1‑71 (MQ=255)
                                       gtcATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCt                              <  1:4295028/70‑1 (MQ=255)
                                       gtcATCTTTCGCTGCGCCAATCGCTTTGCGTACGTCATTAg                                                           <  1:4005168/41‑1 (MQ=255)
                                        tcATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCttt                            <  1:1066643/71‑1 (MQ=255)
                                        tcATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCttt                            <  1:1507555/71‑1 (MQ=255)
                                        tcATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCtgt                            <  1:3211853/71‑3 (MQ=255)
                                          aTCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCttt                            >  1:854449/1‑69 (MQ=255)
                                          aTCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCttt                            >  1:4573221/1‑69 (MQ=255)
                                          aTCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCt                              >  1:2227540/1‑67 (MQ=255)
                                           tCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGt                                 <  1:2344504/63‑1 (MQ=255)
                                           tCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGt                                 <  1:2151191/63‑1 (MQ=255)
                                           tCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGt                                 <  1:728736/63‑1 (MQ=255)
                                           tCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGt                                 <  1:1701597/63‑1 (MQ=255)
                                           tCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGt                                 <  1:3955287/63‑1 (MQ=255)
                                           tCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGt                                 <  1:2209235/63‑1 (MQ=255)
                                           tCTTTCGCTTCGCCAATCGATTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGt                                 <  1:151320/63‑1 (MQ=255)
                                            cTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGcgac                                                       >  1:4467551/1‑40 (MQ=255)
                                            cTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGcgac                                                       >  1:1142341/1‑40 (MQ=255)
                                            cTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTc                                        >  1:3736170/1‑55 (MQ=255)
                                            cTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTc                                        >  1:39396/1‑55 (MQ=255)
                                            cTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTc                                        >  1:2461658/1‑55 (MQ=255)
                                            cTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTc                                        >  1:3437394/1‑55 (MQ=255)
                                            cTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTc                                        >  1:743323/1‑55 (MQ=255)
                                             tttCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGt                       >  1:4187738/1‑71 (MQ=255)
                                             tttCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGt                       >  1:227971/1‑71 (MQ=255)
                                              ttcgcttcgcCAATCGCTTTGCGTACGTCATTAGCGACGAt                                                    >  1:2681437/1‑41 (MQ=255)
                                              ttcgcttcgcCAATCGCTTTGCGTACGTCATTAGCGACGAt                                                    >  1:4552551/1‑41 (MQ=255)
                                                cgcttcgcCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGa                     <  1:4119343/70‑1 (MQ=255)
                                                 gcttcgcCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGAt                    <  1:2063998/70‑1 (MQ=255)
                                                   ttcgcCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCt                 >  1:1221642/1‑71 (MQ=255)
                                                    tcgcCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCt                 >  1:1483666/1‑70 (MQ=255)
                                                     cgcCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTg                >  1:3696813/1‑70 (MQ=255)
                                                     cgcCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTg                >  1:1257251/1‑70 (MQ=255)
                                                     cgcCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTg                >  1:2548571/1‑70 (MQ=255)
                                                     cgcCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTg                >  1:3387320/1‑70 (MQ=255)
                                                         aaTCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGt           <  1:1518852/71‑1 (MQ=255)
                                                         aaTCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCACTTCTTTCAGGTGATCCTGAACGt           <  1:2898700/71‑1 (MQ=255)
                                                          aTCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGtt          >  1:185442/1‑71 (MQ=255)
                                                          aTCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGt           <  1:1998694/70‑1 (MQ=255)
                                                          aTCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGt           >  1:4558850/1‑70 (MQ=255)
                                                           tCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTtg         <  1:2552308/71‑1 (MQ=255)
                                                           tCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTtg         <  1:1529562/71‑1 (MQ=255)
                                                           tCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTtg         <  1:3723952/71‑1 (MQ=255)
                                                           tCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTtg         <  1:3158847/71‑1 (MQ=255)
                                                             gCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTtgtg       >  1:2328301/1‑71 (MQ=255)
                                                             gCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTtgtg       >  1:4365133/1‑71 (MQ=255)
                                                             gCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTtgtg       >  1:1945514/1‑71 (MQ=255)
                                                             gCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTtgtg       >  1:219107/1‑71 (MQ=255)
                                                              cTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTTGTgg      <  1:1677097/71‑1 (MQ=255)
                                                              cTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTTGTgg      <  1:3504009/71‑1 (MQ=255)
                                                              cTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTTGTgg      <  1:3442409/71‑1 (MQ=255)
                                                              cTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTTGTgg      <  1:3019641/71‑1 (MQ=255)
                                                                   cGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTTGTGGATGt  <  1:2661567/70‑1 (MQ=255)
                                                                     |                                                                   
CAGGTCGCGAGACGCGGCGGTCAGGATATCTGCGGTGTCGTCATCTTTCGCTTCGCCAATCGCTTTGCGTACGTCATTAGCGACGATTGCGTAACGGTCAGCCAGTTCTTTCAGGTGATCCTGAACGTTGTGGATGT  >  minE/539364‑539500

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: