Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,485,672 C→G 17.7% S192S (TCG→TCC) ‡ 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,6720CG17.7% 233.4 / 37.1 113S192S (TCG→TCC) ‡yidCcytoplasmic insertase into membrane protein, Sec system
Reads supporting (aligned to +/- strand):  ref base C (48/45);  new base G (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 = 9.96e-01

Only 100 of 113 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)
gctgctTCCGGTATCGAGATGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAAGGTGGAGAtt                                                       >  1:329561/1‑67 (MQ=255)
             tCGAGATGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTtctc                                         >  1:2166800/1‑68 (MQ=255)
             tCGAGATCCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTtctc                                         >  1:180431/1‑68 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGcc                                      >  1:732514/1‑69 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGcc                                      >  1:3114082/1‑69 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGcc                                      >  1:3147639/1‑69 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGcc                                      >  1:3353352/1‑69 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGc                                       >  1:1345762/1‑68 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGc                                       >  1:3019530/1‑68 (MQ=255)
               gagaTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGc                                       >  1:288003/1‑68 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAg                                    >  1:522373/1‑68 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGc                                   >  1:3396427/1‑69 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGc                                   >  1:2230534/1‑69 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGc                                   >  1:623282/1‑69 (MQ=255)
                  aTGCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGc                                   >  1:2559267/1‑69 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:565879/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:1695937/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:2165570/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:593782/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:1846582/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:2069188/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:620832/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTcc                                                    <  1:3012156/50‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:2316142/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:2243591/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:2282532/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:229767/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:1062960/69‑1 (MQ=255)
                    gCGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGt                                 <  1:178034/69‑1 (MQ=255)
                     cGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGtt                                <  1:794734/69‑1 (MQ=255)
                     cGGTGGCAGAGTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGtt                                <  1:2273039/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:2938195/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:19994/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:2941719/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:3331625/1‑67 (MQ=255)
                               gTGATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGt                        >  1:1294191/1‑67 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:3104844/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2377246/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:3133929/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2608909/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:513212/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2475479/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:244539/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1061567/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:229701/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1345180/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1461552/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1961399/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2507461/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1106280/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:537395/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1160764/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:382321/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:3402843/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:3264489/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:1852559/69‑1 (MQ=255)
                                 gATGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2966609/69‑1 (MQ=255)
                                 gATGGAATGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:661073/69‑1 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAAGGTGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:2916499/68‑1 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGt                     >  1:3056402/1‑67 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTa                    <  1:421244/68‑1 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:1913605/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:2304671/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:1037846/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:275362/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:1339033/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:3316739/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:930968/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:1892561/1‑69 (MQ=255)
                                  aTGGATTGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAg                   >  1:677007/1‑69 (MQ=255)
                                       ttGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgac              >  1:1035684/1‑69 (MQ=255)
                                       ttGCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgac              >  1:1074775/1‑69 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTtc                                           >  1:2200526/1‑38 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCATGTTGTAgttgagaa            >  1:295003/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:668228/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:1002396/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:2592541/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:1152541/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:736301/1‑66 (MQ=255)
                                         gCTTCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgagaa            >  1:3045603/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:97425/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:963769/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:818243/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:487089/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:1900610/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:1831995/68‑1 (MQ=255)
                                           ttCAACTGACCAAACGAGGAGATTTCCAGCGGTTTCTCGCCAGCGTTCTGCACGTTGTAgttgacgtt           <  1:1495922/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)
                                                         |                                                              
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: