Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,221,891 C→T 6.1% L156L (CTG→TTG)  yebT → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,221,8910CT 6.1% 333.4 / 7.6 98L156L (CTG→TTG) yebTconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (55/37);  new base T (2/4);  total (57/41)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

GCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCATTGGT  >  minE/1221823‑1221961
                                                                    |                                                                      
gCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                                      >  1:70409/1‑71 (MQ=255)
  atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                                    <  1:250044/71‑2 (MQ=38)
  atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                                    <  1:2652691/71‑2 (MQ=38)
  atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                                    <  1:2802799/71‑2 (MQ=38)
  atgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTggg                                                                    <  1:3955505/71‑2 (MQ=38)
   tgatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                                    >  1:2012182/1‑70 (MQ=255)
    gatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTgggta                                                                  >  1:906849/1‑68 (MQ=37)
    gatgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAGCCGATATTTCGGTTgggta                                                                  >  1:3587952/1‑68 (MQ=37)
     atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                                 >  1:1029857/1‑71 (MQ=255)
     atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                                 >  1:2725512/1‑71 (MQ=255)
     atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                                 >  1:3852980/1‑71 (MQ=255)
     atgCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                                 >  1:802471/1‑71 (MQ=255)
            gtaaaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                                          >  1:4177585/1‑71 (MQ=255)
             taaaggtaaagAGCAGTATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCg                                                            <  1:4320184/68‑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:3794254/1‑55 (MQ=255)
               aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                                       >  1:4073706/1‑55 (MQ=255)
               aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCt                                                                       >  1:3292777/1‑55 (MQ=255)
               aaggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCt                                                        <  1:4350618/70‑1 (MQ=255)
                aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                                      >  1:3657989/1‑71 (MQ=255)
                aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                                      >  1:3518189/1‑71 (MQ=255)
                aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                                      >  1:3089539/1‑71 (MQ=255)
                aggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtga                                                      >  1:1707162/1‑71 (MQ=255)
                 ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                                    >  1:1809656/1‑56 (MQ=255)
                 ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGa                                                                    >  1:2567144/1‑56 (MQ=255)
                 ggtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                                     <  1:3418502/71‑1 (MQ=255)
                  gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                                      <  1:859358/53‑1 (MQ=255)
                  gtaaagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTg                                                                      <  1:4046307/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:3108247/1‑62 (MQ=255)
                     aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGAt                                                          >  1:3959734/1‑62 (MQ=255)
                     aagAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATc                                                 >  1:3429980/1‑71 (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:3476280/1‑39 (MQ=255)
                                     tGTCGCACTCGATACCCAACCGAAATATCGGCTGGACaa                                                                 >  1:1275829/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:191557/1‑71 (MQ=255)
                                          cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc                            <  1:1100457/71‑1 (MQ=255)
                                          cACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc                            >  1:3299012/1‑71 (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:3488793/43‑1 (MQ=255)
                                             tCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCtgat                                                     <  1:3575275/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:1348553/1‑71 (MQ=255)
                                               gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc                       >  1:3863571/1‑71 (MQ=255)
                                               gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc                       <  1:2952791/71‑1 (MQ=255)
                                               gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc                       >  1:4136526/1‑71 (MQ=255)
                                               gATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTc                       >  1:1940407/1‑71 (MQ=255)
                                                 tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                                  <  1:2289012/42‑1 (MQ=255)
                                                 tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                                  <  1:3089195/42‑1 (MQ=255)
                                                 tACCCAACCGAAATATCGGCTGGACAATGGCGATCTgatgat                                                  <  1:1265361/42‑1 (MQ=255)
                                                 tACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCa                                         >  1:741907/1‑51 (MQ=255)
                                                 tACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCa                                         >  1:1649786/1‑51 (MQ=255)
                                                   cccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCt                    >  1:614977/1‑70 (MQ=255)
                                                   cccAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTg                   <  1:979830/71‑1 (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)
                                                              tatCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGTTTCGCTGAACAGCGGTTc        <  1:2176364/71‑1 (MQ=255)
                                                              tatCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGtt         >  1:3391205/1‑70 (MQ=255)
                                                               atCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCa       >  1:810067/1‑71 (MQ=255)
                                                                tCGGCTGGACAATGGCTATCTGATGATCCACCTGCAAGCCCCCGAGCTCGGTTCGCTGAACAGCGGTTc        <  1:347844/69‑1 (MQ=255)
                                                                tCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGAtctc                            >  1:512257/1‑49 (MQ=255)
                                                                tCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCAt      <  1:174165/71‑1 (MQ=255)
                                                                tCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCAt      >  1:2035508/1‑71 (MQ=255)
                                                                tCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCAt      >  1:2444025/1‑71 (MQ=255)
                                                                tCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCAt      >  1:2058439/1‑71 (MQ=255)
                                                                tCGGCTGGACAATGGCGATCTGATCATCCACCTGCAAGCCCCCGATCTCGGTTCGCTg                   >  1:4093425/1‑58 (MQ=255)
                                                                    cTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCATTggg  <  1:779825/71‑2 (MQ=255)
                                                                    cTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCATTGGt  <  1:1293498/71‑1 (MQ=255)
                                                                    |                                                                      
GCATGATGCCGGGTAAAGGTAAAGAGCAGGATCACTTTGTCGCACTCGATACCCAACCGAAATATCGGCTGGACAATGGCGATCTGATGATCCACCTGCAAGCCCCCGATCTCGGTTCGCTGAACAGCGGTTCATTGGT  >  minE/1221823‑1221961

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: