Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 953,140 A→C 21.8% R279R (CGT→CGG)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb953,1400AC21.8% 71.4 / 35.6 55R279R (CGT→CGG) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (20/23);  new base C (7/5);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 5.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.60e-01

GCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCGAAGTAAGTCCACTGGA  >  W3110S.gb/953072‑953210
                                                                    |                                                                      
gcTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTAc                                                                       >  1:1090749/1‑70 (MQ=255)
  ttcttGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGAc                                                                    >  1:1345947/1‑71 (MQ=255)
     ttGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCg                                                                  >  1:1768399/1‑70 (MQ=255)
       gTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGa                                                                 <  1:602485/69‑1 (MQ=255)
       gTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGa                                                                 <  1:616097/69‑1 (MQ=255)
          cGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGTTGGAGGTCCGACCGAAGGa                                                             <  1:3022725/70‑1 (MQ=255)
           ggTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACa                                                           <  1:2182110/71‑1 (MQ=255)
              gATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGaa                                                                >  1:3757805/1‑63 (MQ=255)
               aTCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTCCGACCGAAGGACATTg                                                        >  1:728855/1‑70 (MQ=255)
                tCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTa                                                                        >  1:749179/1‑53 (MQ=255)
                tCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACa                                                           >  1:2374540/1‑66 (MQ=255)
                tCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTg                                                        >  1:231800/1‑69 (MQ=255)
                 cTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTCCGAcc                                                                   >  1:1526003/1‑57 (MQ=255)
                 cTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTCCGAcc                                                                   >  1:1531434/1‑57 (MQ=255)
                 cTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGAc                                                            >  1:1687832/1‑64 (MQ=255)
                  ttGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTCCGACCGAAGGACATTgcag                                                     >  1:519424/1‑70 (MQ=255)
                  ttGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGAc                                                            <  1:2797923/63‑1 (MQ=255)
                           ttCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCt                                            >  1:2555288/1‑70 (MQ=255)
                            tCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTCCGACCGAAGGACATTGCAGCACCg                                                <  1:865568/65‑1 (MQ=255)
                            tCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGAc                                                                    >  1:4466844/1‑45 (MQ=255)
                            tCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTgc                                                       >  1:3338942/1‑58 (MQ=255)
                             cAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTga                                          <  1:2700700/70‑1 (MQ=255)
                             cAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTga                                          <  1:1077653/70‑1 (MQ=255)
                              aGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTgaga                                        <  1:829700/71‑1 (MQ=255)
                                  cACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACa                                                           >  1:925248/1‑48 (MQ=255)
                                   aCGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTgcag                                                     >  1:248101/1‑53 (MQ=255)
                                   aCGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTaa                                   >  1:4456545/1‑71 (MQ=255)
                                      ttCGATGTACACATCCAGGAAGGTGGAGGTCCGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAAcag                                >  1:963034/1‑71 (MQ=255)
                                         gATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAAcagcag                             <  1:1138710/71‑1 (MQ=255)
                                         gATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAAcagcag                             <  1:2565572/71‑1 (MQ=255)
                                           tGTACACATCCAGGAAGGTGGAGGTCCGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGcc                           <  1:4447722/71‑1 (MQ=255)
                                           tGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGcc                           <  1:1126737/71‑1 (MQ=255)
                                              acacATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGcc                           <  1:1850419/68‑1 (MQ=255)
                                               cactTCCAGGAAGGTGGAGGTCCGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAAca                                 >  1:4259324/1‑61 (MQ=255)
                                               cacaTCCAGGAAGGTGGAGGTCCGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAgg                        <  1:2311723/70‑1 (MQ=255)
                                               cacaTCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAgg                        <  1:2841015/70‑1 (MQ=255)
                                               cacaTCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAgg                        <  1:4062777/70‑1 (MQ=255)
                                               cacaTCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAgg                        <  1:1492424/70‑1 (MQ=255)
                                               cacaTCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAgg                        >  1:2873634/1‑70 (MQ=255)
                                               cacaTCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGt                       <  1:2214406/71‑1 (MQ=255)
                                                acaTCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTa                      <  1:2979568/71‑1 (MQ=255)
                                                 caTCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTa                      <  1:3740201/70‑1 (MQ=255)
                                                  aTCCAGGAAGGTGGAGGTACGACCGAAGGACATTgc                                                       <  1:1467217/36‑1 (MQ=255)
                                                     cAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCg                  >  1:2012931/1‑70 (MQ=255)
                                                      aGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCTTTCTGAGACTTAACAGCAGCCAGGTAGCCg                  >  1:2017392/1‑69 (MQ=255)
                                                       ggAAGGTGGAGCTACGACCGAAGGACATTGCAGCACCGTTCTgaca                                        <  1:143804/46‑3 (MQ=255)
                                                        gAAGGTGGAGGTCCGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCGaag               >  1:2532820/1‑70 (MQ=255)
                                                        gAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACt                                      >  1:2647719/1‑47 (MQ=255)
                                                        gAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACt                                      >  1:4457040/1‑47 (MQ=255)
                                                         aaGGTGGAGGTCCGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCGaagt              <  1:4639287/70‑1 (MQ=255)
                                                         aaGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCGaagt              <  1:4604153/70‑1 (MQ=255)
                                                          aGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTa                      >  1:2469191/1‑61 (MQ=255)
                                                             tGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACagcagc                            <  1:1797633/52‑1 (MQ=255)
                                                             tGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCGAAGTAAGTc         <  1:4419476/71‑1 (MQ=255)
                                                               gAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCGAAGTAAGTCCa       >  1:3467331/1‑71 (MQ=255)
                                                                  gTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCGAAGTAAGTCCACTg    <  1:509914/71‑1 (MQ=255)
                                                                    aCGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCGAAGTAAGTCCACTGGa  >  1:2706323/1‑71 (MQ=255)
                                                                    |                                                                      
GCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGAAGGTGGAGGTACGACCGAAGGACATTGCAGCACCGTTCTGAGACTTAACAGCAGCCAGGTAGCCGAAGTAAGTCCACTGGA  >  W3110S.gb/953072‑953210

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: ‑