Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 608,700 G→A 10.5% D518D (GAC→GAT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported608,7000GA10.5% 262.4 / 15.3 95D518D (GAC→GATpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (52/33);  new base A (4/6);  total (56/39)
Fisher's exact test for biased strand distribution p-value = 3.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

ATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTT  >  Exported/608630‑608768
                                                                       |                                                                    
atTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACg                                                                      >  1:1149481/1‑71 (MQ=255)
 tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                                                     >  1:4624953/1‑71 (MQ=255)
 tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                                                     >  1:3464094/1‑71 (MQ=255)
 tttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                                                     >  1:4120457/1‑71 (MQ=255)
  ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                                                     >  1:3549713/1‑70 (MQ=255)
  ttGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGt                                                                     >  1:4482926/1‑70 (MQ=255)
   tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTc                                                                    <  1:4495706/70‑1 (MQ=255)
   tGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACACCATGGTGCGGATAACGTCa                                                                   >  1:1826472/1‑71 (MQ=255)
      ttGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACgg                                                                <  1:709857/71‑1 (MQ=255)
       tGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGt                                                               >  1:3881570/1‑71 (MQ=255)
        gCAGACAGGGAGTCAGCAGCGACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                                                              >  1:1409918/1‑71 (MQ=255)
        gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                                                              >  1:5365467/1‑71 (MQ=255)
        gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                                                              >  1:4431275/1‑71 (MQ=255)
         cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                                                              <  1:70206/70‑1 (MQ=255)
         cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCg                                                             >  1:3893114/1‑71 (MQ=255)
           gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCgcg                                                           <  1:994016/71‑3 (MQ=255)
           gACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTg                                                           <  1:2664318/71‑1 (MQ=255)
            aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTg                                                           >  1:4380746/1‑70 (MQ=255)
            aCAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGc                                                          <  1:259841/71‑1 (MQ=255)
             cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGc                                                          <  1:5063168/70‑1 (MQ=255)
             cAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCa                                                         <  1:2870297/71‑1 (MQ=255)
              aGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAg                                                        >  1:4929491/1‑71 (MQ=255)
                ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCa                                                         >  1:767944/1‑68 (MQ=255)
                ggAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCa                                                         >  1:4325930/1‑68 (MQ=255)
                    tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                  <  1:5488704/71‑1 (MQ=255)
                      agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAg                                                        >  1:4966893/1‑63 (MQ=255)
                      agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                  >  1:3859894/1‑69 (MQ=255)
                      agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                  >  1:2339797/1‑69 (MQ=255)
                      agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCa                                                <  1:5218487/71‑1 (MQ=255)
                             aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagaga                                            >  1:606469/1‑68 (MQ=255)
                             aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagaga                                            >  1:5597367/1‑68 (MQ=255)
                             aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGc                                         >  1:225641/1‑71 (MQ=255)
                              ccgACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCt                                        >  1:2536424/3‑71 (MQ=255)
                               ggACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                       <  1:203083/71‑1 (MQ=255)
                                 aCAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                     >  1:190803/1‑71 (MQ=255)
                                       ccgCGATGCCGCATGCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTg                               >  1:2590246/3‑71 (MQ=255)
                                       cAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                       >  1:449662/1‑63 (MQ=255)
                                       cAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                       >  1:4797998/1‑63 (MQ=255)
                                       cAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                       >  1:4804471/1‑63 (MQ=255)
                                       cAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                       >  1:789345/1‑63 (MQ=255)
                                       cAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                               <  1:4899071/71‑1 (MQ=255)
                                       cAGCGATACCACACGCCATGGTGCGGATAACGTCACGGGCGTGCAGCGCCAt                                                  <  1:4551539/52‑1 (MQ=255)
                                         gCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                  <  1:5505039/50‑1 (MQ=255)
                                         gCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                  <  1:4739138/50‑1 (MQ=255)
                                         gCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                  <  1:774347/50‑1 (MQ=255)
                                         gCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                  <  1:1076884/50‑1 (MQ=255)
                                         gCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTa                             >  1:1757715/1‑71 (MQ=255)
                                          cGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                  <  1:3946373/49‑1 (MQ=255)
                                           gATACCACACGCCATGGTGCGGATAACGTCACGGTCGTg                                                           >  1:5614179/1‑39 (MQ=255)
                                           gATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCaga                                              >  1:2844041/1‑52 (MQ=255)
                                           gATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCa                                                >  1:761863/1‑50 (MQ=255)
                                           gATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTAc                            >  1:4218048/1‑70 (MQ=255)
                                           gATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACt                           >  1:4981778/1‑71 (MQ=255)
                                            aTACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTc                                      >  1:2933757/1‑59 (MQ=255)
                                             tACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAgc                                                       >  1:350219/1‑41 (MQ=255)
                                               ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagag                                             >  1:2049538/1‑49 (MQ=255)
                                               ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagag                                             >  1:1411553/1‑49 (MQ=255)
                                               ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                               <  1:1775790/63‑1 (MQ=255)
                                               ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                               <  1:4109552/63‑1 (MQ=255)
                                               ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                        <  1:1880320/70‑1 (MQ=255)
                                               ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                        <  1:4063157/70‑1 (MQ=255)
                                               ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTc                       <  1:4524482/71‑1 (MQ=255)
                                                cacaCGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                        <  1:3882643/69‑1 (MQ=255)
                                                cacaCGCCATCGTGCGGATAACGTCACGGTCGTGCAGCGCCATCaga                                              >  1:4626461/1‑47 (MQ=255)
                                                 acacGCCATGGTGCTGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGt                     <  1:5588762/71‑1 (MQ=255)
                                                 acacGCCATGGTGCGGATGACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGt                     <  1:5425087/71‑1 (MQ=255)
                                                 acacGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagag                                             >  1:2945116/1‑47 (MQ=255)
                                                 acacGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagag                                             >  1:5594075/1‑47 (MQ=255)
                                                  catgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg                    <  1:4083297/68‑1 (MQ=255)
                                                  catgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg                    <  1:3102444/68‑1 (MQ=255)
                                                  catgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg                    <  1:3557553/68‑1 (MQ=255)
                                                  cacGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg                    >  1:3314966/1‑71 (MQ=255)
                                                   acGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGc                   <  1:3358513/71‑1 (MQ=255)
                                                    tgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATc                                                 <  1:1138975/39‑1 (MQ=9)
                                                    tgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATc                                                 <  1:965006/39‑1 (MQ=9)
                                                    tgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATc                                                 <  1:126045/39‑1 (MQ=9)
                                                    tgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa                  >  1:3838314/2‑71 (MQ=255)
                                                    tgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCa                  >  1:568748/2‑71 (MQ=255)
                                                      ccATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATg                <  1:4580196/71‑1 (MQ=255)
                                                       cATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGt               <  1:4603465/71‑1 (MQ=255)
                                                        aTGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg                    <  1:3682641/65‑1 (MQ=255)
                                                        aTGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg                    <  1:5161873/65‑1 (MQ=255)
                                                           gTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCt                                        >  1:4085135/1‑42 (MQ=255)
                                                           gTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTc                                      >  1:4631848/1‑44 (MQ=255)
                                                           gTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTg           <  1:1283842/71‑1 (MQ=255)
                                                            tGCGGATAACGTCACGGTCGTGCAGCGCCATCagaga                                            >  1:3987840/1‑37 (MQ=255)
                                                               ggATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAg             >  1:697820/1‑65 (MQ=255)
                                                                 aTAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga      >  1:1649975/1‑70 (MQ=255)
                                                                  tAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                     >  1:3158117/1‑38 (MQ=255)
                                                                  tAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                     >  1:1718547/1‑38 (MQ=255)
                                                                  tAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg    >  1:2485028/1‑71 (MQ=255)
                                                                   aaCGTCACGGTCGTGCAGCGCCATCAGAGAGTCTTc                                      >  1:1685677/1‑36 (MQ=255)
                                                                   aaCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt   >  1:3993825/1‑71 (MQ=255)
                                                                   aaCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt   >  1:5158606/1‑71 (MQ=255)
                                                                     aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt  <  1:392572/71‑1 (MQ=255)
                                                                      cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt  >  1:3245952/1‑70 (MQ=255)
                                                                       |                                                                    
ATTTGATTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTT  >  Exported/608630‑608768

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: