Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 193,541 A→G 28.0% V142V (GTA→GTG yaeQ → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE193,5410AG28.0% 138.6 / 58.2 133V142V (GTA→GTGyaeQconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (44/51);  new base G (25/12);  total (69/64)
Fisher's exact test for biased strand distribution p-value = 3.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 134 total aligned reads displayed.

TGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTG  >  minE/193486‑193600
                                                       |                                                           
tgCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCg                                                <  1:2567305/69‑1 (MQ=255)
tgCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCg                                                <  1:1305051/69‑1 (MQ=255)
tgCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCg                                                <  1:2772439/69‑1 (MQ=255)
  cAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGAt                                              <  1:2921213/69‑1 (MQ=255)
  cAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGAt                                              <  1:24839/69‑1 (MQ=255)
  cAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGAt                                              <  1:151048/69‑1 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTa                                          >  1:176244/1‑68 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:2752304/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:3198483/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:3030052/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:2897254/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:2158222/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:23479/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:391002/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:472404/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:60292/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:11639/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:2702906/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTAc                                         >  1:1073935/1‑69 (MQ=255)
       tGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCATTTGCCGATCGTAc                                         >  1:685934/1‑69 (MQ=255)
            aTCTTTCCGTCTGGTATCTGGACGACGAGCAGCTGGCGAAAGTGAGCGGCTTTGCCGATCGt                                           <  1:826673/62‑1 (MQ=255)
            aTCTTTCCGTCTGGTATCTGGACGACGAGCAGCTGGCGAAAGTGAGCGGCTTTGCCGATCGt                                           <  1:3189650/62‑1 (MQ=255)
            aTCTTTCCGTCTGGTATCTGGACGACGAGCAGCTGGCGAAAGTGAGCGGCTTTGCCGATCGt                                           <  1:2643505/62‑1 (MQ=255)
            aTCTTTCCGTCTGGTATCTGGACGACGAGCAGCTGGCGAAAGTGAGCGGCTTTGCCGATCGt                                           <  1:1882635/62‑1 (MQ=255)
            aTCTTTCCGTCTGGTATCTGGACGACGAGCAGCTGGCGAAAGTGAGCGGCTTTGCCGATCGt                                           <  1:342669/62‑1 (MQ=255)
                         gTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAAc                       >  1:3140582/1‑69 (MQ=255)
                         gTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAAc                       >  1:2206127/1‑69 (MQ=255)
                           aTCTGTACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:1575139/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:1829103/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:670843/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:1247239/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:2733746/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:1397120/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:2854944/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:2412294/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:2235593/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:3067780/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:2041051/1‑69 (MQ=255)
                           aTCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:3248398/1‑69 (MQ=255)
                           aTCTGGACGATCAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGa                     >  1:1739696/1‑69 (MQ=255)
                              tGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCAt                                      >  1:3370014/1‑49 (MQ=255)
                              tGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCAt                                      >  1:442698/1‑49 (MQ=255)
                              tGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCAt                                      >  1:3366171/1‑49 (MQ=255)
                                  cGATGAACAACTGGCGATAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:761528/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTTAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:1058117/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:174677/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:155679/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:3402571/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:3018684/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:1944169/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:2131150/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:1673246/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:1112715/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:101675/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:2412463/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:2422401/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:925189/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:3314727/67‑1 (MQ=255)
                                  cGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACAATGACGCTGCAGGCAACGATTCAg                <  1:422183/67‑1 (MQ=255)
                                       aacaacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCAt                                      <  1:1902372/40‑1 (MQ=255)
                                       aacaacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCAt                                      <  1:943832/40‑1 (MQ=255)
                                       aacaacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCAt                                      <  1:3127798/40‑1 (MQ=255)
                                       aacaacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCAt                                      <  1:2070891/40‑1 (MQ=255)
                                       aacaacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCAt                                      <  1:942562/40‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:1382722/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:536248/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:580346/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:642431/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:1151321/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:726942/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:2703302/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:3337668/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:3105751/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:311003/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:2003810/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:2872693/69‑1 (MQ=255)
                                         caacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:2158057/69‑1 (MQ=255)
                                          aacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:2939773/68‑1 (MQ=255)
                                          aacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:577579/68‑1 (MQ=255)
                                          aacTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTg       <  1:513243/68‑1 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCa                 >  1:1553783/1‑56 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTAt     >  1:1795310/1‑68 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:2603917/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:2865937/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:907737/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:1202135/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:2479772/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:2479810/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:1373443/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:496116/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:3140624/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:1671320/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:3158303/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:3363667/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGTGTTATc    >  1:3356817/1‑69 (MQ=255)
                                            cTGGCGAAAGTGAGCGGCTTTGCCGATCGTACAATGACGCTGCAGGCAACGATTCAGGAt             >  1:2984291/1‑60 (MQ=255)
                                              ggCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTg  >  1:347132/1‑69 (MQ=255)
                                              ggCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTg  >  1:3318168/1‑69 (MQ=255)
                                                 gAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                <  1:2114640/52‑1 (MQ=255)
                                                       aaGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAg                >  1:40859/1‑46 (MQ=255)
                                                       |                                                           
TGCAGTTTGCCAATCTTTCCGTCTGGTATCTGGACGATGAACAACTGGCGAAAGTAAGCGCCTTTGCCGATCGTACCATGACGCTGCAGGCAACGATTCAGGATGGCGTGATCTG  >  minE/193486‑193600

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: