Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,330,617 G→A 21.2% A323A (GCC→GCT cpsG ← phosphomannomutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,330,6170GA21.2% 211.9 / 41.5 145A323A (GCC→GCTcpsGphosphomannomutase
Reads supporting (aligned to +/- strand):  ref base G (50/64);  new base A (18/13);  total (68/77)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 145 total aligned reads displayed.

GTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCC  >  minE/1330566‑1330682
                                                   |                                                                 
gTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:3113859/1‑69 (MQ=255)
gTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:412660/1‑69 (MQ=255)
gTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:429457/1‑69 (MQ=255)
gTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:2616485/1‑69 (MQ=255)
gTAAGCGAAATCACGGAAGTAATGGTGGCCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCAt                                                  >  1:1390572/1‑69 (MQ=255)
                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           <  1:492084/56‑1 (MQ=255)
                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           <  1:2443155/56‑1 (MQ=255)
                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           >  1:1468083/1‑56 (MQ=255)
                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           <  1:3189745/56‑1 (MQ=255)
                    aaTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           >  1:3005646/1‑56 (MQ=255)
                     aTGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCt                                           <  1:645253/55‑1 (MQ=255)
                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:2832139/1‑69 (MQ=255)
                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:27379/1‑69 (MQ=255)
                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:2498442/1‑69 (MQ=255)
                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:2222746/1‑69 (MQ=255)
                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:1473748/1‑69 (MQ=255)
                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:2981893/1‑69 (MQ=255)
                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTcc                          >  1:3115308/1‑69 (MQ=255)
                        gtggGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTc                           >  1:1764713/1‑68 (MQ=255)
                                tCATTTCGCCACCAT‑‑ATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTTCGGTTTTCGAc                <  1:244677/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCg                  >  1:929843/1‑68 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2365334/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:324172/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 >  1:2293168/1‑69 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:3246428/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:3267664/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2458292/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2632478/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2639810/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1035737/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:3207122/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:294780/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2920587/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:2795751/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1798146/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1185038/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:600266/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:823559/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:90702/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 >  1:165521/1‑69 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1745072/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:114423/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1819249/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1103832/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1895164/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:1923026/69‑1 (MQ=255)
                                 cATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:3271967/69‑1 (MQ=255)
                                  aTTTCGCCACCTTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTGTCGa                 <  1:3311737/68‑1 (MQ=255)
                                  aTTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 <  1:197060/68‑1 (MQ=255)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCttt                                         >  1:595270/1‑42 (MQ=38)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:164412/1‑45 (MQ=255)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:3063721/1‑45 (MQ=255)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:1176114/1‑45 (MQ=255)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:1012224/1‑45 (MQ=255)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:507133/1‑45 (MQ=255)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:2701262/1‑45 (MQ=255)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:2145298/1‑45 (MQ=255)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:2610870/1‑45 (MQ=255)
                                    ttCGCCACCGTAGATAGCGTCTTCCTTGCGCATACGTTCTTtaat                                      >  1:328369/1‑45 (MQ=255)
                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:3184265/67‑1 (MQ=255)
                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:342751/67‑1 (MQ=255)
                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:3421052/67‑1 (MQ=255)
                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:961302/67‑1 (MQ=255)
                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:1752381/67‑1 (MQ=255)
                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:1049239/67‑1 (MQ=255)
                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:1164308/67‑1 (MQ=255)
                                    ttCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:1606605/67‑1 (MQ=255)
                                    ttCGCCACCATAGATGGCGCCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:890639/67‑1 (MQ=255)
                                     tCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACAtt             >  1:1620985/1‑69 (MQ=255)
                                     tCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACAtt             >  1:1387190/1‑69 (MQ=255)
                                     tCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACAtt             >  1:1900984/1‑69 (MQ=255)
                                          accaTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGttt                     <  1:52958/56‑1 (MQ=255)
                                          accaTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGttt                     <  1:1773336/56‑1 (MQ=255)
                                          accaTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGttt                     <  1:1164997/56‑1 (MQ=255)
                                          accaTAGATGGCGGCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGttt                     <  1:1969828/56‑1 (MQ=255)
                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGAc                <  1:2767097/58‑1 (MQ=255)
                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGt     >  1:2058057/1‑69 (MQ=255)
                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGt     >  1:1865226/1‑69 (MQ=255)
                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGt     >  1:644441/1‑69 (MQ=255)
                                             aTAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGt     >  1:343835/1‑69 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:695880/69‑1 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:766822/69‑1 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:772426/69‑1 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:920561/69‑1 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:2724464/69‑1 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:307092/69‑1 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:3435025/69‑1 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:205735/69‑1 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:2306310/69‑1 (MQ=255)
                                               aGATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:3200649/69‑1 (MQ=255)
                                                gATGGCGTCTTCCTTGGGCATACGTTCTTTAATAAAgg                                 >  1:2170510/1‑38 (MQ=255)
                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCg                               >  1:306757/1‑40 (MQ=255)
                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCTGTTTTCGa                 >  1:1147570/1‑54 (MQ=255)
                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 >  1:2183660/1‑54 (MQ=255)
                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGa                 >  1:801744/1‑54 (MQ=255)
                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGcc  >  1:1204799/1‑69 (MQ=255)
                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGcc  >  1:2874256/1‑69 (MQ=255)
                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGcc  >  1:3232131/1‑69 (MQ=255)
                                                gATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGcc  >  1:1093804/1‑69 (MQ=255)
                                                gATAGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGGGTGc   <  1:802731/68‑1 (MQ=255)
                                                   |                                                                 
GTAAGCGAAATCACGGAAGTAATGGTGGGCGCTCATTTCGCCACCATAGATGGCGTCTTCCTTGCGCATACGTTCTTTAATAAAGGCGTGTCCGGTTTTCGACATTACCGGCGTGCC  >  minE/1330566‑1330682

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: