Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 608,415 A→G 14.4% Y613Y (TAT→TAC pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported608,4150AG14.4% 201.6 / 22.3 83Y613Y (TAT→TACpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (32/39);  new base G (7/5);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCT  >  Exported/608347‑608485
                                                                    |                                                                      
gggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATa                                                                      <  1:1332168/71‑1 (MQ=255)
     aGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACaa                                                                 >  1:4628740/1‑71 (MQ=255)
       cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                                               <  1:3273572/71‑1 (MQ=255)
       cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                                                               >  1:5220220/1‑71 (MQ=255)
        aCCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                              <  1:1479695/71‑1 (MQ=255)
         cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                              >  1:2916497/1‑70 (MQ=255)
           cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                           <  1:1196942/71‑1 (MQ=255)
           cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                           <  1:4386550/71‑1 (MQ=255)
                cGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                                       >  1:3633108/1‑70 (MQ=255)
                cGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                                      <  1:899285/71‑1 (MQ=255)
                  aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                           <  1:2862864/64‑1 (MQ=255)
                  aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                           <  1:5010659/64‑1 (MQ=255)
                  aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGAt                                                    >  1:3351370/1‑71 (MQ=255)
                    cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATgg                                                  >  1:3361145/1‑71 (MQ=255)
                      gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                                 >  1:1669343/1‑70 (MQ=255)
                      gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                                 >  1:2058523/1‑70 (MQ=255)
                      gcgcgcCAGCACGACGACCGTCTGGGGTGTTAACCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt                                                 >  1:2981590/1‑70 (MQ=255)
                        gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                                                                    <  1:2340756/49‑1 (MQ=255)
                        gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                                                                    <  1:3983295/49‑1 (MQ=255)
                        gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                              <  1:210484/55‑1 (MQ=255)
                        gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                                                              <  1:4115035/55‑1 (MQ=255)
                         cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg                                              >  1:2490241/1‑70 (MQ=255)
                          gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                            >  1:781163/1‑71 (MQ=255)
                              gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                                                   <  1:1390336/44‑1 (MQ=255)
                              gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg                                                      <  1:392881/57‑1 (MQ=255)
                              gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                        >  1:2419006/1‑71 (MQ=255)
                              gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                        >  1:1605881/1‑71 (MQ=255)
                              gCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                        >  1:1599268/1‑71 (MQ=255)
                               cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATacac                                                                   <  1:3516012/43‑1 (MQ=255)
                               cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCa                                               >  1:2961006/1‑63 (MQ=255)
                               cACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGAc                                       <  1:5056468/71‑1 (MQ=255)
                                 cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                                     <  1:1131484/71‑1 (MQ=255)
                                 cgacgaCCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTg                                     <  1:985849/71‑1 (MQ=255)
                                  gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTATaa                                                         >  1:964193/1‑50 (MQ=255)
                                  gacgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACt                                      >  1:779657/1‑69 (MQ=255)
                                   acgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                           <  1:5523749/47‑1 (MQ=255)
                                    cgacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGt                                  >  1:1411661/1‑71 (MQ=255)
                                     gacCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTc                                 <  1:2463970/71‑1 (MQ=255)
                                      acCGGCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCg                                <  1:998088/71‑1 (MQ=255)
                                       ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTa                                                            >  1:2208333/1‑42 (MQ=255)
                                       ccGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCgg                               <  1:1327464/71‑1 (MQ=255)
                                       ccGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                               >  1:3422674/1‑71 (MQ=255)
                                        cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                              <  1:3416414/71‑1 (MQ=255)
                                        cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                              <  1:858281/71‑1 (MQ=255)
                                        cGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCggg                              <  1:3959433/71‑1 (MQ=255)
                                         gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                          >  1:5010318/1‑58 (MQ=255)
                                         gTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                          >  1:1857874/1‑58 (MQ=255)
                                            tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                            <  1:3540171/53‑1 (MQ=255)
                                            tGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                            <  1:3858501/53‑1 (MQ=255)
                                             ggggTGTTACCCGTTTTCTTACCATACACAACGTTAg                                                           <  1:1368028/37‑1 (MQ=255)
                                                   ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGt                                                       >  1:5545369/1‑35 (MQ=12)
                                                   ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGa                                                     <  1:18036/37‑1 (MQ=12)
                                                   ttACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGta                                 >  1:2900815/1‑56 (MQ=255)
                                                   ttACCCGTTTTCTTACCATACACAACGTTAGAAGt                                                       <  1:3801679/35‑1 (MQ=255)
                                                   ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa                                            <  1:1265460/46‑1 (MQ=255)
                                                   ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACg                   >  1:4274323/1‑71 (MQ=255)
                                                   ttACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGCAt                            <  1:4533783/62‑1 (MQ=255)
                                                    tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAgg                               <  1:4226166/58‑1 (MQ=255)
                                                    tACCGGTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGt                       <  1:3955736/66‑1 (MQ=255)
                                                     aCCGGTTTTCTGGCCGTACACCACGTTAGAAGTGAta                                                   >  1:4473590/1‑36 (MQ=12)
                                                     aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                                     >  1:3691654/1‑35 (MQ=255)
                                                     aCCCGTTTTCTTACCATACACAACGTTAGAAGTGa                                                     >  1:1242288/1‑35 (MQ=255)
                                                     aCCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                  >  1:115970/1‑70 (MQ=255)
                                                         gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                  >  1:1062633/1‑66 (MQ=255)
                                                         gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACgg                  >  1:3634530/1‑66 (MQ=255)
                                                         gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt             >  1:2879597/1‑71 (MQ=255)
                                                         gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt             >  1:2541054/1‑71 (MQ=255)
                                                         gTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgt             >  1:4710746/1‑71 (MQ=255)
                                                         gTTTTCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACTGGGTAGGGACGGCGTTGCGATAGgt             >  1:2102100/1‑71 (MQ=255)
                                                          ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                          <  1:5667625/41‑1 (MQ=255)
                                                          ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                          <  1:2378405/41‑1 (MQ=255)
                                                          ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACa                                          <  1:38625/41‑1 (MQ=255)
                                                          ttttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCAc                    <  1:4317281/63‑1 (MQ=255)
                                                           tttCTGGCCGTACACCACGTTAGAAGTGATAGTCAGAATCGACt                                      >  1:4245046/1‑44 (MQ=12)
                                                            ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGgtgt           <  1:4831545/70‑1 (MQ=255)
                                                            ttCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGtgtg          >  1:1910693/1‑71 (MQ=255)
                                                            ttCTTACCATACACAACGTTAGAAGAGATGGTCAGAACAGACTGAGTCGGGATAg                          >  1:3934166/1‑55 (MQ=255)
                                                             tCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGa                                        <  1:2128709/40‑1 (MQ=255)
                                                                 aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCa        <  1:3568046/68‑1 (MQ=255)
                                                                 aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt      <  1:4449503/70‑1 (MQ=255)
                                                                 aCCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGt      >  1:1920407/1‑70 (MQ=255)
                                                                  ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATATCGTCACGGTAGGTGTGCAGttt    <  1:82175/71‑1 (MQ=255)
                                                                  ccATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAg       <  1:5049511/68‑1 (MQ=255)
                                                                    aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGt                 <  1:3255532/56‑1 (MQ=255)
                                                                    aTACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCt  >  1:1650878/1‑71 (MQ=255)
                                                                    |                                                                      
GGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAACAGACTGAGTCGGGATAGCGTCACGGTAGGTGTGCAGTTTCT  >  Exported/608347‑608485

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: