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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported608,7060GA26.3% 196.9 / 72.0 118D516D (GAC→GATpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (54/33);  new base A (13/18);  total (67/51)
Fisher's exact test for biased strand distribution p-value = 6.00e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

TTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAG  >  Exported/608636‑608776
                                                                       |                                                                      
ttGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACgg                                                                        <  1:709857/71‑1 (MQ=255)
 tGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGt                                                                       >  1:3881570/1‑71 (MQ=255)
  gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                                                                      >  1:4431275/1‑71 (MQ=255)
  gCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                                                                      >  1:5365467/1‑71 (MQ=255)
   cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTc                                                                      <  1:70206/70‑1 (MQ=255)
   cAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCg                                                                     >  1:3893114/1‑71 (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)
              tCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                          <  1:5488704/71‑1 (MQ=255)
                agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAg                                                                >  1:4966893/1‑63 (MQ=255)
                agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                          >  1:2339797/1‑69 (MQ=255)
                agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                          >  1:3859894/1‑69 (MQ=255)
                agcagcAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCa                                                        <  1:5218487/71‑1 (MQ=255)
                       aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagaga                                                    >  1:5597367/1‑68 (MQ=255)
                       aCGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGc                                                 >  1:225641/1‑71 (MQ=255)
                        ccgACAGGCCCGCGATGCCGCATGCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCt                                                >  1:2536424/3‑71 (MQ=255)
                           aCAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                             >  1:190803/1‑71 (MQ=255)
                                 cAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                               >  1:4797998/1‑63 (MQ=255)
                                 cAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCtt                                               >  1:449662/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:4739138/50‑1 (MQ=255)
                                   gCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                          <  1:1076884/50‑1 (MQ=255)
                                   gCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCAt                                                          <  1:5505039/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)
                                       tACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAgc                                                               >  1:350219/1‑41 (MQ=255)
                                         ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagag                                                     >  1:1411553/1‑49 (MQ=255)
                                         ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagag                                                     >  1:2049538/1‑49 (MQ=255)
                                         ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTg                                       <  1:4109552/63‑1 (MQ=255)
                                         ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                                <  1:1880320/70‑1 (MQ=255)
                                         ccACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGt                                <  1:4063157/70‑1 (MQ=255)
                                          cacaCGCCATCGTGCGGATAACGTCACGGTCGTGCAGCGCCATCaga                                                      >  1:4626461/1‑47 (MQ=255)
                                           acacGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagag                                                     >  1:5594075/1‑47 (MQ=255)
                                           acacGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCagag                                                     >  1:2945116/1‑47 (MQ=255)
                                            catgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg                            <  1:3557553/68‑1 (MQ=255)
                                            catgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTg                            <  1:4083297/68‑1 (MQ=255)
                                            cacGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTg                            >  1:3314966/1‑71 (MQ=255)
                                             acGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGc                           <  1:3358513/71‑1 (MQ=255)
                                              tgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATc                                                         <  1:965006/39‑1 (MQ=9)
                                              tgCCATAGTGCGATAGACATCACGATCGTGCAGCGCCATc                                                         <  1:1138975/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:1718547/1‑38 (MQ=255)
                                                            tAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCg                                             >  1:3158117/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)
                                                               acatCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt          <  1:4782496/68‑1 (MQ=255)
                                                               aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt          <  1:392572/71‑1 (MQ=255)
                                                                cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt          >  1:3245952/1‑70 (MQ=255)
                                                                 atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata            >  1:1286047/2‑66 (MQ=255)
                                                                 atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata            >  1:3689863/2‑66 (MQ=255)
                                                                 atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt          <  1:1976578/68‑1 (MQ=255)
                                                                 atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTca        <  1:3606150/70‑1 (MQ=255)
                                                                 atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTca        <  1:5568922/70‑1 (MQ=255)
                                                                 atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTca        <  1:2054779/70‑1 (MQ=255)
                                                                  tCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt          <  1:5486016/68‑1 (MQ=255)
                                                                  tCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTa                      <  1:4340230/56‑1 (MQ=255)
                                                                  tCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTa                      <  1:783134/56‑1 (MQ=255)
                                                                  tCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt          <  1:2958806/68‑1 (MQ=255)
                                                                   cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                 >  1:5182136/1‑60 (MQ=255)
                                                                   cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagt      >  1:3647712/1‑71 (MQ=255)
                                                                   cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgat             >  1:5504193/1‑64 (MQ=255)
                                                                   cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTc         <  1:2432898/68‑1 (MQ=255)
                                                                   cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTc         <  1:378611/68‑1 (MQ=255)
                                                                    aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg     >  1:2192726/1‑71 (MQ=255)
                                                                    aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg     <  1:2905936/71‑1 (MQ=255)
                                                                    aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg     >  1:3045651/1‑71 (MQ=255)
                                                                    aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCg                              >  1:2028101/1‑46 (MQ=255)
                                                                    aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg     >  1:3443817/1‑69 (MQ=255)
                                                                    aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg     >  1:956321/1‑69 (MQ=255)
                                                                    aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg     >  1:323897/1‑69 (MQ=255)
                                                                    aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg     >  1:2214593/1‑69 (MQ=255)
                                                                       gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag  >  1:2970207/1‑71 (MQ=255)
                                                                       gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag  >  1:3317045/1‑71 (MQ=255)
                                                                       gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag  >  1:3504231/1‑71 (MQ=255)
                                                                       |                                                                      
TTGCAGACAGGGAGTCAGCAGCAACGGACAGACCAGCGATACCACACGCCATGGTGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAG  >  Exported/608636‑608776

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: