Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 608,274 A→G 17.2% A660A (GCT→GCC pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported608,2740AG17.2% 260.0 / 41.6 116A660A (GCT→GCCpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (51/45);  new base G (16/4);  total (67/49)
Fisher's exact test for biased strand distribution p-value = 4.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 117 total aligned reads displayed.

TCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTG  >  Exported/608204‑608342
                                                                      |                                                                    
tCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTa                                                                      <  1:199760/71‑1 (MQ=255)
  ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAg                                                                     >  1:1839091/1‑70 (MQ=255)
  ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAg                                                                     >  1:1325124/1‑70 (MQ=255)
  ttACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAg                                                                     >  1:2734750/1‑70 (MQ=255)
   tACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                                   >  1:4288506/1‑71 (MQ=255)
    aCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                                   >  1:4504611/1‑70 (MQ=255)
    aCGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCg                                                                   >  1:1604319/1‑70 (MQ=255)
      gAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAg                                                                     <  1:4489631/66‑1 (MQ=255)
        aCTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAg                                                               >  1:4495514/1‑70 (MQ=255)
         cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                             <  1:3001695/71‑1 (MQ=255)
         cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                             <  1:3063704/71‑1 (MQ=255)
         cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                             <  1:2603901/71‑1 (MQ=255)
         cTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCa                                                             <  1:3577571/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:1411211/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:5430160/70‑1 (MQ=255)
                        ccAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGc                                                       <  1:4068465/62‑1 (MQ=255)
                            tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                           >  1:3162435/1‑70 (MQ=255)
                            tGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACg                                           >  1:2595823/1‑70 (MQ=255)
                              cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACgg                                          >  1:726501/1‑69 (MQ=255)
                              cGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                        <  1:184186/71‑1 (MQ=255)
                                ttCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTa                                                                 <  1:1140240/44‑1 (MQ=255)
                                   ggAACGATTGAGAAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCa                                    >  1:3755348/1‑70 (MQ=255)
                                   ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                   <  1:5653632/71‑1 (MQ=255)
                                   ggAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCag                                   <  1:1390292/71‑1 (MQ=255)
                                    gAACGATTGAGAAGGTGTACGAGATCCCATCTTTGGCg                                                                   <  1:2045071/38‑1 (MQ=21)
                                     aaCGATAGAGAATGTGTAGGAGATACCATCTTTAg                                                                     <  1:1369106/35‑1 (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:4597779/1‑49 (MQ=255)
                                       cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                                     >  1:4943854/1‑49 (MQ=255)
                                       cGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAg                                                     >  1:2378580/1‑49 (MQ=255)
                                        gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                        <  1:4060079/61‑1 (MQ=255)
                                        gATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGaa                                        <  1:5255274/61‑1 (MQ=255)
                                         aTAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                             <  1:1327759/71‑1 (MQ=255)
                                          ttgagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGc                             >  1:5678493/3‑70 (MQ=255)
                                          tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                             >  1:1762916/1‑70 (MQ=255)
                                          tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                             >  1:3348759/1‑70 (MQ=255)
                                          tAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGc                             <  1:4135293/70‑1 (MQ=255)
                                           agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                           >  1:2969125/1‑71 (MQ=255)
                                           agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                           >  1:4378990/1‑71 (MQ=255)
                                           agagAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                           >  1:5573347/1‑71 (MQ=255)
                                            gagaAGGTGTACGAGATCCCATCTTTGGCGTAGGTGAACGGCAGTTTCGc                                               <  1:4579883/50‑1 (MQ=255)
                                              gaAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTg                        <  1:920718/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:3013444/1‑48 (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)
                                                       ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                   <  1:3020729/35‑1 (MQ=255)
                                                       ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                   <  1:3027479/35‑1 (MQ=255)
                                                       ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                   <  1:1009061/35‑1 (MQ=255)
                                                       ggAGATACCATCTTTAGCGTAAGCAAACGGCAGtt                                                   <  1:1787219/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: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:4826911/1‑58 (MQ=255)
                                                        gagaTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTa                           >  1:3596066/1‑58 (MQ=255)
                                                           aTACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGt                                      <  1:3959556/44‑1 (MQ=255)
                                                             aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                             >  1:498555/1‑35 (MQ=255)
                                                             aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                             >  1:2090346/1‑35 (MQ=255)
                                                             aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCaa                                             >  1:2515299/1‑35 (MQ=255)
                                                             aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGa                                         >  1:585966/1‑39 (MQ=255)
                                                             aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc         >  1:5691076/1‑71 (MQ=255)
                                                              ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc         >  1:3321072/1‑70 (MQ=255)
                                                              ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc         >  1:3762863/1‑70 (MQ=255)
                                                              ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg             <  1:1633976/66‑1 (MQ=255)
                                                              ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc         >  1:526394/1‑70 (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)
                                                                 tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGATCACGAc      >  1:432341/1‑70 (MQ=255)
                                                                  cTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAcca    >  1:5093733/1‑70 (MQ=255)
                                                                   tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                  >  1:4230043/1‑40 (MQ=255)
                                                                   tttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCaga                                  <  1:1066885/40‑1 (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)
                                                                      |                                                                    
TCTTACGAACTTCGTCGTCTTTACCCAGTGCGTTCGGAACGATAGAGAAGGTGTAGGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTG  >  Exported/608204‑608342

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: