Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 609,336 A→G 18.2% V306V (GTT→GTC) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported609,3360AG18.2% 143.1 / 23.9 69V306V (GTT→GTC) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (33/23);  new base G (4/8);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

CAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGA  >  Exported/609266‑609406
                                                                      |                                                                      
cAGAGAACAGTTCATCGTATTCCTGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                                        <  1:5031927/71‑1 (MQ=255)
cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                                        <  1:278814/71‑1 (MQ=255)
  gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                                                       >  1:878281/1‑69 (MQ=255)
  gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                                                                       >  1:2931471/1‑69 (MQ=255)
    gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa                                                                    <  1:1885619/71‑1 (MQ=255)
      aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt                                                                  >  1:1719359/1‑71 (MQ=255)
      aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt                                                                  >  1:4089161/1‑71 (MQ=255)
       cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  >  1:2999037/1‑70 (MQ=255)
            cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCGACCATTTCCtac                                                            >  1:3684893/1‑69 (MQ=255)
            cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                                            >  1:3693472/1‑71 (MQ=255)
            cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc                                                            >  1:1981534/1‑71 (MQ=255)
               caaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgc                                                          >  1:1560493/4‑70 (MQ=255)
                gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGctt                                                        >  1:4107251/1‑71 (MQ=255)
                aaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcc                                                         <  1:939224/68‑2 (MQ=255)
                   ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGcttctt                                                     >  1:5283728/1‑71 (MQ=255)
                    tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  <  1:4064226/57‑1 (MQ=255)
                    tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  <  1:5318468/57‑1 (MQ=255)
                     ccGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCaa                                                                       >  1:1607554/1‑51 (MQ=255)
                      cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                  <  1:4578514/71‑1 (MQ=255)
                      cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                  <  1:236672/71‑1 (MQ=255)
                        gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                             >  1:4909761/3‑58 (MQ=255)
                        gtgTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTg                                                             >  1:1729511/3‑58 (MQ=255)
                        gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                                                                        <  1:3155037/47‑1 (MQ=255)
                        gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAttt                                                                 <  1:4530147/54‑1 (MQ=255)
                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:1205687/60‑3 (MQ=255)
                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:1301603/60‑3 (MQ=255)
                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:339874/60‑3 (MQ=255)
                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:3122764/60‑3 (MQ=255)
                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:4251922/60‑3 (MQ=255)
                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:4567928/60‑3 (MQ=255)
                           tACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcct                                                        <  1:2205378/60‑3 (MQ=255)
                             cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:1787213/1‑50 (MQ=255)
                             cGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:299515/1‑50 (MQ=255)
                                aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:1682521/1‑47 (MQ=255)
                                aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:2102882/1‑47 (MQ=255)
                                aGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:4362599/1‑47 (MQ=255)
                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  <  1:811417/44‑1 (MQ=255)
                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                                                                  <  1:242077/44‑1 (MQ=255)
                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                       >  1:299446/1‑71 (MQ=255)
                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                       <  1:3169799/71‑1 (MQ=255)
                                 ggAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                       >  1:1523225/1‑71 (MQ=255)
                                  gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCt                                         >  1:4128862/1‑68 (MQ=255)
                                  gAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                                      >  1:5601422/1‑71 (MQ=255)
                                    aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTg                                       >  1:5283690/1‑68 (MQ=255)
                                    aGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGc                                      >  1:2070553/1‑69 (MQ=255)
                                       gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATc                                          <  1:4481714/62‑1 (MQ=255)
                                       gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                                 >  1:990197/1‑71 (MQ=255)
                                       gAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCt                                 <  1:1970934/71‑1 (MQ=255)
                                          ccATACGCAGTTTCATGACCAGGTGGTCAACCATTTc                                                                >  1:3614193/1‑37 (MQ=255)
                                            aTACGCAGTTTCTTGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt                                                  <  1:30087/49‑1 (MQ=255)
                                               cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc                         <  1:4365751/71‑1 (MQ=255)
                                               cGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTc                         <  1:297722/71‑1 (MQ=255)
                                                   gTTTCATGACCAGGTGGTCTACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGt                     <  1:5111201/71‑1 (MQ=255)
                                                    tttCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCttt                               <  1:960013/60‑1 (MQ=255)
                                                      tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                                        >  1:4065374/1‑49 (MQ=255)
                                                      tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCtt                                        >  1:422164/1‑49 (MQ=255)
                                                      tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:1026504/70‑1 (MQ=255)
                                                      tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:1665705/70‑1 (MQ=255)
                                                      tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   >  1:728656/1‑70 (MQ=255)
                                                      tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:2058404/70‑1 (MQ=255)
                                                      tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:333195/70‑1 (MQ=255)
                                                      tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   >  1:3027650/1‑70 (MQ=255)
                                                      tCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTc                   <  1:4254341/70‑1 (MQ=255)
                                                         tGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATg               >  1:178792/1‑71 (MQ=255)
                                                             cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                             >  1:2169634/1‑53 (MQ=255)
                                                             cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg                            >  1:2332258/1‑54 (MQ=255)
                                                             cAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAg                            >  1:1121064/1‑54 (MQ=255)
                                                               ggtggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCa                             >  1:680557/1‑51 (MQ=255)
                                                                gtggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACAt        >  1:5601343/1‑71 (MQ=255)
                                                                  ggtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCa                                    >  1:2663644/1‑41 (MQ=255)
                                                                   gtCAACCATTTCCTGCTCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATcc      <  1:1769834/70‑1 (MQ=255)
                                                                   gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATc       >  1:4725418/1‑69 (MQ=255)
                                                                   gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATc       >  1:2417452/1‑69 (MQ=255)
                                                                   gtCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCa     <  1:4826755/71‑1 (MQ=255)
                                                                     cAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTa             <  1:5674027/61‑1 (MQ=255)
                                                                      aaCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAgg   >  1:4826953/1‑70 (MQ=255)
                                                                      aaCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAgg   <  1:3425056/70‑1 (MQ=255)
                                                                      aaCCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGa  <  1:1648544/71‑1 (MQ=255)
                                                                      aaCCATTTACTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGa  <  1:2019610/71‑1 (MQ=255)
                                                                      |                                                                      
CAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTTCGGTGATCTTGCCAGCTTTCAGGTCACGTTCGATGTACACATCCAGGA  >  Exported/609266‑609406

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: