Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,048,438 C→T 45.6% V639V (GTG→GTA glnE ← fused deadenylyltransferase and adenylyltransferase for glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,048,4380CT45.6% 15.7 / 76.7 92V639V (GTG→GTAglnEfused deadenylyltransferase and adenylyltransferase for glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base C (30/20);  new base T (20/22);  total (50/42)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGCCTGCTGAACGACGGCATCTATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTGTTTGAACTG  >  minE/2048370‑2048500
                                                                    |                                                              
cACGCCTGCTGAACGACGGCATCTATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTtcatc                                                                <  1:1211350/69‑1 (MQ=255)
cACGCCTGCTGAACGACGGCATCTATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTtcatc                                                                <  1:2677948/69‑1 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCCGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:1035152/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGt                                                          >  1:557717/1‑68 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGt                                                          >  1:2763407/1‑68 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:3362151/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:161020/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:1795185/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:1972423/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:2142928/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:2667972/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:3051575/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:3053368/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:1017763/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:375993/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:54384/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:707662/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:774058/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:79342/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:1486834/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:1145093/1‑69 (MQ=255)
       gctgAACGACGGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTa                                                         >  1:1106486/1‑69 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCTGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:2046597/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTATGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:2186120/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGctt                                                     <  1:828415/67‑3 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGctt                                                     <  1:824536/67‑3 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:2177871/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:76363/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:736548/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:1096839/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:2564342/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:1408550/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:2856938/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:3312304/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:1274165/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:3115017/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:1492154/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:1674021/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:1999568/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:880724/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:1003756/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:1963672/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:185212/67‑1 (MQ=255)
             cgacgGCATCAATCATGGCTTCTGCCAGCCAGGTTAAGTGATCGCTCACTTTCATTACCGGTAGCGt                                                     <  1:859703/67‑1 (MQ=255)
                   cATCTATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGa                                             <  1:778655/69‑1 (MQ=255)
                   cATCTATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGa                                             <  1:3063976/69‑1 (MQ=255)
                   cATCTATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGa                                             <  1:2194756/69‑1 (MQ=255)
                   cATCTATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGa                                             <  1:15399/69‑1 (MQ=255)
                   cATCTATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGa                                             <  1:2834993/69‑1 (MQ=255)
                       tATCATGGCTTCCGCCAGCCATGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGatat                                          >  1:1645787/1‑68 (MQ=255)
                       tATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATc                                         >  1:1365893/1‑69 (MQ=255)
                       tATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATc                                         >  1:2826521/1‑69 (MQ=255)
                       tATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATc                                         >  1:2830049/1‑69 (MQ=255)
                       tATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATc                                         >  1:2744846/1‑69 (MQ=255)
                       tATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATc                                         >  1:1994869/1‑69 (MQ=255)
                       tATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATc                                         >  1:2885224/1‑69 (MQ=255)
                       tATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATc                                         >  1:1729657/1‑69 (MQ=255)
                       tATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATc                                         >  1:1049681/1‑69 (MQ=255)
                       tATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATc                                         >  1:2034765/1‑69 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATGGCTCACTATCATCACCGGTAGCGTACCGGCGATATCCGCt                                     >  1:1252922/1‑69 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCgg                                                           >  1:434717/1‑47 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCgg                                                           >  1:384548/1‑47 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCgg                                                           >  1:354815/1‑47 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCgg                                                           >  1:2866634/1‑47 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCgg                                                           >  1:3064092/1‑47 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCgg                                                           >  1:2643899/1‑47 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCt                                     >  1:2176015/1‑69 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCt                                     >  1:2257296/1‑69 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCt                                     >  1:1947898/1‑69 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCt                                     >  1:678449/1‑69 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCt                                     >  1:2662691/1‑69 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCt                                     >  1:869957/1‑69 (MQ=255)
                           aTGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCt                                     >  1:2974410/1‑69 (MQ=255)
                                ttCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTgcggc                                <  1:2757371/69‑1 (MQ=255)
                                        gccagGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTa                        <  1:580190/69‑1 (MQ=255)
                                        gccagGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTa                        <  1:1219251/69‑1 (MQ=255)
                                        gccagGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTa                        <  1:1875750/69‑1 (MQ=255)
                                               tAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTg                 >  1:491751/1‑69 (MQ=255)
                                               tAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTg                 >  1:173804/1‑69 (MQ=255)
                                               tAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTg                 >  1:318627/1‑69 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:1832603/69‑1 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:82963/69‑1 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:978826/69‑1 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:14951/69‑1 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:993837/69‑1 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:2418706/69‑1 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:2443956/69‑1 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:523916/69‑1 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:3007611/69‑1 (MQ=255)
                                                     aTCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTg           <  1:3154785/69‑1 (MQ=255)
                                                              tttCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTGTTTGAAc    >  1:1600080/1‑67 (MQ=255)
                                                              tttCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTGTTTGAACt   <  1:1157333/68‑1 (MQ=255)
                                                              tttCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTGTTTGAACTg  >  1:2186882/1‑69 (MQ=255)
                                                              tttCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTGTTTGAACTg  >  1:400660/1‑69 (MQ=255)
                                                                    |                                                              
CACGCCTGCTGAACGACGGCATCTATCATGGCTTCCGCCAGCCAGGTTAAGTGATCGCTCACTTTCATCACCGGTAGCGTACCGGCGATATCCGCTGCGGCGATGCGTAACAGCTGCGCCTGTTTGAACTG  >  minE/2048370‑2048500

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: