Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,485,674 A→T 17.5% S192T (TCG→ACG) ‡ yidC ← cytoplasmic insertase into membrane protein, Sec system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,485,6740AT17.5% 225.3 / 33.1 113S192T (TCG→ACG) ‡yidCcytoplasmic insertase into membrane protein, Sec system
Reads supporting (aligned to +/- strand):  ref base A (48/45);  new base T (7/13);  total (55/58)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

GCTGCTTCCGGTATCGAGATGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAGTTGACGTTGACAGCGTA  >  minE/2485615‑2485734
                                                           |                                                            
gctgctTCCGGTATCGAGATGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAAGGTGGAGAtt                                                       >  1:2003316/1‑67 (MQ=255)
gctgctTCCGGTATCGAGATGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAAGGTGGAGAtt                                                       >  1:1776584/1‑67 (MQ=255)
gctgctTCCGGTATCGAGATGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAAGGTGGAGAtt                                                       >  1:1628672/1‑67 (MQ=255)
             tCGAGATGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTtctc                                         >  1:2701595/1‑68 (MQ=255)
             tCGAGATCCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTtctc                                         >  1:180431/1‑68 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGcc                                      >  1:3114082/1‑69 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGcc                                      >  1:1859009/1‑69 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGcc                                      >  1:3353352/1‑69 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGcc                                      >  1:3147639/1‑69 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGcc                                      >  1:732514/1‑69 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGc                                       >  1:288003/1‑68 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGc                                       >  1:1345762/1‑68 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGc                                       >  1:3019530/1‑68 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAg                                    >  1:522373/1‑68 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGc                                   >  1:2559267/1‑69 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGc                                   >  1:3396427/1‑69 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGc                                   >  1:2230534/1‑69 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGc                                   >  1:623282/1‑69 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:565879/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:1846582/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:3012156/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:1695937/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:2165570/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:2069188/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:620832/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:178034/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:2243591/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:2282532/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:1062960/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:362398/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:2316142/69‑1 (MQ=255)
                     cGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGtt                                <  1:2273039/69‑1 (MQ=255)
                     cGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGtt                                <  1:794734/69‑1 (MQ=255)
                            agagTGATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCGGCg                                                >  1:549505/1‑46 (MQ=255)
                            agagTGATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCAc                          >  1:1869480/1‑68 (MQ=255)
                            agagTGATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACg                         >  1:1863269/1‑69 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:19994/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:3208421/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:2941719/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:2938195/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:585082/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:3331625/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:1294191/1‑67 (MQ=255)
                                 gATGGATTGCTTCTACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCTGCGTTCTGCACGTTGTa                    <  1:1627918/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:3133929/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:3104844/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1061567/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2608909/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1345180/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:244539/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1961399/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:513212/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:229701/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2377246/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1461552/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2475479/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:3402843/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1160764/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1106280/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:3264489/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1852559/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:382321/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2966609/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:537395/69‑1 (MQ=255)
                                 gATGGAATGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:661073/69‑1 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGt                     >  1:3056402/1‑67 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:3316739/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:677007/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:2436163/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:3007732/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:275362/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:1037846/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:1339033/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:930968/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:2304671/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:2245465/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:1892561/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:1913605/1‑69 (MQ=255)
                                       ttGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgac              >  1:1035684/1‑69 (MQ=255)
                                       ttGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgac              >  1:1074775/1‑69 (MQ=255)
                                       ttGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgac              >  1:1986660/1‑69 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTtc                                           >  1:2200526/1‑38 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCATGTTGTAgttgagaa            >  1:295003/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:1152541/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:2592541/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:668228/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:736301/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:1002396/1‑66 (MQ=255)
                                           ttCAACTGATCAAACGAGGAGATTTCCAGGGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:924256/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:2887374/46‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:814410/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:963769/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:97425/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:818243/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:487089/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:1030172/68‑1 (MQ=255)
                                            tCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:1886051/67‑1 (MQ=255)
                                                     cAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAGTTGACGTTGACAGCGTa  <  1:2676677/67‑1 (MQ=255)
                                                     cAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAGTTGACGTTGACAGCGTa  <  1:2333703/67‑1 (MQ=255)
                                                          gaggagATTTCCAGCTGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1578834/44‑1 (MQ=255)
                                                           |                                                            
GCTGCTTCCGGTATCGAGATGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAGTTGACGTTGACAGCGTA  >  minE/2485615‑2485734

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: