Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,249,100 A→G 5.1% D97G (GAC→GGC)  thiF → thiamin (thiazole moiety) biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,249,1000AG 5.1% 286.9 / 5.3 98D97G (GAC→GGC) thiFthiamin (thiazole moiety) biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base A (47/46);  new base G (2/3);  total (49/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 102 total aligned reads displayed.

TTACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGCCGAT  >  minE/2249031‑2249170
                                                                     |                                                                      
ttACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGa                                                                        >  1:310457/1‑70 (MQ=255)
ttACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGa                                                                        >  1:3728011/1‑70 (MQ=255)
 tACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACa                                                                      >  1:1931028/1‑71 (MQ=255)
 tACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACa                                                                      >  1:1533246/1‑71 (MQ=255)
  aCCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACa                                                                      <  1:3486575/70‑1 (MQ=255)
  aCCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACAt                                                                     >  1:2225614/1‑71 (MQ=255)
   ccACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACa                                                                      <  1:377496/69‑1 (MQ=255)
         aaGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAAc                                                                <  1:2678762/69‑1 (MQ=255)
          aGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCaa                                                                 <  1:2955225/67‑1 (MQ=38)
          aGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCaa                                                                 <  1:1262891/67‑1 (MQ=38)
          aGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCaa                                                                 <  1:2776724/67‑1 (MQ=38)
             tatCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                          <  1:2109898/71‑1 (MQ=255)
             tatCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                          <  1:3253787/71‑1 (MQ=255)
             tatCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCAGACATTCAACTGACAg                                                          <  1:3601901/71‑1 (MQ=255)
              atcgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                          >  1:2002520/1‑70 (MQ=255)
               tcgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCa                                                                  <  1:1193030/61‑1 (MQ=255)
               tcgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCa                                                                  <  1:1612525/61‑1 (MQ=255)
                cgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCAt                                                       >  1:1006459/1‑71 (MQ=255)
                cgatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCAGACATTCAACTGACAGCAt                                                       >  1:1439610/1‑71 (MQ=255)
                 gatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                          <  1:511637/67‑1 (MQ=255)
                 gatcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCAtt                                                      <  1:2838668/71‑1 (MQ=255)
                  atcgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAg                                                          <  1:3313781/66‑1 (MQ=255)
                    cgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCAtt                                                      >  1:1704172/1‑68 (MQ=255)
                    cgCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCAt                                                       >  1:3391501/1‑67 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACAtt                                                                    >  1:1407061/1‑49 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTacaa                                                  >  1:2209311/1‑67 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTacaa                                                  >  1:1053957/1‑67 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACaacaac                                              >  1:2571391/1‑71 (MQ=255)
                         gAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACaacaac                                              >  1:59185/1‑71 (MQ=255)
                          aaaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACg                                             >  1:2066369/1‑71 (MQ=255)
                          aaaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACg                                             <  1:155787/71‑1 (MQ=255)
                          aaaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACg                                             <  1:2914453/71‑1 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                            >  1:2633378/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                            <  1:2202934/71‑1 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                            >  1:3552400/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                            >  1:3510594/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                            >  1:761964/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                            >  1:2568380/1‑71 (MQ=255)
                           aaTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                            >  1:2329542/1‑71 (MQ=255)
                            aTCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATGTCGGGATTCAACt                                                               >  1:471242/1‑51 (MQ=21)
                             tCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGtt                                          <  1:221539/71‑1 (MQ=255)
                             tCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGt                                           >  1:3226309/1‑70 (MQ=255)
                                cAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAAc                                       <  1:2550780/71‑1 (MQ=255)
                                 aGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACg                                      >  1:76835/1‑71 (MQ=255)
                                  ggTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACgg                                     >  1:3213508/1‑71 (MQ=255)
                                  ggTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACgg                                     >  1:1984309/1‑71 (MQ=255)
                                   gTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTAcaacaa                                               <  1:1062202/60‑1 (MQ=255)
                                   gTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTAcaacaa                                               <  1:2836354/60‑1 (MQ=255)
                                   gTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                            <  1:1235658/63‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACTGTTAACGGGt                                   <  1:1535744/71‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                   <  1:3089535/71‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                   <  1:259178/71‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                   <  1:2972876/71‑1 (MQ=255)
                                    tCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                   <  1:1666646/71‑1 (MQ=255)
                                      aGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGa                                 <  1:1009428/71‑1 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                                 >  1:3092991/1‑38 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                                 >  1:3606443/1‑38 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                                 >  1:1781061/1‑38 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                                 >  1:2267382/1‑38 (MQ=255)
                                       gCCAACAGCGACTGACACAGTTGAATCCCGACATTCaa                                                                 >  1:2376599/1‑38 (MQ=255)
                                        ccAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGa                                 >  1:746357/1‑69 (MQ=255)
                                             aGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGtt                           >  1:3680797/1‑70 (MQ=255)
                                             aGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGt                            >  1:2039817/1‑69 (MQ=255)
                                               cGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaa                        <  1:1750771/71‑1 (MQ=255)
                                               cGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaa                        <  1:1236793/71‑1 (MQ=255)
                                               cGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaa                        <  1:849242/71‑1 (MQ=255)
                                               cGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaa                         >  1:403044/1‑70 (MQ=255)
                                                gACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaa                        <  1:469085/70‑1 (MQ=255)
                                                 aCTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAg                      <  1:1206511/71‑1 (MQ=255)
                                                 aCTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAg                      <  1:2266739/71‑1 (MQ=255)
                                                   tGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                   <  1:893385/56‑1 (MQ=255)
                                                   tGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGt                                   <  1:2833306/56‑1 (MQ=255)
                                                    gACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGtt                                          <  1:2544291/48‑1 (MQ=255)
                                                      cacaGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCg                 <  1:3322120/71‑1 (MQ=255)
                                                      cacaGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCg                 >  1:2780328/1‑71 (MQ=255)
                                                      cacaGTTGAATCCCGACAGTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCg                 >  1:3370693/1‑71 (MQ=255)
                                                        caGTTGAATCCCGACATTCAACTGACAGCATTACAACAACgg                                            >  1:2991207/1‑42 (MQ=255)
                                                        caGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGt               >  1:1798135/1‑71 (MQ=255)
                                                        caGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGt               >  1:3383977/1‑71 (MQ=255)
                                                         aGTTGAATGTCGGGATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt              >  1:2863709/1‑71 (MQ=38)
                                                         aGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAg                      >  1:3271254/1‑63 (MQ=255)
                                                         aGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt              <  1:1157989/71‑1 (MQ=255)
                                                         aGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt              <  1:1282952/71‑1 (MQ=255)
                                                           ttGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGt                            >  1:2441946/1‑55 (MQ=255)
                                                           ttGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGt                            >  1:2437659/1‑55 (MQ=255)
                                                             gAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCa           <  1:1188539/70‑1 (MQ=255)
                                                             gAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCa           <  1:2431431/70‑1 (MQ=255)
                                                              aaTCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCAc          <  1:671049/70‑1 (MQ=255)
                                                                tCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTaaaa                       <  1:2404064/55‑1 (MQ=255)
                                                                tCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCgg                <  1:2973394/62‑1 (MQ=255)
                                                                 cccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACggg       >  1:2650473/1‑70 (MQ=255)
                                                                 cccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGc      >  1:3343496/1‑71 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt              >  1:374197/1‑62 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt              >  1:3087601/1‑62 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGtt              >  1:1560025/1‑62 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGcc     <  1:2241826/71‑1 (MQ=255)
                                                                  ccGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGcc     >  1:862905/1‑71 (MQ=255)
                                                                   cGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCa           <  1:302635/64‑1 (MQ=255)
                                                                     aCATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGCCGAt  >  1:2870498/1‑71 (MQ=255)
                                                                     aCATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGCCGAt  >  1:1423083/1‑71 (MQ=255)
                                                                     |                                                                      
TTACCACTGAAGATATCGATCGCCCGAAATCGCAGGTCAGCCAACAGCGACTGACACAGTTGAATCCCGACATTCAACTGACAGCATTACAACAACGGTTAACGGGTGAGGCGTTAAAAGATGCGGTTGCACGGGCCGAT  >  minE/2249031‑2249170

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: