Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 628,135 C→T 40.4% Y195Y (TAC→TAT kdsB → 3‑deoxy‑manno‑octulosonate cytidylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE628,1350CT40.4% 40.7 / 57.0 84Y195Y (TAC→TATkdsB3‑deoxy‑manno‑octulosonate cytidylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (32/18);  new base T (12/22);  total (44/40)
Fisher's exact test for biased strand distribution p-value = 1.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-01

CGTTGGCGATAACTTCCTGCGTCATCTTGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCGTGT  >  minE/628069‑628191
                                                                  |                                                        
cGTTGGCGATAACTTCCTGCGTCATCTTGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACg                                                         >  1:2742729/1‑68 (MQ=255)
                       aTCTTGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAg                                 <  1:1799739/69‑1 (MQ=255)
                       aTCTTGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAg                                 <  1:3300917/69‑1 (MQ=255)
                       aTCTTGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAg                                 <  1:2766910/69‑1 (MQ=255)
                       aTCTTGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAg                                 <  1:2828298/69‑1 (MQ=255)
                         cTTGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGaa                               <  1:3377016/69‑1 (MQ=255)
                         cTTGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGaa                               <  1:1312679/69‑1 (MQ=255)
                          ttGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGaa                               <  1:1244373/68‑1 (MQ=255)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:3308963/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:1880536/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:1880087/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:1742714/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:1708459/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:361037/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:1573196/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:636601/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:2565806/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:69766/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGtt                                                       >  1:2289208/1‑35 (MQ=37)
                                  tATGGCTACCGAGCAGGCTTTATCCGTCGTTATGt                                                        >  1:528515/1‑35 (MQ=38)
                                      gCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAAt                    >  1:663387/1‑67 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:2959569/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:3090586/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:2357298/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:953271/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:346333/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:348214/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:446810/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:530301/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:1438181/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:687548/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcagc         >  1:757834/1‑69 (MQ=255)
                                               cAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGagcag          >  1:239853/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtg   >  1:413389/1‑69 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtg   >  1:2617713/1‑69 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtg   >  1:1495765/1‑69 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtg   >  1:194395/1‑69 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtg   >  1:2000621/1‑69 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtg   >  1:2193526/1‑69 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:2436881/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:1523286/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:586039/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:1169669/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:2948880/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:2051863/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:73348/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:2414729/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:219530/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgt    >  1:2499301/1‑68 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCg     >  1:1334146/1‑67 (MQ=255)
                                                     ttATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGATTCgt    >  1:1731961/1‑68 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGATATGTTAGAGCAGCTTCgtgt  <  1:191890/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:1922496/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:1651411/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:455903/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:468467/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:516021/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:2715316/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:1368077/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:2579488/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:104196/67‑1 (MQ=255)
                                                        tCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCgtgt  <  1:2807491/67‑1 (MQ=255)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATTGAAATg                   <  1:3156/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:596453/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:1448877/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:639120/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:1364422/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:1165964/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:1075679/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:41270/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:375868/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:369146/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:1685770/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:174154/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:3408650/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:1826496/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:3103583/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:3076512/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:2998318/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:2964726/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:1974372/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:2647942/49‑1 (MQ=25)
                                                         ccGTCGTTATGTTACCTGGCAGCCAAGTCCGTTAGAACATATCGAAATg                   <  1:2591946/49‑1 (MQ=25)
                                                           gtcgtTATGTTACCTGTCAGCCAAGTCCGTTAGAACTTATCGAAATg                   <  1:2452928/47‑1 (MQ=255)
                                                                  |                                                        
CGTTGGCGATAACTTCCTGCGTCATCTTGGTATTTATGGCTACCGTGCAGGCTTTATCCGTCGTTACGTCAACTGGCAGCCAAGTCCGTTAGAACACATCGAAATGTTAGAGCAGCTTCGTGT  >  minE/628069‑628191

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: