Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 608,943 T→A 29.2% A437A (GCA→GCT)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported608,9430TA29.2% 233.7 / 116.5 144A437A (GCA→GCT) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (50/52);  new base A (24/18);  total (74/70)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

Only 100 of 146 total aligned reads displayed.

CAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGT  >  Exported/608876‑609010
                                                                   |                                                                   
cgaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcg                                                                  >  1:5607572/3‑70 (MQ=255)
cAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCa                                                                  >  1:2086234/1‑71 (MQ=255)
 gaCCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgc                                                                 >  1:4024147/2‑69 (MQ=255)
  aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                                <  1:2942887/71‑1 (MQ=255)
  aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                                <  1:1813025/71‑1 (MQ=255)
  aCCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcg                                                                <  1:4610772/71‑1 (MQ=255)
    cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                                              <  1:3899536/71‑1 (MQ=255)
    cTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCAcgcg                                                              >  1:2418760/1‑71 (MQ=255)
    cTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcgtg                                                              <  1:604767/71‑3 (MQ=255)
      gCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCa                                                            >  1:763935/1‑71 (MQ=255)
       cATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAc                                                           >  1:664257/1‑71 (MQ=255)
         ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc                                                          >  1:3300909/1‑70 (MQ=255)
         ttttCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCAcc                                                          >  1:3640368/1‑70 (MQ=255)
            tCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCg                                                         >  1:4919861/1‑68 (MQ=255)
             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                      <  1:3217052/70‑1 (MQ=255)
             cAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                      <  1:5057062/70‑1 (MQ=255)
                 ttttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                 >  1:4708645/1‑71 (MQ=255)
                 ttctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTg                                                 >  1:5586913/4‑71 (MQ=255)
                  tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                 >  1:2813492/1‑70 (MQ=255)
                  tttCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTg                                                 >  1:3031028/1‑70 (MQ=255)
                   ctCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGc                                                <  1:5301396/69‑1 (MQ=255)
                    tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                              >  1:2057723/1‑71 (MQ=255)
                    tCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                              >  1:1163239/1‑71 (MQ=255)
                    tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                    >  1:3078734/1‑65 (MQ=255)
                    tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                    >  1:4648764/1‑65 (MQ=255)
                    tCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCa                                               <  1:2517823/70‑1 (MQ=255)
                     cGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCAt                                              <  1:3273245/70‑1 (MQ=255)
                         aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                         <  1:898099/71‑1 (MQ=255)
                         aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                         >  1:2669789/1‑71 (MQ=255)
                         aaCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                         >  1:2402722/1‑71 (MQ=255)
                          aCGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGtt                                         >  1:3459613/1‑70 (MQ=255)
                           ccccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                               <  1:1155144/61‑1 (MQ=255)
                            cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                      >  1:3756992/2‑71 (MQ=255)
                            cccgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAc                                      >  1:1557742/2‑71 (MQ=255)
                             ccgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTAc                                      >  1:5557805/1‑70 (MQ=255)
                              cgccgTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaag                                                      <  1:4665827/53‑1 (MQ=255)
                              cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTgct                                         >  1:1081659/1‑64 (MQ=255)
                              cgccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                     >  1:4221099/1‑70 (MQ=255)
                                ccgTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaat                                  <  1:2375314/71‑2 (MQ=255)
                                    taaTTGCTTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCa                                               <  1:438023/52‑1 (MQ=255)
                                    tGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaa                                                       >  1:3720903/1‑46 (MQ=255)
                                     gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                              <  1:4810051/70‑1 (MQ=255)
                                     gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                              <  1:4320014/70‑1 (MQ=255)
                                     gATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                              <  1:2818379/70‑1 (MQ=255)
                                     aaTTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAt                                              >  1:4177436/2‑54 (MQ=255)
                                      aTTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                             <  1:3114830/70‑1 (MQ=255)
                                       ttGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaag                                                      <  1:3908913/44‑1 (MQ=12)
                                        tGCGTACAGCATGGTTTTCCCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTa                                       >  1:3366110/1‑58 (MQ=255)
                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                                         >  1:4149446/1‑38 (MQ=9)
                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCAcca                                                         >  1:1466690/1‑38 (MQ=9)
                                         gCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCaaagaa                                                    >  1:4326718/1‑44 (MQ=12)
                                         gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                         >  1:5340523/1‑71 (MQ=255)
                                         gCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                         >  1:3865652/1‑71 (MQ=255)
                                          catACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                    >  1:2913451/3‑43 (MQ=255)
                                          cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCgaagaa                                                    >  1:2067710/1‑43 (MQ=255)
                                          cGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCgg                         <  1:493001/70‑1 (MQ=255)
                                            tACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                      <  1:4018357/71‑1 (MQ=255)
                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                   >  1:136434/1‑40 (MQ=255)
                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                   >  1:5458941/1‑40 (MQ=255)
                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAAc                                                   >  1:3498844/1‑40 (MQ=255)
                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                    >  1:5663115/1‑71 (MQ=255)
                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                    >  1:5572364/1‑71 (MQ=255)
                                              cAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                    >  1:2240311/1‑71 (MQ=255)
                                               aGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAc                   <  1:1128710/71‑1 (MQ=255)
                                               aGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTAc                   <  1:196137/71‑1 (MQ=255)
                                                 cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTa                    <  1:5202882/68‑1 (MQ=255)
                                                 cATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACg                  >  1:2583007/1‑70 (MQ=255)
                                                   tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcag               >  1:814231/1‑71 (MQ=255)
                                                   tGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACAATCATCGGGCTTACg                  <  1:508464/68‑1 (MQ=255)
                                                     gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                            <  1:2207997/56‑1 (MQ=255)
                                                     gTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc              >  1:1569553/1‑70 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                             <  1:576990/54‑1 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatca                             <  1:3517010/54‑1 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCtt                     <  1:1943784/62‑1 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACgcagc              >  1:4480111/1‑69 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa            <  1:4812571/71‑1 (MQ=255)
                                                      ttttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaa            <  1:5109958/71‑1 (MQ=255)
                                                       tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                              >  1:5134384/1‑52 (MQ=255)
                                                       tttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag           >  1:3055222/1‑71 (MQ=255)
                                                        ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag           >  1:2663018/1‑70 (MQ=255)
                                                        ttCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaag           >  1:2391923/1‑70 (MQ=255)
                                                        ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg           >  1:5184458/1‑68 (MQ=255)
                                                        ttCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACgcagcagg           >  1:5207214/1‑68 (MQ=255)
                                                         tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                      <  1:5587630/58‑1 (MQ=255)
                                                         tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc          <  1:2032018/70‑1 (MQ=255)
                                                         tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc          <  1:5547664/70‑1 (MQ=255)
                                                         tCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc          <  1:5274704/70‑1 (MQ=255)
                                                           gCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAgcaagcaa        <  1:4224293/70‑1 (MQ=255)
                                                             cAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaatc                                 <  1:1905524/43‑3 (MQ=19)
                                                              aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc    >  1:2480333/1‑71 (MQ=255)
                                                              aGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGc    >  1:4947648/1‑71 (MQ=255)
                                                              aGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGc    <  1:4923163/71‑1 (MQ=255)
                                                               ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatc                              >  1:1341037/1‑44 (MQ=255)
                                                               ggTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCagcaagc          >  1:4642752/1‑64 (MQ=255)
                                                               ggTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCaa                                   >  1:1476382/1‑39 (MQ=16)
                                                                gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                      <  1:4766030/51‑1 (MQ=255)
                                                                gTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCt                      <  1:4772410/51‑1 (MQ=255)
                                                                gTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGt  <  1:2925223/71‑1 (MQ=255)
                                                                 tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAAc                                  <  1:1089070/38‑1 (MQ=255)
                                                                 tttGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGatcat                            <  1:1834487/44‑1 (MQ=255)
                                                                   |                                                                   
CAACCTGCATTTTCAGTTTTTCGTCAACGCCGCCGTTGATTGCGTACAGCATGGTTTTCGCCAGGTTTGCACGCGCACCGAAGAACTGCATTTGTTTACCAACGATCATCGGGCTTACGCAGCAAGCAATAGCGT  >  Exported/608876‑609010

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