Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 581,798 T→C 44.4% T612T (ACT→ACC clpA → ATPase and specificity subunit of ClpA‑ClpP ATP‑dependent serine protease, chaperone activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE581,7980TC44.4% 24.1 / 103.5 124T612T (ACT→ACCclpAATPase and specificity subunit of ClpA‑ClpP ATP‑dependent serine protease, chaperone activity
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base C (30/25);  total (63/61)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 124 total aligned reads displayed.

AGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCA  >  minE/581741‑581859
                                                         |                                                             
aGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATc                                                    <  1:1395385/69‑1 (MQ=255)
           gTAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCtt                                          <  1:1035535/68‑1 (MQ=255)
           gTAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCtt                                          <  1:2229355/68‑1 (MQ=255)
            tAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTa                                         <  1:918588/68‑1 (MQ=255)
            tAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTa                                         <  1:2753636/68‑1 (MQ=255)
            tAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTa                                         <  1:1794739/68‑1 (MQ=255)
            aaaCGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCtt                                          <  1:3298817/66‑1 (MQ=255)
                   gtgCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAg                                 >  1:3294242/1‑69 (MQ=255)
                   gtgCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAg                                 >  1:3119136/1‑69 (MQ=255)
                   gtgCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAg                                 >  1:2895657/1‑69 (MQ=255)
                   gtgCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAg                                 >  1:28031/1‑69 (MQ=255)
                   gtgCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAg                                 >  1:2469342/1‑69 (MQ=255)
                   gtgCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAg                                 >  1:1369212/1‑69 (MQ=255)
                      cTGGTGATGCCCACCAACGCCGGGGTACGGGAAACTg                                                              >  1:2155386/1‑37 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTg                                                              >  1:1071749/1‑37 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTg                                                              >  1:2011223/1‑37 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTg                                                              >  1:1053796/1‑37 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTg                                                              >  1:2236669/1‑37 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTg                                                              >  1:2485825/1‑37 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTg                                                              >  1:420230/1‑37 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGAt                              <  1:674416/69‑1 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGAt                              <  1:659372/69‑1 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGAt                              <  1:2796444/69‑1 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGAt                              <  1:2965256/69‑1 (MQ=255)
                      cTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGAt                              <  1:1906688/69‑1 (MQ=255)
                        ggTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGAt                              <  1:202022/67‑1 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCtt                                          <  1:1225352/48‑1 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCa                       >  1:3036559/1‑67 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:3381981/69‑1 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     >  1:2398979/1‑69 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     >  1:2219953/1‑69 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     >  1:1990209/1‑69 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     >  1:1883744/1‑69 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     >  1:1313720/1‑69 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:1326975/69‑1 (MQ=255)
                               accaccAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     >  1:827506/1‑69 (MQ=255)
                                ccactaACGCGGGGTTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:3295978/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     >  1:2552793/1‑68 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:1947632/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:330719/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:1138260/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:508367/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:2957405/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:2899323/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:1352358/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:2760456/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:1469850/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:1799122/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:193612/1‑69 (MQ=255)
                                ccactaACGCGGGGGTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACcg                    >  1:780280/1‑69 (MQ=255)
                                ccaccaACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:1074978/68‑1 (MQ=255)
                                   ccaACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTa                                         <  1:2044708/45‑1 (MQ=255)
                                   ccaACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:853260/65‑1 (MQ=255)
                                           ggggTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGaa         <  1:3252351/69‑1 (MQ=255)
                                           ggggTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGaa         <  1:1152389/69‑1 (MQ=255)
                                            gggTACGTGAAACCTAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:2626416/1‑66 (MQ=255)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATc                                       >  1:2076625/1‑38 (MQ=38)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:1841462/1‑66 (MQ=255)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:2084157/1‑66 (MQ=255)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:1183935/1‑66 (MQ=255)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:1652518/1‑66 (MQ=255)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:1030763/1‑66 (MQ=255)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:3412027/1‑66 (MQ=255)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:3370239/1‑66 (MQ=255)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:884507/1‑66 (MQ=255)
                                            gggTACGTGAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATgg           >  1:772875/1‑66 (MQ=255)
                                               tACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:562924/53‑1 (MQ=255)
                                               tACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:2435395/53‑1 (MQ=255)
                                               tACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:1655730/53‑1 (MQ=255)
                                               tACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:137757/53‑1 (MQ=255)
                                               tACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:2763173/53‑1 (MQ=255)
                                               tACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:277199/53‑1 (MQ=255)
                                               tACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:118867/53‑1 (MQ=255)
                                               tACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:524094/53‑1 (MQ=255)
                                               tACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCAcc                     <  1:635004/53‑1 (MQ=255)
                                                 cGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAgaga     >  1:1045743/1‑67 (MQ=255)
                                                   tgAAACCGATCTCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:1799439/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGTAAGAGATCa  <  1:682438/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:723792/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:690442/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:559402/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:2667248/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:31846/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:1116254/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:151428/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:1062810/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:1988932/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:2280815/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:2671235/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:2741981/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:2795861/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:3094084/67‑1 (MQ=255)
                                                   tgAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:331988/67‑1 (MQ=255)
                                                   agAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:1414360/67‑1 (MQ=255)
                                                    gAAACCGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  <  1:2881001/67‑1 (MQ=255)
                                                     aaaCTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  >  1:253716/1‑66 (MQ=255)
                                                     aaaCTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  >  1:2903927/1‑66 (MQ=255)
                                                     aaaCTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  >  1:2986293/1‑66 (MQ=255)
                                                     aaaCTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  >  1:1736398/1‑66 (MQ=255)
                                                     aaaCTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCa  >  1:1326647/1‑66 (MQ=255)
                                                         |                                                             
AGCAGACTTCCGTAACGTGGTGCTGGTGATGACCACCAACGCCGGGGTACGGGAAACTGAGCGCAAATCCATTGGTCTTATCCACCAGGATAACAGCACCGATGCGATGGAAGAGATCA  >  minE/581741‑581859

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: