Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,221,883 A→T 6.0% K153I (AAA→ATA)  yebT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,221,8830AT 6.0% 290.3 / 8.3 99K153I (AAA→ATA) yebTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (51/42);  new base T (2/4);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 104 total aligned reads displayed.

AACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGC  >  minE/1221813‑1221950
                                                                      |                                                                   
aaCTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGaa                                                                     <  1:3085703/71‑1 (MQ=255)
  cTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGaaa                                                                    >  1:4257227/1‑70 (MQ=255)
  cTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAat                                                                   <  1:2595770/71‑1 (MQ=255)
  cTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAat                                                                   <  1:3931379/71‑1 (MQ=255)
  cTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAat                                                                   <  1:351015/71‑1 (MQ=255)
     tatCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAatat                                                                 >  1:2924238/1‑70 (MQ=255)
     tatCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATc                                                                <  1:2798219/71‑1 (MQ=255)
     tatCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATc                                                                >  1:1528017/1‑71 (MQ=255)
          gCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                           >  1:70409/1‑71 (MQ=255)
            atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                         <  1:250044/71‑2 (MQ=38)
            atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                         <  1:3955505/71‑2 (MQ=38)
            atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                         <  1:2652691/71‑2 (MQ=38)
            atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                         <  1:2802799/71‑2 (MQ=38)
             tgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                         >  1:2012182/1‑70 (MQ=255)
              gatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTgggta                                                       >  1:3587952/1‑68 (MQ=37)
              gatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTgggta                                                       >  1:906849/1‑68 (MQ=37)
               atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                      >  1:2725512/1‑71 (MQ=255)
               atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                      >  1:3852980/1‑71 (MQ=255)
               atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                      >  1:802471/1‑71 (MQ=255)
               atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                      >  1:1029857/1‑71 (MQ=255)
                      gtaaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                               >  1:4177585/1‑71 (MQ=255)
                       taaaggtaaagAGCAGTATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                                 <  1:4320184/68‑1 (MQ=255)
                       taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAat                                                                   <  1:2924015/50‑1 (MQ=255)
                       taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                                 <  1:3559872/68‑1 (MQ=255)
                       taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                                 <  1:3412172/68‑1 (MQ=255)
                       taaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATc                                              >  1:2849333/1‑71 (MQ=255)
                         aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                            >  1:4073706/1‑55 (MQ=255)
                         aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                            >  1:3292777/1‑55 (MQ=255)
                         aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                            >  1:3794254/1‑55 (MQ=255)
                         aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCt                                             <  1:4350618/70‑1 (MQ=255)
                          aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                           >  1:3518189/1‑71 (MQ=255)
                          aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                           >  1:1707162/1‑71 (MQ=255)
                          aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                           >  1:3089539/1‑71 (MQ=255)
                          aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                           >  1:3657989/1‑71 (MQ=255)
                           ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                         >  1:2567144/1‑56 (MQ=255)
                           ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                         >  1:1809656/1‑56 (MQ=255)
                           ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                          <  1:3418502/71‑1 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGaaa                                                                    >  1:440773/1‑44 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGaaa                                                                    >  1:1159695/1‑44 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGaaa                                                                    >  1:2771463/1‑44 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                           <  1:4046307/53‑1 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                           <  1:859358/53‑1 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                           <  1:2656646/53‑1 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                      >  1:2414956/1‑58 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgatg                                         <  1:452570/71‑1 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgatg                                         <  1:1226216/71‑1 (MQ=255)
                            gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                          >  1:4003162/1‑70 (MQ=255)
                             taaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGAGATCtgatg                                         >  1:2503251/1‑70 (MQ=255)
                              aaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                               >  1:3519554/1‑63 (MQ=255)
                               aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                               >  1:3959734/1‑62 (MQ=255)
                               aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                               >  1:3108247/1‑62 (MQ=255)
                               aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATc                                      >  1:3429980/1‑71 (MQ=255)
                                agagCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCg                                                               <  1:2750500/45‑1 (MQ=255)
                                agagCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCg                                                               <  1:1423413/45‑1 (MQ=255)
                                agagCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCg                                                               <  1:868985/45‑1 (MQ=255)
                                   gCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCAc                                   <  1:2187188/70‑1 (MQ=255)
                                    cAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCt                                 <  1:1151494/71‑1 (MQ=255)
                                    cAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCt                                 <  1:1546748/71‑1 (MQ=255)
                                     aGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTg                                >  1:3115701/1‑71 (MQ=255)
                                      ggATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGc                               <  1:3742953/71‑1 (MQ=255)
                                        aTCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGc                               <  1:2464709/69‑1 (MQ=255)
                                             tttGTCGCACTCGATACCCAACCGAAATATCGGCTGGAc                                                        <  1:1780079/39‑1 (MQ=255)
                                             tttGTCGCACTCGATACCCAACCGAAATATCGGCTGGAc                                                        <  1:806404/39‑1 (MQ=255)
                                               tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                      >  1:1275829/1‑39 (MQ=255)
                                               tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                      >  1:3476280/1‑39 (MQ=255)
                                               tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                               >  1:1969534/1‑46 (MQ=255)
                                               tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGccc                         <  1:1301739/68‑1 (MQ=255)
                                                   gCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtct                  <  1:1955496/71‑1 (MQ=255)
                                                    cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc                 >  1:3299012/1‑71 (MQ=255)
                                                    cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc                 >  1:191557/1‑71 (MQ=255)
                                                    cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc                 <  1:1100457/71‑1 (MQ=255)
                                                      cTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCgg               >  1:2097516/1‑71 (MQ=255)
                                                      cTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCgg               <  1:2811734/71‑1 (MQ=255)
                                                      cTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCgg               >  1:77109/1‑71 (MQ=255)
                                                       tCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                          <  1:3575275/43‑1 (MQ=255)
                                                       tCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                          <  1:3488793/43‑1 (MQ=255)
                                                        cGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCAc                                   >  1:2631649/1‑49 (MQ=255)
                                                        cGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGtt             >  1:2608283/1‑71 (MQ=255)
                                                         gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAg                            >  1:2740496/1‑55 (MQ=255)
                                                         gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc            >  1:1940407/1‑71 (MQ=255)
                                                         gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc            >  1:4136526/1‑71 (MQ=255)
                                                         gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc            >  1:1348553/1‑71 (MQ=255)
                                                           tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                       <  1:3089195/42‑1 (MQ=255)
                                                           tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                       <  1:1265361/42‑1 (MQ=255)
                                                           tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                       <  1:2289012/42‑1 (MQ=255)
                                                           tACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCa                              >  1:1649786/1‑51 (MQ=255)
                                                           tACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCa                              >  1:741907/1‑51 (MQ=255)
                                                             cccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCt         >  1:614977/1‑70 (MQ=255)
                                                             cccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTg        <  1:4404829/71‑1 (MQ=255)
                                                              ccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtct                  <  1:2667694/60‑1 (MQ=255)
                                                              ccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtct                  <  1:1911608/60‑1 (MQ=255)
                                                               cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa      >  1:1774182/1‑71 (MQ=255)
                                                               cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa      >  1:2272879/1‑71 (MQ=255)
                                                               cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa      >  1:2990037/1‑71 (MQ=255)
                                                               cAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa      >  1:3159878/1‑71 (MQ=255)
                                                                aaCCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAAc     <  1:2911884/71‑1 (MQ=255)
                                                                 aCCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa      >  1:3417214/1‑69 (MQ=255)
                                                                 aCCGAAATATAGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGaa      >  1:1017712/1‑69 (MQ=255)
                                                                  ccGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAg   <  1:814123/71‑1 (MQ=255)
                                                                   cGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGc  <  1:495221/71‑1 (MQ=255)
                                                                      |                                                                   
AACTATATCGGCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGC  >  minE/1221813‑1221950

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: