Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,104,147 C→G 17.1% T49T (ACC→ACG kbaZ → tagatose 6‑phosphate aldolase 1, kbaZ subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,104,1470CG17.1% 284.2 / 45.5 123T49T (ACC→ACGkbaZtagatose 6‑phosphate aldolase 1, kbaZ subunit
Reads supporting (aligned to +/- strand):  ref base C (43/59);  new base G (16/5);  total (59/64)
Fisher's exact test for biased strand distribution p-value = 7.29e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 128 total aligned reads displayed.

TCCGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTC  >  minE/2104078‑2104217
                                                                     |                                                                      
 ccGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTc                                                                      >  1:22334/1‑71 (MQ=255)
  cGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTc                                                                      >  1:2228587/1‑70 (MQ=255)
       gTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATc                                                                 <  1:4838487/70‑1 (MQ=255)
        tGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAg                                                               <  1:1855657/71‑1 (MQ=255)
          cTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGt                                                             <  1:2121608/71‑1 (MQ=255)
             gAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAg                                                           <  1:4578151/70‑1 (MQ=255)
             gAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGa                                                          >  1:1036656/1‑71 (MQ=255)
              aaGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGAc                                                         <  1:2397512/71‑1 (MQ=255)
                 cTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATc                                                                 <  1:5275941/60‑1 (MQ=255)
                 cTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                      <  1:5440345/71‑1 (MQ=255)
                    cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                   >  1:1552048/1‑71 (MQ=255)
                    cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                   <  1:2613918/71‑1 (MQ=255)
                    cAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                   >  1:4476016/1‑71 (MQ=255)
                       tCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggc                                                <  1:183417/71‑1 (MQ=255)
                       tCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggc                                                >  1:3645520/1‑71 (MQ=255)
                       tCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggc                                                <  1:1224130/71‑1 (MQ=255)
                         cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGAc                                                         >  1:3255240/1‑60 (MQ=255)
                         cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGAc                                                         >  1:1308532/1‑60 (MQ=255)
                         cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGAc                                                         >  1:4230862/1‑60 (MQ=255)
                         cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcgg                                              >  1:3361451/1‑71 (MQ=255)
                         cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcg                                               <  1:2347729/70‑1 (MQ=255)
                         cGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcg                                               <  1:2947930/70‑1 (MQ=255)
                          gCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGt                                             <  1:390094/71‑1 (MQ=255)
                            tACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                           >  1:1221334/1‑71 (MQ=255)
                              cGCCAGTGCAAACCAAACGCCGCTACTGATTGAAGCAACCTCCAATCAGGTAGAc                                                         <  1:2169827/55‑1 (MQ=255)
                                  aGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCgg                                     <  1:4261789/71‑1 (MQ=255)
                                    tGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATc                                                                 <  1:217903/41‑1 (MQ=255)
                                    tGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATc                                                                 <  1:3107623/41‑1 (MQ=255)
                                    tGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAgg                                                              <  1:4713340/44‑1 (MQ=255)
                                    tGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGaa                                   <  1:2162299/71‑1 (MQ=255)
                                      caaaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTa                                                            >  1:3983490/1‑44 (MQ=255)
                                      caaaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTa                                                            >  1:874942/1‑44 (MQ=255)
                                       aaaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGa                                <  1:961907/71‑1 (MQ=255)
                                        aaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAg                                                           <  1:447841/43‑1 (MQ=255)
                                        aaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAg                                                           <  1:833901/43‑1 (MQ=255)
                                        aaccaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                           >  1:1162522/1‑59 (MQ=255)
                                         accaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcg                                                  >  1:268734/1‑51 (MQ=255)
                                         accaaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcg                              >  1:3833679/1‑71 (MQ=255)
                                           caaacGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgcc                            >  1:4833883/1‑71 (MQ=255)
                                              acGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtata                                         >  1:5337673/1‑55 (MQ=255)
                                               cGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgcccgcc                        >  1:2964781/1‑71 (MQ=255)
                                                gCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCccgccg                       <  1:5108881/71‑1 (MQ=255)
                                                  cGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                     >  1:688022/1‑71 (MQ=255)
                                                   gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                    <  1:5237288/71‑1 (MQ=255)
                                                   gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                    <  1:4663429/71‑1 (MQ=255)
                                                   gTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                    <  1:4670184/71‑1 (MQ=255)
                                                    ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                      <  1:819041/36‑1 (MQ=255)
                                                    ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                      <  1:2823204/36‑1 (MQ=255)
                                                    ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAg                                                      <  1:872318/36‑1 (MQ=255)
                                                    ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGt                                             <  1:5040655/45‑1 (MQ=255)
                                                    ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtata                                         <  1:5157905/49‑1 (MQ=255)
                                                    ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAttt                   <  1:4517395/71‑1 (MQ=255)
                                                    ttACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAtt                    <  1:571206/70‑1 (MQ=255)
                                                     tgcTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTta                                           <  1:5381246/44‑1 (MQ=255)
                                                     tgcTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGActtt                  >  1:1452742/3‑71 (MQ=255)
                                                     tACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCg                                      <  1:1515090/51‑1 (MQ=255)
                                                      aCTGATTGAAGCAACCTCCAATCAGGTAGACCAga                                                     >  1:3201554/1‑34 (MQ=255)
                                                      aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTcggcg                                               >  1:3768371/1‑41 (MQ=255)
                                                      aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTc                                                   >  1:4137225/1‑37 (MQ=255)
                                                      aCTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtata                                         >  1:3605181/1‑47 (MQ=255)
                                                       cTGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTtata                                         >  1:4144947/1‑46 (MQ=255)
                                                        tGATTGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGActt                   <  1:4642488/67‑1 (MQ=255)
                                                          aTTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATAcc                                       <  1:4982771/45‑1 (MQ=255)
                                                          aTTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                  <  1:4570683/50‑1 (MQ=255)
                                                          aTTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                  <  1:1889668/50‑1 (MQ=255)
                                                           ttGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaa            >  1:2350933/1‑69 (MQ=255)
                                                           ttGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaa            >  1:2560392/1‑69 (MQ=255)
                                                           ttGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                  <  1:2531918/49‑1 (MQ=255)
                                                           ttGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGc            >  1:5301465/1‑71 (MQ=255)
                                                            tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCgg                                     >  1:2220899/1‑45 (MQ=255)
                                                            tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat           >  1:4235724/1‑68 (MQ=255)
                                                            tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat           >  1:5112299/1‑68 (MQ=255)
                                                            tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat           >  1:2493561/1‑68 (MQ=255)
                                                            tGAAGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTcgcgaat           >  1:3330137/1‑68 (MQ=255)
                                                            tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtat                                          <  1:4668275/40‑1 (MQ=255)
                                                            tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTtat                                          <  1:4066572/40‑1 (MQ=255)
                                                            tGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACTCCCGCCGATTTTCGCgg             >  1:189999/1‑69 (MQ=255)
                                                             gAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                           >  1:2316795/1‑38 (MQ=255)
                                                               aGCAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGTAATGACACCGGCAGACTTTCGCGAAtttg        <  1:1525802/71‑1 (MQ=255)
                                                               aGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                           >  1:1862434/1‑36 (MQ=255)
                                                               aGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                           >  1:1416292/1‑36 (MQ=255)
                                                               aGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTta                                           >  1:1977596/1‑36 (MQ=255)
                                                                 cAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGaa                                   >  1:1464212/1‑42 (MQ=255)
                                                                 cAACGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGaa                                   >  1:2724506/1‑42 (MQ=255)
                                                                  aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGa                                >  1:491660/1‑44 (MQ=255)
                                                                  aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGa                                >  1:3224518/1‑44 (MQ=255)
                                                                  aaCGTCAAACCAGGTCAATCAATTTGGCGGTTATACCGGAATGACACCGGCAGACTTTCGCGAAtttgttt     <  1:5011078/71‑1 (MQ=255)
                                                                  aaCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt     <  1:5089189/71‑1 (MQ=255)
                                                                  aaCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgttt     <  1:3228111/71‑1 (MQ=255)
                                                                   aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGAcgc                             >  1:982795/1‑46 (MQ=255)
                                                                   aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                     <  1:4117836/54‑1 (MQ=255)
                                                                   aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGAt                     <  1:959601/54‑1 (MQ=255)
                                                                   aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgtttg    <  1:1969733/71‑1 (MQ=255)
                                                                   aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgt       >  1:3450725/1‑68 (MQ=255)
                                                                   aCCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttgt       >  1:4174387/1‑68 (MQ=255)
                                                                    ccTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                  <  1:5441795/40‑1 (MQ=255)
                                                                    ccTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAAt                                  <  1:5025624/40‑1 (MQ=255)
                                                                    ccTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTttgtttgt   <  1:4330402/71‑1 (MQ=255)
                                                                     cTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCtttg        <  1:5583957/65‑1 (MQ=255)
                                                                     cTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTc  <  1:3011683/71‑1 (MQ=255)
                                                                     |                                                                      
TCCGCTGGTGCTGGAAGCTGCAATCCGCTACGCCAGTGCAAACCAAACGCCGTTACTGATTGAAGCAACCTCCAATCAGGTAGACCAGTTCGGCGGTTATACCGGAATGACGCCCGCCGATTTTCGCGGCTTTGTTTGTC  >  minE/2104078‑2104217

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: