Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 53,750 C→T 35.9% K346K (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
*minE53,7500CT35.9% 106.4 / 90.1 150K346K (AAG→AAAimpexported protein required for envelope biosynthesis and integrity
Reads supporting (aligned to +/- strand):  ref base C (42/54);  new base T (22/32);  total (64/86)
Fisher's exact test for biased strand distribution p-value = 7.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

Only 100 of 157 total aligned reads displayed.

GAAGTTTTGCACCGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACAC  >  minE/53684‑53818
                                                                  |                                                                    
gAAGTTTTGCACCGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACtt                                                                    >  1:1694744/1‑69 (MQ=255)
gAAGTTTTGCACCGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACtt                                                                    >  1:2381701/1‑69 (MQ=255)
gAAGTTTTGCACCGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACtt                                                                    >  1:2149549/1‑69 (MQ=255)
   gTTTTGCACCGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGt                                                                  <  1:855655/68‑1 (MQ=255)
      ttGCACCGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTAt                                                               >  1:2746527/1‑68 (MQ=255)
           ccGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAAt                                                         <  1:2467716/69‑1 (MQ=255)
           ccGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAAt                                                         <  1:3337994/69‑1 (MQ=255)
           ccGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAAt                                                         <  1:3054471/69‑1 (MQ=255)
           ccGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAAt                                                         <  1:353282/69‑1 (MQ=255)
           ccGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAAt                                                         <  1:2183345/69‑1 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:2567812/1‑68 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:2631282/1‑68 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:986253/68‑1 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:215365/68‑1 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:2912273/68‑1 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:2313863/68‑1 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:2288229/1‑68 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:15113/1‑68 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   >  1:591515/1‑68 (MQ=255)
                  gCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGa                                                   <  1:754405/68‑1 (MQ=255)
                    cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:1312218/1‑69 (MQ=255)
                    cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:2499399/1‑69 (MQ=255)
                    cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:2819090/1‑69 (MQ=255)
                    cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:3125072/1‑69 (MQ=255)
                    cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:1772788/1‑69 (MQ=255)
                    cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:960673/1‑69 (MQ=255)
                    cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:1496126/1‑69 (MQ=255)
                    cAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGt                                                >  1:1672675/1‑69 (MQ=255)
                           gAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAgg                                         <  1:1021380/69‑1 (MQ=255)
                           gAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAgg                                         <  1:343134/69‑1 (MQ=255)
                           gAATTTTTGCGTTGCGTAGCCGGCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAgg                                         <  1:2075375/69‑1 (MQ=255)
                                       tgtgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2647594/66‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:869080/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:697937/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2313475/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:1017530/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:65872/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3380199/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2977069/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:335326/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3144948/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3179502/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3223025/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:3234362/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2270104/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:2166252/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:1125598/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:192655/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:145483/69‑1 (MQ=255)
                                       tgcgtAGCCGTCGCTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:1178663/69‑1 (MQ=255)
                                        gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:550011/68‑1 (MQ=255)
                                        gcgtAGCCGTCGGTACTGGAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACttt                             <  1:631914/68‑1 (MQ=255)
                                        gcgtAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTg                            >  1:3060091/1‑69 (MQ=255)
                                                   aGTACTGGAACCGTACTTGTTATCGAAATCTTTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCgc                  >  1:620848/1‑67 (MQ=255)
                                                   aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:696079/1‑69 (MQ=255)
                                                   aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:2014728/1‑69 (MQ=255)
                                                   aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:1216823/1‑69 (MQ=255)
                                                   aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:1515753/1‑69 (MQ=255)
                                                   aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:319390/1‑69 (MQ=255)
                                                   aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:250882/1‑69 (MQ=255)
                                                   aGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGAc                 >  1:1868940/1‑69 (MQ=255)
                                                     tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:1245649/69‑1 (MQ=255)
                                                     tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:529168/69‑1 (MQ=255)
                                                     tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:1193483/69‑1 (MQ=255)
                                                     tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:1757760/69‑1 (MQ=255)
                                                     tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:3126507/69‑1 (MQ=255)
                                                     tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:1008402/69‑1 (MQ=255)
                                                     tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:293543/69‑1 (MQ=255)
                                                     tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:200293/69‑1 (MQ=255)
                                                     tACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTACTCGACGt               <  1:2571563/69‑1 (MQ=255)
                                                      aCTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               <  1:1387042/68‑1 (MQ=255)
                                                       cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt              >  1:899636/1‑68 (MQ=255)
                                                       cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt              >  1:193140/1‑68 (MQ=255)
                                                       cTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGtt              >  1:1965963/1‑68 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTTGGATCGCTGACTTTg                            >  1:834786/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAgg                                         >  1:662616/1‑39 (MQ=38)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:1919040/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:2055972/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:2022502/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:1337653/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:432295/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:387248/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:1452030/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:2765329/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:1569004/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:3195544/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:2820069/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTGAAGTAGCTAGGATCGCTGACTTTg                            >  1:1770557/1‑52 (MQ=255)
                                                         ggAACCGTATTTGTTATCGAAGTCATTCAAGTAGCTAGGATCGCAGAct                               >  1:2600137/1‑48 (MQ=255)
                                                         ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:2698557/68‑1 (MQ=255)
                                                         ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:132714/68‑1 (MQ=255)
                                                         ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:2525945/68‑1 (MQ=255)
                                                         ggAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGAGCGCTGACCTTGGTGTAGTCGACGTTGa            <  1:3141628/68‑1 (MQ=255)
                                                             ccGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGc       <  1:2248658/69‑1 (MQ=255)
                                                             ccGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGc       <  1:2608109/69‑1 (MQ=255)
                                                                  cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGt               >  1:806029/1‑56 (MQ=255)
                                                                  cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:1606981/1‑69 (MQ=255)
                                                                  cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:1836879/1‑69 (MQ=255)
                                                                  cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:2069906/1‑69 (MQ=255)
                                                                  cTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCacac  >  1:1378874/1‑69 (MQ=255)
                                                                  |                                                                    
GAAGTTTTGCACCGCATAGCCAACGCTGAATTTTTGCGTTGCGTAGCCGTCAGTACTGGAACCGTACTTGTTATCGAAATCATTGAAGTAGCTAGGATCGCTGACCTTGGTGTAGTCGACGTTGAAACGCCACAC  >  minE/53684‑53818

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: