Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,314,255 T→C 18.4% H170H (CAT→CAC hisB → fused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,314,2550TC18.4% 282.3 / 57.7 125H170H (CAT→CAChisBfused histidinol‑phosphatase and imidazoleglycerol‑phosphate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (47/55);  new base C (14/9);  total (61/64)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 126 total aligned reads displayed.

TACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAA  >  minE/1314187‑1314323
                                                                    |                                                                    
tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCAt                                                                      >  1:2750104/1‑69 (MQ=255)
tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATg                                                                     >  1:781509/1‑70 (MQ=255)
tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATg                                                                     >  1:136050/1‑70 (MQ=255)
tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATg                                                                     >  1:2781301/1‑70 (MQ=255)
tACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATg                                                                     >  1:1745829/1‑70 (MQ=255)
  cGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTa                                                                   <  1:3034328/70‑1 (MQ=255)
  cGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGt                                                                    <  1:629060/69‑1 (MQ=255)
     ccGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGc                                                               <  1:1628333/71‑1 (MQ=255)
       gcgAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGt                                                             <  1:1917839/71‑1 (MQ=255)
          aaaCCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTaa                                                           <  1:1936855/70‑1 (MQ=255)
          aaaCCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTaa                                                           <  1:2910389/70‑1 (MQ=255)
            aCCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAAt                                                          >  1:3772008/1‑69 (MQ=255)
             cccTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTaa                                                           >  1:4130530/1‑67 (MQ=255)
              ccTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCa                                                      <  1:825592/71‑1 (MQ=255)
               cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                                                    <  1:2877820/56‑1 (MQ=255)
               cTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGt                                                                    <  1:2304491/56‑1 (MQ=255)
                 gAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                                                           >  1:2338797/1‑63 (MQ=255)
                   aCTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAag                                                   <  1:757001/69‑1 (MQ=255)
                    cTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAgaga                                                >  1:3933083/1‑71 (MQ=255)
                    cTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGAGTAATACCAAAgaga                                                >  1:2493918/1‑71 (MQ=255)
                      ggCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACg                                              >  1:3823304/1‑71 (MQ=255)
                       gCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGCCg                                              >  1:3711809/1‑70 (MQ=255)
                        ccAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCa                                            <  1:2115196/71‑1 (MQ=255)
                         cAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACg                                              >  1:2823416/1‑68 (MQ=255)
                          aaTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTACTACCAAAGAGACg                                              >  1:2598809/1‑67 (MQ=255)
                           aTGATTGGCGAGCAACTCACCAGACTTGACCGTTACGCTCATGTAGTGCGTaa                                                           >  1:92977/1‑53 (MQ=255)
                           aTGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGAt                                         <  1:2897383/71‑1 (MQ=255)
                            tGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATg                                                                     >  1:4137841/1‑42 (MQ=255)
                            tGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTa                                                                   >  1:2069914/1‑44 (MQ=255)
                            tGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAg                                                                  >  1:3985364/1‑45 (MQ=255)
                            tGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAg                                                                  >  1:2833321/1‑45 (MQ=255)
                              aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTaa                                                           >  1:1024984/1‑50 (MQ=255)
                              aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTg                                       <  1:2355388/70‑1 (MQ=255)
                              aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTg                                       <  1:3927231/70‑1 (MQ=255)
                              aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTg                                       <  1:3857566/70‑1 (MQ=255)
                              aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                                      >  1:1687928/1‑71 (MQ=255)
                              aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGa                                      >  1:1275704/1‑71 (MQ=255)
                              aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGt                                                             <  1:1695105/48‑1 (MQ=255)
                              aTTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTaa                                                           >  1:785557/1‑50 (MQ=255)
                               ttGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTg                                       <  1:1182371/69‑1 (MQ=255)
                               ttGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                     >  1:3451305/1‑71 (MQ=255)
                               ttGGCGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                     >  1:897130/1‑71 (MQ=255)
                                tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                     >  1:3183148/1‑70 (MQ=255)
                                tGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                     >  1:2162095/1‑70 (MQ=255)
                                 ggCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTg                                       <  1:487166/67‑1 (MQ=255)
                                 ggCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTg                                       <  1:113065/67‑1 (MQ=255)
                                 ggCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTg                                       <  1:3658369/67‑1 (MQ=255)
                                 ggCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTg                                       <  1:355644/67‑1 (MQ=255)
                                 ggCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGt                                   >  1:145848/1‑71 (MQ=255)
                                  gCGAGCAACTCACCAGACGTGACCGTTACGCTCATg                                                                     <  1:1426739/36‑1 (MQ=255)
                                   cGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATa                                                         >  1:2018943/1‑47 (MQ=255)
                                   cGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATa                                                         >  1:2611087/1‑47 (MQ=255)
                                   cGAGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTc                                 <  1:2224189/71‑1 (MQ=255)
                                     aGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAg                               <  1:4066028/71‑1 (MQ=255)
                                     aGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAg                               <  1:3776668/71‑1 (MQ=255)
                                     aGCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAg                               <  1:1951248/71‑1 (MQ=255)
                                      gCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCaaa                                                    >  1:2827190/1‑49 (MQ=255)
                                      gCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                     >  1:452290/1‑64 (MQ=255)
                                      gCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAgg                              <  1:2772566/71‑1 (MQ=255)
                                      gCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAgg                              <  1:3250726/71‑1 (MQ=255)
                                      gCAACTCACCAGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAgg                              >  1:2699141/1‑71 (MQ=255)
                                        aaCTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgtg                            <  1:2210775/71‑1 (MQ=255)
                                         aCTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgtgt                           <  1:598174/71‑1 (MQ=255)
                                         aCTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGgtgt                           <  1:2503197/71‑1 (MQ=255)
                                          cacACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                          >  1:1998222/3‑71 (MQ=255)
                                          cTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAgaga                                                >  1:2461868/1‑49 (MQ=255)
                                          cTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAgaga                                                <  1:4007307/49‑1 (MQ=255)
                                          cTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCa                                            <  1:2480623/53‑1 (MQ=255)
                                          cTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                          >  1:2566771/1‑71 (MQ=255)
                                          cTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                          <  1:1496360/71‑1 (MQ=255)
                                          cTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGtgtg                          >  1:3923354/1‑71 (MQ=255)
                                            cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACg                                              >  1:1584824/1‑49 (MQ=255)
                                            cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGc                        >  1:325090/1‑71 (MQ=255)
                                            cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGc                        >  1:1292661/1‑71 (MQ=255)
                                            cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGc                        >  1:3216261/1‑71 (MQ=255)
                                            cACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGc                        >  1:320075/1‑71 (MQ=255)
                                             aCCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTgg                         <  1:3153055/69‑1 (MQ=255)
                                             aCCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTgg                         <  1:4143712/69‑1 (MQ=255)
                                             aCCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCt                       >  1:3767103/1‑71 (MQ=255)
                                                aGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGAc                                     >  1:4008615/1‑54 (MQ=255)
                                                aGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGa                    <  1:3900635/71‑1 (MQ=255)
                                                aGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGa                    <  1:1834016/71‑1 (MQ=255)
                                                aGACGTGACCGTTACGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGa                    >  1:3643766/1‑71 (MQ=255)
                                                  aCGTGACCGTTACGCTCATGTAGTGCGTAATACCaaa                                                    <  1:2122438/37‑1 (MQ=255)
                                                  aCGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATc                  <  1:2712540/71‑1 (MQ=255)
                                                   cGTGACCGTTACGCTCATGTAGTGCGTAATACCAAag                                                   <  1:878925/37‑1 (MQ=255)
                                                   cGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATc                  >  1:399600/1‑70 (MQ=255)
                                                   cGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATc                  <  1:3697173/70‑1 (MQ=255)
                                                       aCCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGTCTGGATCGTGaa             <  1:1125249/71‑1 (MQ=255)
                                                          gTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACg                                    <  1:3055327/45‑1 (MQ=255)
                                                          gTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAgg                              <  1:306051/51‑1 (MQ=255)
                                                          gTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAgg                              <  1:3175574/51‑1 (MQ=255)
                                                             aCGCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTgg        >  1:2095505/1‑70 (MQ=255)
                                                               gCTCACGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgcag     >  1:188870/1‑71 (MQ=255)
                                                                cTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgcagc    <  1:1802161/71‑1 (MQ=255)
                                                                cTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgcagc    <  1:1439660/71‑1 (MQ=255)
                                                                cTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGgcagc    <  1:3316506/71‑1 (MQ=255)
                                                                 tCATGTAGTTCGTAATACCAAAGAGACGCAGATTGACGtt                                  <  1:1267507/40‑1 (MQ=255)
                                                                   aTGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCaa  <  1:2987186/70‑1 (MQ=255)
                                                                   aTGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCaa  <  1:3378783/70‑1 (MQ=255)
                                                                    |                                                                    
TACGACCGCGAAACCCTGAACTGGCCAATGATTGGCGAGCAACTCACCAGACGTGACCGTTACGCTCATGTAGTGCGTAATACCAAAGAGACGCAGATTGACGTTCAGGTGTGGCTGGATCGTGAAGGTGGCAGCAA  >  minE/1314187‑1314323

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: