Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 54,086 C→T 30.8% K234K (AAG→AAA imp ← exported protein required for envelope biosynthesis and integrity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE54,0860CT30.8% 78.5 / 38.9 84K234K (AAG→AAAimpexported protein required for envelope biosynthesis and integrity
Reads supporting (aligned to +/- strand):  ref base C (26/32);  new base T (5/21);  total (31/53)
Fisher's exact test for biased strand distribution p-value = 2.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAA  >  minE/54023‑54143
                                                               |                                                         
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:160594/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:313040/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:2805350/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:1111194/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:729339/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:2379285/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:1294935/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:3040957/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:2901126/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:3001086/1‑69 (MQ=255)
cATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGc                                                      >  1:1448008/1‑69 (MQ=255)
        gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCggg                                               <  1:1337146/68‑1 (MQ=255)
        gcgcGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCggg                                               <  1:2601500/68‑1 (MQ=255)
              tGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTTGTGTACTTGGCGTTCGGGATCAAGa                                        <  1:2814336/69‑1 (MQ=255)
               gTTCCAGTAATATGGCAGGTAGAACTCATAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa                                       <  1:2460011/69‑1 (MQ=255)
               gTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa                                       <  1:992725/69‑1 (MQ=255)
               gTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa                                       <  1:2760891/69‑1 (MQ=255)
               gTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa                                       <  1:812859/69‑1 (MQ=255)
               gTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa                                       <  1:731245/69‑1 (MQ=255)
               gTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGaa                                       <  1:1531723/69‑1 (MQ=255)
                     gTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCaa                                          >  1:1950849/1‑60 (MQ=255)
                     gTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCaga                                 >  1:1695386/1‑69 (MQ=255)
                     gTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCaga                                 >  1:2575683/1‑69 (MQ=255)
                     gTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCag                                  >  1:818643/1‑68 (MQ=255)
                     gTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCa                                   >  1:2981791/1‑67 (MQ=255)
                        atatGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGc                               >  1:2205206/1‑68 (MQ=255)
                        atatGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCg                              >  1:3211814/1‑69 (MQ=255)
                        atatGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCg                              >  1:2598022/1‑69 (MQ=255)
                        atatGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCg                              >  1:1508778/1‑69 (MQ=255)
                        atatGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCg                              >  1:1723776/1‑69 (MQ=255)
                           tGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACg                           <  1:1086374/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTTGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:2232267/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGTTCAAGAAACCAGAGCGACGtt                         <  1:1422658/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:104418/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:389466/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:634834/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:804507/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:78268/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:602677/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:536535/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:3242406/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:3279948/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:1031323/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:2303440/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:1216517/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:2397895/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:2420747/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:2510832/69‑1 (MQ=255)
                             gCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:2526719/69‑1 (MQ=255)
                              cAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTATTTGGCGTTCGGGTTCAAGAAACCAGAGCGACGtt                         <  1:1767386/68‑1 (MQ=255)
                                  tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTc                                                  >  1:1108060/1‑39 (MQ=255)
                                  tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTc                                                  >  1:3320548/1‑39 (MQ=255)
                                  tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGAt                                             >  1:595702/1‑44 (MQ=255)
                                  tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCa                                           >  1:3018547/1‑46 (MQ=255)
                                  tAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTca                    >  1:989566/1‑69 (MQ=255)
                                     aaCTCAAAGTAGTTGGTTGTGGTGTATTTGGCGTTCGGTATCAAGAAACCAGAGCGACttt                         <  1:1115190/61‑1 (MQ=255)
                                       cTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTcacccac               >  1:29638/1‑69 (MQ=255)
                                       cTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTcacccac               >  1:532765/1‑69 (MQ=255)
                                       cTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTcacccac               >  1:1528282/1‑69 (MQ=255)
                                       cTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTcacccac               >  1:986875/1‑69 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:1604022/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:649578/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:2626637/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:1967200/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:2673706/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:2439173/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:2783271/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:3093598/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:2283941/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:2252995/68‑1 (MQ=255)
                                            agtagtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGc           <  1:2127890/68‑1 (MQ=255)
                                               agtTGGTGGTGGTGTATTTGGCGTTCGGGATCAAGAAACCAGAGCGACGtt                         <  1:854054/51‑1 (MQ=255)
                                               agtTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAAc        >  1:1850011/1‑68 (MQ=255)
                                                  tggtggtggTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCa    <  1:3262893/69‑1 (MQ=255)
                                                  tggtggtggTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCa    <  1:334180/69‑1 (MQ=255)
                                                  tggtggtggTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCa    <  1:3368284/69‑1 (MQ=255)
                                                  tggtggtggTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCa    <  1:531569/69‑1 (MQ=255)
                                                  tggtggtggTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCa    <  1:792451/69‑1 (MQ=255)
                                                  tggtggtggTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCa    <  1:3092808/69‑1 (MQ=255)
                                                  tggtggtggTGTACTTGGCGTTCGGGAGCAAGAACCCAGAGCGACGTTTGTCACCCACCGGCAACTGCa    <  1:3007227/69‑1 (MQ=255)
                                                     tggtggtgTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCaaa  <  1:1523011/68‑1 (MQ=255)
                                                     tggtggtgTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCaaa  <  1:2427734/68‑1 (MQ=255)
                                                     tggtggtgTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCaaa  <  1:1346620/68‑1 (MQ=255)
                                                     tggtggtgTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCaaa  <  1:44233/68‑1 (MQ=255)
                                                               |                                                         
CATATTTGGCGCGATGTTCCAGTAATATGGCAGGTAGAACTCAAAGTAGTTGGTGGTGGTGTACTTGGCGTTCGGGATCAAGAAACCAGAGCGACGTTTGTCACCCACCGGCAACTGCAAA  >  minE/54023‑54143

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: