Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,570,756 G→A 16.2% V221I (GTC→ATC) ‡ yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,570,7560GA16.2% 275.4 / 37.9 117V221I (GTC→ATC) ‡yibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/54);  new base A (12/7);  total (56/61)
Fisher's exact test for biased strand distribution p-value = 2.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 119 total aligned reads displayed.

CGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTG  >  minE/2570686‑2570822
                                                                      |                                                                  
cGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTg                                                                    <  1:2334190/71‑1 (MQ=255)
cGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTg                                                                    <  1:4627215/71‑1 (MQ=255)
 ggCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGt                                                                   >  1:2578877/1‑71 (MQ=255)
 ggCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGt                                                                   >  1:1530045/1‑71 (MQ=255)
      gCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTc                                                              <  1:2452704/71‑1 (MQ=255)
      gCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTc                                                              <  1:3401954/71‑1 (MQ=255)
          tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                          >  1:5307078/1‑71 (MQ=255)
          tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                          >  1:19330/1‑71 (MQ=255)
          tCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgcgc                                                          >  1:508860/1‑71 (MQ=255)
            gCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCgcgcgc                                                        >  1:3657184/1‑71 (MQ=255)
                aTTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTcgc                                                            <  1:1117837/63‑1 (MQ=255)
                aTTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAg                                                    >  1:4646603/1‑71 (MQ=255)
                 ttGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   <  1:3267752/71‑1 (MQ=255)
                  tgtgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   >  1:3519118/1‑70 (MQ=255)
                    tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGATCAGACTGTCATTCGCGCGCCGAGCaa                                                 <  1:4146578/70‑1 (MQ=255)
                    tgAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                 <  1:2134918/70‑1 (MQ=255)
                      aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   <  1:3576077/66‑1 (MQ=255)
                      aGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATg                                               <  1:5068480/70‑1 (MQ=255)
                       gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                             <  1:4609709/71‑1 (MQ=255)
                       gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                             <  1:4739520/71‑1 (MQ=255)
                       gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                             <  1:2443825/71‑1 (MQ=255)
                       gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                             <  1:249362/71‑1 (MQ=255)
                       gCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                             <  1:1826526/71‑1 (MQ=255)
                         ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGt                                                                   <  1:2248416/47‑1 (MQ=255)
                         ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAAt                                                >  1:750550/1‑66 (MQ=255)
                         ttAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTa                                           >  1:2094604/1‑71 (MQ=255)
                          tAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTAc                                          >  1:2671339/1‑71 (MQ=255)
                           aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                 >  1:3434175/1‑47 (MQ=255)
                           aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCa                                                                 >  1:4540073/1‑47 (MQ=255)
                           aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                         >  1:3452316/1‑71 (MQ=255)
                           aaGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                         >  1:5260611/1‑71 (MQ=255)
                           aaGAGCGCAACTTACTGAAGCAAAATAAAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACg                                         >  1:4828358/1‑71 (MQ=255)
                            agagCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGt                                        <  1:5429040/71‑1 (MQ=255)
                              agCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTa                                      >  1:4617417/1‑71 (MQ=255)
                               gcgcAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTAc                                     >  1:217300/1‑71 (MQ=255)
                                  cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   <  1:989349/54‑1 (MQ=255)
                                  cAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACt                                    >  1:2620716/1‑69 (MQ=255)
                                   aaCTTACTGAAGCAAAATATAACCTTGAGTAGACTGTCATTCGCGCGCCGAGCAATGGCt                                            <  1:1854577/60‑1 (MQ=255)
                                   aaCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                            <  1:2430738/60‑1 (MQ=255)
                                   aaCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCt                                            <  1:5539252/60‑1 (MQ=255)
                                    aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCa                                                  <  1:2932552/53‑1 (MQ=255)
                                    aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgt                                <  1:1579372/71‑2 (MQ=255)
                                    aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                >  1:766332/1‑71 (MQ=255)
                                    aCTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                >  1:5504789/1‑71 (MQ=255)
                                     cTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                 >  1:2418141/1‑53 (MQ=255)
                                      ttACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   <  1:2555669/50‑1 (MQ=255)
                                        aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                 >  1:2963613/1‑50 (MQ=255)
                                        aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCaa                                                 >  1:3004292/1‑50 (MQ=255)
                                        aCTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCATGTa                              <  1:3341062/69‑1 (MQ=255)
                                         cTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTa                              <  1:1242373/68‑1 (MQ=255)
                                          tGAAGCAAAATATAACCTTGAGCAGACTGTCATTcg                                                             <  1:4180285/36‑1 (MQ=255)
                                           gAAGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATggtt                                            <  1:4658593/52‑3 (MQ=255)
                                            aaGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATc                        >  1:5325897/1‑71 (MQ=255)
                                            aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                        >  1:5628654/1‑70 (MQ=255)
                                            aaGCAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGAtt                        >  1:2660324/1‑70 (MQ=255)
                                              gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                       <  1:617446/38‑1 (MQ=255)
                                              gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                       <  1:4456269/38‑1 (MQ=255)
                                              gCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGcc                                                       <  1:1806351/38‑1 (MQ=255)
                                               cAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGc                                             <  1:5421041/47‑1 (MQ=255)
                                               cAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCg                      <  1:930966/70‑1 (MQ=255)
                                                aaaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGAc                                     >  1:4144447/1‑54 (MQ=255)
                                                aaaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgt                     <  1:212453/70‑2 (MQ=255)
                                                 aaaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGtt                                       <  1:3532462/51‑1 (MQ=255)
                                                 aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGAt                                                <  1:1039785/42‑1 (MQ=12)
                                                 aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGAt                                                <  1:2680923/42‑1 (MQ=12)
                                                 aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                   <  1:3208415/71‑1 (MQ=255)
                                                 aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                   >  1:1253642/1‑71 (MQ=255)
                                                  aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   <  1:3725302/38‑1 (MQ=255)
                                                  aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   <  1:3718608/38‑1 (MQ=255)
                                                  aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCa                  <  1:1698509/71‑1 (MQ=255)
                                                  aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTa                              >  1:3456040/1‑59 (MQ=255)
                                                  aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                  >  1:2733471/1‑70 (MQ=255)
                                                  aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                  >  1:2880584/1‑70 (MQ=255)
                                                  aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                  >  1:2610651/1‑70 (MQ=255)
                                                   atatAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                 >  1:1906756/1‑71 (MQ=255)
                                                    tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGc                                                   <  1:3082228/36‑1 (MQ=255)
                                                    tataACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAgg                >  1:3649775/1‑71 (MQ=255)
                                                      taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa              >  1:1397815/1‑71 (MQ=255)
                                                      taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa              >  1:4119278/1‑71 (MQ=255)
                                                      taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa              <  1:3749721/71‑1 (MQ=255)
                                                      taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTa              >  1:5628176/1‑71 (MQ=255)
                                                        aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa            <  1:2085988/71‑1 (MQ=255)
                                                         ccTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgtc                    >  1:3467639/1‑60 (MQ=255)
                                                          cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                <  1:2274309/49‑1 (MQ=255)
                                                          cTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAgg                                <  1:2222215/49‑1 (MQ=255)
                                                          cTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTa          >  1:1271195/1‑71 (MQ=255)
                                                           ttGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGa                          >  1:562770/1‑54 (MQ=255)
                                                            tGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACg        >  1:2174325/1‑71 (MQ=255)
                                                             gAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGa                          >  1:2387129/1‑52 (MQ=255)
                                                              aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGc       <  1:1028155/70‑1 (MQ=255)
                                                              aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGc       <  1:2295133/70‑1 (MQ=255)
                                                                 aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt   <  1:3218372/71‑1 (MQ=255)
                                                                 aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt   <  1:3225794/71‑1 (MQ=255)
                                                                 aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt   <  1:428766/71‑1 (MQ=255)
                                                                  gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGc                     >  1:359309/1‑52 (MQ=255)
                                                                  gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa            <  1:1560887/61‑1 (MQ=255)
                                                                  gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg  <  1:2670348/71‑1 (MQ=255)
                                                                  gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg  >  1:3098634/1‑71 (MQ=255)
                                                                    cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc    >  1:917187/1‑67 (MQ=255)
                                                                    cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc    >  1:3899776/1‑67 (MQ=255)
                                                                      |                                                                  
CGGCGAGCAATCGCAGATTGTGAGCTTAAGAGCGCAACTTACTGAAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTG  >  minE/2570686‑2570822

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: