Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,570,800 C→T 17.8% I235I (ATC→ATT yibH → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,570,8000CT17.8% 247.4 / 41.2 112I235I (ATC→ATTyibHhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (45/47);  new base T (14/6);  total (59/53)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

Only 100 of 113 total aligned reads displayed.

AAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGG  >  minE/2570730‑2570869
                                                                      |                                                                     
   cAAAATATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCg                                                                     <  1:930966/70‑1 (MQ=255)
    aaaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgt                                                                    <  1:212453/70‑2 (MQ=255)
     aaaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtcc                                                                  >  1:1253642/1‑71 (MQ=255)
      aaTATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCa                                                                 <  1:1698509/71‑1 (MQ=255)
      aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                                                                 >  1:2880584/1‑70 (MQ=255)
      aaTATAACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGtccg                                                                 >  1:2610651/1‑70 (MQ=255)
       atatAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAg                                                                >  1:1906756/1‑71 (MQ=255)
          taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                                             >  1:1397815/1‑71 (MQ=255)
          taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                                             <  1:3749721/71‑1 (MQ=255)
          taACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTa                                                             >  1:4119278/1‑71 (MQ=255)
          taACCTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTa                                                             >  1:5628176/1‑71 (MQ=255)
            aCCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                                                           <  1:2085988/71‑1 (MQ=255)
             ccTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCgtc                                                                   >  1:3467639/1‑60 (MQ=255)
              cTTGAGCAGACGATAGTGCGTGCGCCGAGCGATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTa                                                         >  1:1271195/1‑71 (MQ=255)
                tGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACg                                                       >  1:2174325/1‑71 (MQ=255)
                  aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGc                                                      <  1:2295133/70‑1 (MQ=255)
                  aGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGc                                                      <  1:1028155/70‑1 (MQ=255)
                     aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt                                                  <  1:428766/71‑1 (MQ=255)
                     aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt                                                  <  1:3218372/71‑1 (MQ=255)
                     aGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCt                                                  <  1:3225794/71‑1 (MQ=255)
                      gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGc                                                                    >  1:359309/1‑52 (MQ=255)
                      gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACa                                                           <  1:1560887/61‑1 (MQ=255)
                      gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                                 >  1:3098634/1‑71 (MQ=255)
                      gACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                                 <  1:2670348/71‑1 (MQ=255)
                        cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                                   >  1:3899776/1‑67 (MQ=255)
                        cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                                   >  1:5506930/1‑67 (MQ=255)
                        cTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGc                                                   >  1:917187/1‑67 (MQ=255)
                           tCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCtt                                             <  1:124460/70‑1 (MQ=255)
                             aTTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc                                          >  1:2726988/1‑71 (MQ=255)
                              ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         <  1:3667828/71‑1 (MQ=255)
                              ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         <  1:4652072/71‑1 (MQ=255)
                              ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         <  1:3303062/71‑1 (MQ=255)
                              ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         <  1:2016074/71‑1 (MQ=255)
                              ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         <  1:3065483/71‑1 (MQ=255)
                              ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         <  1:1277911/71‑1 (MQ=255)
                              ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         <  1:5155473/71‑1 (MQ=255)
                              ttCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         <  1:96965/71‑1 (MQ=255)
                              ttCGCGCGCCGAGCAAAGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         <  1:1593534/71‑1 (MQ=255)
                               tCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                                        >  1:4405786/1‑71 (MQ=255)
                                cgcgcgCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCg                                         >  1:4819656/1‑69 (MQ=255)
                                   gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                                            >  1:2616414/1‑47 (MQ=255)
                                   gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                                            >  1:2029958/1‑47 (MQ=255)
                                   gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                                            >  1:3886713/1‑47 (MQ=255)
                                   gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                                            >  1:974992/1‑47 (MQ=255)
                                   gcgcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTAc                                                            >  1:3157193/1‑47 (MQ=255)
                                    cgcCGAGCCATGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTatg                                                       >  1:5221741/1‑49 (MQ=255)
                                     gcagAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                  >  1:5095792/4‑71 (MQ=255)
                                     gcCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                  >  1:5027913/1‑71 (MQ=255)
                                        gAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                               <  1:4523482/71‑1 (MQ=255)
                                        gAGCAATGGCTACGTTACTCAGGTACTGAACCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTcc                                 >  1:2462052/1‑69 (MQ=255)
                                         aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                                 >  1:13360/1‑52 (MQ=255)
                                         aGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                                 >  1:1658442/1‑52 (MQ=255)
                                          gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                               <  1:3771478/69‑1 (MQ=255)
                                          gCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCgg                               <  1:2791965/69‑1 (MQ=255)
                                             aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg                          <  1:4950043/71‑1 (MQ=255)
                                             aTGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATg                          >  1:544441/1‑71 (MQ=255)
                                             aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                  <  1:2249279/63‑1 (MQ=255)
                                             aTGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                  <  1:5461260/63‑1 (MQ=255)
                                              tGGTTATGTGACCCAGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATgg                         >  1:622311/1‑71 (MQ=255)
                                              tGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGa                            <  1:198856/68‑1 (MQ=255)
                                               ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                                >  1:4733665/1‑63 (MQ=255)
                                               ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCg                                >  1:951436/1‑63 (MQ=255)
                                               ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt                        >  1:4296935/1‑71 (MQ=255)
                                               ggCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGt                        >  1:5129572/1‑71 (MQ=255)
                                                gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                       <  1:2328269/71‑1 (MQ=255)
                                                gCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                       >  1:3392954/1‑71 (MQ=255)
                                                 cTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCt                      <  1:4172809/71‑1 (MQ=255)
                                                    cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTc                    <  1:4657486/70‑1 (MQ=255)
                                                    cGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCa                   >  1:2203160/1‑71 (MQ=255)
                                                     gTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGcc                                          <  1:4329799/47‑1 (MQ=255)
                                                      ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                                        <  1:3557122/48‑1 (MQ=255)
                                                      ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTc                                  <  1:2937265/54‑1 (MQ=255)
                                                      ttACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATc                 >  1:508147/1‑71 (MQ=255)
                                                        aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTg                                      >  1:2601556/1‑48 (MQ=255)
                                                        aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTg                                      >  1:2121406/1‑48 (MQ=255)
                                                        aCTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATcc                >  1:415778/1‑70 (MQ=255)
                                                         cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                                 <  1:3699059/36‑1 (MQ=255)
                                                         cTCAGGTACTGATCCGCCCAGGTACATACGCAGCTg                                                 <  1:879172/36‑1 (MQ=255)
                                                            aGGTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCg                                >  1:126606/1‑50 (MQ=255)
                                                            aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGc                                        >  1:4093893/1‑42 (MQ=255)
                                                            aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTc                       >  1:4223055/1‑59 (MQ=255)
                                                            aGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg           >  1:2878984/1‑71 (MQ=255)
                                                             ggTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg           >  1:976414/1‑70 (MQ=255)
                                                             ggTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAg           >  1:3890918/1‑70 (MQ=255)
                                                              gTACTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCa         <  1:2235105/71‑1 (MQ=255)
                                                              gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                     <  1:554302/59‑1 (MQ=255)
                                                              gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                     <  1:4955373/59‑1 (MQ=255)
                                                              gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                     <  1:2622058/59‑1 (MQ=255)
                                                              gTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCtt                     <  1:261161/59‑1 (MQ=255)
                                                                aCTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCag        >  1:5457882/1‑69 (MQ=255)
                                                                 cTGATTCGTCCGGGTACCTATGCCGCGTCGCTGCCGCTACGTCCGGTGATGGTGTTTATACCCGATCaga       <  1:1211790/70‑3 (MQ=255)
                                                                 cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaa      <  1:3691402/71‑1 (MQ=255)
                                                                 cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa       <  1:3172040/70‑1 (MQ=255)
                                                                 cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa       <  1:2417283/70‑1 (MQ=255)
                                                                 cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa       <  1:1625402/70‑1 (MQ=255)
                                                                 cTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaa       <  1:1549848/70‑1 (MQ=255)
                                                                   gATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCa         >  1:1108039/1‑66 (MQ=255)
                                                                    aTCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCaaaaa     <  1:5462102/69‑1 (MQ=255)
                                                                     tCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCAt                  >  1:2717923/1‑55 (MQ=255)
                                                                     tCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACgg  >  1:4744885/1‑71 (MQ=255)
                                                                      |                                                                     
AAGCAAAATATAACCTTGAGCAGACTGTCATTCGCGCGCCGAGCAATGGCTACGTTACTCAGGTACTGATCCGCCCAGGTACATACGCAGCTGCCTTGCCGCTGCGTCCGGTGATGGTCTTCATCCCCGAGCAAAAACGG  >  minE/2570730‑2570869

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: