Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 609,073 T→A 13.6% D394V (GAC→GTC)  pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported609,0730TA13.6% 202.0 / 22.9 81D394V (GAC→GTC) pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base A (7/4);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 3.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.60e-01

ATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTT  >  Exported/609004‑609143
                                                                     |                                                                      
aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:1082422/71‑1 (MQ=255)
aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:5588921/71‑1 (MQ=255)
aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       >  1:2302400/1‑71 (MQ=255)
aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:1551045/71‑1 (MQ=255)
 tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCg                                                                      <  1:5383566/71‑1 (MQ=255)
   gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt                                                                    <  1:2751295/71‑1 (MQ=255)
   gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa                                                                     <  1:940828/70‑1 (MQ=255)
   gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGAt                                                                    >  1:2137280/1‑71 (MQ=255)
      tAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGa                                                                 <  1:3301044/71‑1 (MQ=255)
         tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                                              <  1:1115463/68‑1 (MQ=255)
         tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                                 >  1:1970741/4‑68 (MQ=255)
         tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                                 >  1:2152665/4‑68 (MQ=255)
         tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                                                 >  1:3663707/4‑68 (MQ=255)
           aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:3179307/1‑71 (MQ=255)
           aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:445539/1‑71 (MQ=255)
           aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:705335/1‑71 (MQ=255)
           aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:3379225/1‑71 (MQ=255)
            tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:3116851/1‑70 (MQ=255)
            tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:3110437/1‑70 (MQ=255)
             cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:2781174/1‑69 (MQ=255)
              gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                                                          <  1:845802/70‑1 (MQ=255)
              gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAg                                                         >  1:5548385/1‑71 (MQ=255)
                tgttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                                       >  1:456444/1‑71 (MQ=255)
                         ccGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAt                                                 <  1:1150458/68‑1 (MQ=255)
                          cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                              <  1:2198931/70‑1 (MQ=255)
                          cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                                              <  1:5337889/70‑1 (MQ=255)
                            gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                           <  1:3077275/71‑1 (MQ=255)
                            gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                           <  1:5231204/71‑1 (MQ=255)
                            gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGAAACTTTAGCGGCGAATTTCtt                                            >  1:249047/1‑70 (MQ=255)
                             aCGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                                          <  1:3909367/71‑1 (MQ=255)
                              cGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                           >  1:1036772/1‑69 (MQ=255)
                              cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                                              >  1:1559307/1‑50 (MQ=12)
                              cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTg                                           >  1:5175063/1‑69 (MQ=255)
                              cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGa                                          <  1:1731496/70‑1 (MQ=255)
                                cATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                                       >  1:1500114/1‑55 (MQ=255)
                                cATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGAAgg                                       <  1:1224327/71‑2 (MQ=255)
                                 aTCAGATCGTCGTTTTCATACTGCAAAGAGGATGTGACGATAGACACCTGCGCGGCATATTTTTTGaa                                         >  1:3717520/1‑68 (MQ=255)
                                  tCAGTTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                     <  1:3847743/71‑1 (MQ=255)
                                  tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:611434/37‑1 (MQ=255)
                                  tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:16782/37‑1 (MQ=255)
                                  tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                           <  1:2152486/65‑1 (MQ=255)
                                  tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                     <  1:4737692/71‑1 (MQ=255)
                                  tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                     <  1:4503201/71‑1 (MQ=255)
                                  tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                     <  1:2073074/71‑1 (MQ=255)
                                   cAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGtt                                      >  1:4818782/1‑69 (MQ=255)
                                    aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                                       <  1:3710638/35‑1 (MQ=255)
                                    aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGa                                                                     >  1:5544640/1‑37 (MQ=255)
                                    aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                                         >  1:2111126/1‑65 (MQ=255)
                                    aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                                         >  1:1385118/1‑65 (MQ=255)
                                    aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                                          >  1:2873580/1‑64 (MQ=255)
                                      gTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                                                          <  1:5636241/46‑1 (MQ=255)
                                       tCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                                <  1:4179864/71‑1 (MQ=255)
                                        cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                                            >  1:3452629/1‑42 (MQ=255)
                                        cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                                >  1:5694392/1‑70 (MQ=255)
                                         aTCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGc                               <  1:5530304/70‑1 (MQ=255)
                                          tCGTTCTCATTCTGCAGAGAAGAGGTGTCGATGGACACTTTAGGGGCGAATttgttg                                           <  1:3695301/57‑1 (MQ=255)
                                          tCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAg                             <  1:1416829/71‑1 (MQ=255)
                                          tCGTTCTCATACTGCAGAGAAGAGCTGTCGATGGa                                                                 <  1:4657500/35‑1 (MQ=255)
                                            gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGGGGCAGtt                           >  1:145938/1‑71 (MQ=255)
                                            gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtt                           >  1:3218500/1‑71 (MQ=255)
                                             ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt                          <  1:701365/71‑1 (MQ=255)
                                             ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt                          <  1:5655573/71‑1 (MQ=255)
                                               ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                  >  1:2721395/1‑45 (MQ=255)
                                               ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                                                  >  1:3134255/1‑45 (MQ=255)
                                               ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                         <  1:5499503/70‑1 (MQ=255)
                                               ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                         >  1:1640732/1‑70 (MQ=255)
                                               ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt                         <  1:986105/70‑1 (MQ=255)
                                                tcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                                           <  1:2005458/51‑1 (MQ=255)
                                                tcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc                                     <  1:2449680/57‑1 (MQ=255)
                                                 cATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                        >  1:3973005/1‑69 (MQ=255)
                                                 cATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAATTCTTGAAGTTCAGCGGCAGttttt                        >  1:3973945/1‑69 (MQ=255)
                                                  aTACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                        <  1:2829398/68‑1 (MQ=255)
                                                  aTACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttttt                        <  1:4208587/68‑1 (MQ=255)
                                                        cAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGt                                       >  1:3659735/1‑47 (MQ=255)
                                                        cAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACca                 >  1:4997603/1‑69 (MQ=255)
                                                         agagAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAg              <  1:2205120/71‑1 (MQ=255)
                                                          gagaAGAGGTGTCGATGGGCACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGa             <  1:2248125/71‑1 (MQ=255)
                                                          gagaAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCa                              >  1:439162/1‑54 (MQ=255)
                                                          gagaAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCa                              >  1:3008616/1‑54 (MQ=255)
                                                           agaagaGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAg              <  1:1388042/69‑1 (MQ=255)
                                                              agagGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATgg         <  1:143769/71‑1 (MQ=255)
                                                                  gtgtCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGt        <  1:944340/68‑1 (MQ=255)
                                                                   tgtCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg                                <  1:5172726/43‑1 (MQ=255)
                                                                    gtCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGt   >  1:5575644/1‑71 (MQ=255)
                                                                     tCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGtt  >  1:3521916/1‑71 (MQ=255)
                                                                     tCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGt   <  1:2777022/70‑1 (MQ=255)
                                                                     |                                                                      
ATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTT  >  Exported/609004‑609143

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