Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 3,259,237 G→A 10.8% S613F (TCT→TTT)  tdcE ← pyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090483,259,2370GA10.8% 189.5 / 11.4 65S613F (TCT→TTT) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (28/29);  new base A (4/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

CGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCA  >  AP009048/3259170‑3259306
                                                                   |                                                                     
cgcgAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAAAAg                                                                    >  1:2145069/1‑71 (MQ=255)
cgcgAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAAAAg                                                                    >  1:1938242/1‑71 (MQ=255)
     aCGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGa                                                                 >  1:808313/1‑69 (MQ=255)
      cGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATa                                                               <  1:2924262/70‑1 (MQ=255)
      cGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAAAAGTGATa                                                               <  1:3036283/70‑1 (MQ=255)
        gtgtTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGt                                                             <  1:3034022/70‑1 (MQ=255)
           ttCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAg                                                          >  1:1347236/1‑70 (MQ=255)
              cGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGaa                                                                     >  1:3192118/1‑56 (MQ=255)
              cGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATc                                                      <  1:1912871/71‑1 (MQ=255)
              cGGCGCGAAGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATc                                                      >  1:66305/1‑71 (MQ=255)
               ggCGCTACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCg                                                     <  1:807962/71‑1 (MQ=255)
               ggCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTaaaa                                                                     >  1:1481/1‑55 (MQ=255)
               ggCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAg                                                                    <  1:847858/56‑1 (MQ=255)
               ggCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAg                                                                    <  1:2969618/56‑1 (MQ=255)
               ggCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAAAAGTGATAGTCAGAATc                                                      >  1:1566787/1‑70 (MQ=255)
               agcgcGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCg                                                     <  1:632086/70‑1 (MQ=255)
                gcgcGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGa                                                                 <  1:3432862/58‑1 (MQ=255)
                  gcgACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGa                                                    <  1:561125/69‑1 (MQ=255)
                   cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCg                                                     <  1:3326378/67‑1 (MQ=255)
                     acgacCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTggg                                               <  1:2312995/71‑1 (MQ=255)
                      cgacCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGt                                              <  1:347907/71‑1 (MQ=255)
                       gacCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTa                                             <  1:74401/71‑1 (MQ=255)
                         ccGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                                           >  1:1573275/1‑71 (MQ=255)
                          cGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                                           >  1:4564942/1‑70 (MQ=255)
                            tCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGAc                                        >  1:3334604/1‑71 (MQ=255)
                             ccGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGa                                         >  1:3545860/1‑69 (MQ=255)
                              cGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATc                                                      >  1:585849/1‑55 (MQ=255)
                              cGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTggg                                               >  1:980575/1‑62 (MQ=255)
                              cGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACgg                                      >  1:169880/1‑71 (MQ=255)
                                 cGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAggg                                          <  1:2377998/64‑1 (MQ=255)
                                 cGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                                           >  1:131499/1‑63 (MQ=255)
                                 cGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                                           >  1:154810/1‑63 (MQ=255)
                                  gTATTACCGGTTTTCTGGCCGTACACCACGTTAGaa                                                                     <  1:2126616/36‑1 (MQ=255)
                                     ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGAc                                                   <  1:3644380/51‑1 (MQ=255)
                                     ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCg                               <  1:155517/71‑1 (MQ=255)
                                       aCCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGaa                                                        >  1:444483/1‑44 (MQ=255)
                                        ccGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATa                            <  1:184668/71‑1 (MQ=255)
                                        ccGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGCTAGGGACGGCGTTGCGATa                            <  1:3841607/71‑1 (MQ=255)
                                        ccGGTTTTCTGGCCGTACACCACGTTAAAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGAt                             <  1:2541551/70‑1 (MQ=255)
                                          ggTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCg                                    >  1:1419061/1‑61 (MQ=255)
                                           gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGGCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGt                         <  1:611177/71‑1 (MQ=255)
                                            ttttCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGtt                        <  1:711161/71‑1 (MQ=255)
                                              ttCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTgg                      >  1:2725093/1‑71 (MQ=255)
                                               tCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTgg                      >  1:2058033/1‑70 (MQ=255)
                                               tCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTgg                      <  1:2317649/70‑1 (MQ=255)
                                                cTGGCCTTACACCACGTTATAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGtt                        <  1:2717938/67‑1 (MQ=255)
                                                   gCCGTACACCACGTTAGAAGTCATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCa                    <  1:4189793/68‑1 (MQ=255)
                                                    ccGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAgg                          <  1:3033317/61‑1 (MQ=255)
                                                    ccGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAgcgc                >  1:3315082/1‑71 (MQ=255)
                                                     cGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAg                           <  1:1136754/59‑1 (MQ=255)
                                                      gTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCtt              >  1:744044/1‑71 (MQ=255)
                                                       tACACCACGTTAGAAGTGATAGTCAGAATCGACTg                                                 >  1:3443269/1‑35 (MQ=255)
                                                         caccacGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGAt                             >  1:693623/1‑53 (MQ=255)
                                                         caccacGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAg                   >  1:3065142/1‑63 (MQ=255)
                                                         caccacGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTaa           <  1:1828131/71‑1 (MQ=255)
                                                         caccacGTTAAAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCt               >  1:4113023/1‑67 (MQ=255)
                                                           ccacGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAAtt         <  1:3829206/71‑1 (MQ=255)
                                                             acGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGtt                        <  1:3075266/54‑1 (MQ=255)
                                                             acGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAAttt        >  1:3269375/1‑70 (MQ=255)
                                                              cGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTg                                 >  1:3932308/1‑44 (MQ=255)
                                                              cGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATa                            >  1:3773618/1‑49 (MQ=255)
                                                              cGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATa                            >  1:3452584/1‑49 (MQ=255)
                                                              cGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTttc      >  1:4602791/1‑71 (MQ=255)
                                                               gTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTttct     <  1:719681/71‑1 (MQ=255)
                                                                ttAAAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTttctt    <  1:1295317/71‑1 (MQ=255)
                                                                  aGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCa  >  1:66667/1‑71 (MQ=255)
                                                                   |                                                                     
CGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGGTTGGCAGCGCTTTAATTTTCTTCA  >  AP009048/3259170‑3259306

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: