Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 608,280 A→G 17.4% A658A (GCT→GCC) ‡ pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported608,2800AG17.4% 235.4 / 39.2 109A658A (GCT→GCC) ‡pflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (46/44);  new base G (16/3);  total (62/47)
Fisher's exact test for biased strand distribution p-value = 1.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 112 total aligned reads displayed.

ACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGT  >  Exported/608212‑608350
                                                                    |                                                                      
aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg                                                                       >  1:4495514/1‑70 (MQ=255)
 cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:4424500/71‑1 (MQ=255)
 cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:3577571/71‑1 (MQ=255)
 cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:3063704/71‑1 (MQ=255)
 cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:3001695/71‑1 (MQ=255)
 cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:2603901/71‑1 (MQ=255)
 cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:601071/71‑1 (MQ=255)
    cgtcgtCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                                   <  1:747403/70‑1 (MQ=255)
     gtcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACg                                                                 >  1:4515711/1‑71 (MQ=255)
      tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACgg                                                                >  1:384793/1‑71 (MQ=255)
      tcgtcTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTATAGCGTAAGCAAACgg                                                                <  1:1381964/71‑1 (MQ=255)
             tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                                   >  1:2424067/1‑61 (MQ=255)
             tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                                   >  1:2353965/1‑61 (MQ=255)
             tACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:350380/69‑1 (MQ=255)
              aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGt                                                            <  1:1894733/67‑1 (MQ=255)
              aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                         <  1:4161336/70‑1 (MQ=255)
              aCCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                         <  1:5430160/70‑1 (MQ=255)
                   gTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                   >  1:5103641/1‑71 (MQ=255)
                    tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                  <  1:666583/71‑1 (MQ=255)
                    tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                   >  1:3162435/1‑70 (MQ=255)
                    tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                                   >  1:2595823/1‑70 (MQ=255)
                      cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                                  >  1:726501/1‑69 (MQ=255)
                           ggAACGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                            >  1:3755348/1‑70 (MQ=255)
                           ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           <  1:5653632/71‑1 (MQ=255)
                           ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           <  1:1390292/71‑1 (MQ=255)
                             aaCGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                            >  1:4607131/1‑68 (MQ=255)
                              aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCa                                                              >  1:3865439/1‑49 (MQ=255)
                              aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAg                                                             >  1:256866/1‑50 (MQ=255)
                              aCGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                          >  1:513774/1‑67 (MQ=255)
                              aCGATTGAGAAGGTGTACGAGATCCCATATTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCaat                                          >  1:4234740/1‑67 (MQ=255)
                               cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                                             >  1:2378580/1‑49 (MQ=255)
                               cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                                             >  1:4597779/1‑49 (MQ=255)
                               cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                                             >  1:4943854/1‑49 (MQ=255)
                                gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                <  1:5255274/61‑1 (MQ=255)
                                gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                                <  1:4060079/61‑1 (MQ=255)
                                 aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     <  1:1327759/71‑1 (MQ=255)
                                  ttgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                                     >  1:5678493/3‑70 (MQ=255)
                                  tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                                     <  1:2004201/38‑1 (MQ=255)
                                  tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     <  1:4135293/70‑1 (MQ=255)
                                  tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     >  1:1762916/1‑70 (MQ=255)
                                  tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     >  1:3348759/1‑70 (MQ=255)
                                   agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                   >  1:5573347/1‑71 (MQ=255)
                                   agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                   >  1:4378990/1‑71 (MQ=255)
                                   agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                   >  1:2969125/1‑71 (MQ=255)
                                    gagaAGGTGTAGGAGATACCATCTTTAGCGTAAGCaaa                                                                   >  1:2964523/1‑38 (MQ=255)
                                    gagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGc                                                       <  1:4579883/50‑1 (MQ=255)
                                      gaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                                <  1:920718/71‑1 (MQ=255)
                                       aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                               <  1:3033685/71‑1 (MQ=255)
                                       aaGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGc                               <  1:3026919/71‑1 (MQ=255)
                                         ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                    >  1:3370390/1‑48 (MQ=255)
                                         ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                    >  1:1242477/1‑48 (MQ=255)
                                         ggTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAAc                                                    >  1:1318659/1‑48 (MQ=255)
                                          gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                                     <  1:15115/62‑1 (MQ=255)
                                          gtgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCAcc                            >  1:397857/1‑71 (MQ=255)
                                           tgtAGTAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:4083921/39‑1 (MQ=255)
                                           tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                            <  1:5068583/54‑1 (MQ=255)
                                           tgtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCa                                            <  1:2288181/54‑1 (MQ=255)
                                            gtAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCtt                          >  1:267150/1‑71 (MQ=255)
                                             tAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                          >  1:324205/1‑54 (MQ=255)
                                              aGGAGATACCATCTTAAGCGAAAGCAAACGGCAGTTTAGc                                                       <  1:177471/40‑1 (MQ=38)
                                               ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:3020729/35‑1 (MQ=255)
                                               ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:5603934/35‑1 (MQ=255)
                                               ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:1787219/35‑1 (MQ=255)
                                               ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:3027479/35‑1 (MQ=255)
                                               ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                           <  1:1009061/35‑1 (MQ=255)
                                               ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTa                                                         >  1:1595621/1‑37 (MQ=255)
                                               ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCagaga                                        >  1:5311233/1‑54 (MQ=255)
                                               ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                       >  1:415944/1‑71 (MQ=255)
                                               cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGttt                                                          >  1:3768728/2‑36 (MQ=12)
                                               cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGttt                                                          >  1:1605421/2‑36 (MQ=12)
                                               cgagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                            <  1:813430/49‑1 (MQ=255)
                                                gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                      >  1:4114886/1‑68 (MQ=255)
                                                gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                      >  1:1693457/1‑68 (MQ=255)
                                                gagaTCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCtttgcg                      <  1:834402/71‑4 (MQ=255)
                                                gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                                   >  1:3596066/1‑58 (MQ=255)
                                                   aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                              <  1:3959556/44‑1 (MQ=255)
                                                     aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                     >  1:2515299/1‑35 (MQ=255)
                                                     aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                     >  1:2090346/1‑35 (MQ=255)
                                                     aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                                     >  1:498555/1‑35 (MQ=255)
                                                     aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                 >  1:5691076/1‑71 (MQ=255)
                                                      ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                     <  1:1633976/66‑1 (MQ=255)
                                                       cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa               >  1:489970/1‑71 (MQ=255)
                                                       cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa               >  1:196632/1‑71 (MQ=255)
                                                       cATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                          <  1:3901048/44‑1 (MQ=255)
                                                        aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc              <  1:1258203/71‑1 (MQ=255)
                                                          cTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAcca            >  1:5093733/1‑70 (MQ=255)
                                                           tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                          <  1:1066885/40‑1 (MQ=255)
                                                           tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                          >  1:4230043/1‑40 (MQ=255)
                                                           tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                           >  1:3917764/1‑39 (MQ=255)
                                                            ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt           >  1:3836571/1‑70 (MQ=255)
                                                            ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg          >  1:1192396/1‑71 (MQ=255)
                                                            ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCCTg          >  1:4201025/1‑71 (MQ=255)
                                                             tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc             <  1:635675/67‑1 (MQ=255)
                                                             tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt           <  1:2513253/69‑1 (MQ=255)
                                                             tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg          >  1:3385579/1‑70 (MQ=255)
                                                                cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc      <  1:3595044/71‑1 (MQ=255)
                                                                 gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg     >  1:5471197/1‑71 (MQ=255)
                                                                    aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc              >  1:3140718/1‑59 (MQ=255)
                                                                    aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc              >  1:4945308/1‑59 (MQ=255)
                                                                    aGCAAACGGCAGTTTAGAAACGGAAGTCAGAGAGGCTACTCCACCTTTCTGGTAACGACCGTGCATCGGGt  >  1:4693906/1‑71 (MQ=255)
                                                                    |                                                                      
ACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGT  >  Exported/608212‑608350

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: