Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 609,165 G→A 10.5% T363T (ACC→ACT pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported609,1650GA10.5% 152.2 / 7.2 57T363T (ACC→ACTpflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (25/26);  new base A (4/2);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCA  >  Exported/609095‑609234
                                                                      |                                                                     
aTTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATg                                                                       >  1:2855274/1‑71 (MQ=255)
aTTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATg                                                                       >  1:237818/1‑71 (MQ=255)
    cTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGgtg                                                                    <  1:3231649/70‑1 (MQ=255)
    cTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGgtg                                                                    <  1:4803976/70‑1 (MQ=255)
          gTTCAGCGGCAGTTTTTCAGACCACAGAATGGGCATGTTCGGTTCCGGAGACGGACCCATGGGGTACAggg                                                             <  1:630309/71‑1 (MQ=255)
           ttCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGgt                                                            <  1:1283183/71‑1 (MQ=255)
             cAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGgtg                                                           <  1:3811779/70‑1 (MQ=255)
                   cAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAgg                                                     >  1:2101799/1‑70 (MQ=255)
                      tttttCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACa                                                                <  1:5682322/56‑1 (MQ=255)
                       ttttCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTAc                                                                 >  1:5387739/1‑54 (MQ=255)
                        tttCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACa                                                                >  1:3332314/1‑54 (MQ=255)
                         ttCCGACCAAAGAATGGTCAGGTTAGGTTCCGGTGCCGGCCCCATAGTGTGCAg                                                               <  1:1147016/54‑1 (MQ=255)
                         ttCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACgg                                               <  1:5219778/70‑1 (MQ=255)
                          tCAGACCACAGAATGGTCATGTTCGGCTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGa                                              <  1:3006246/70‑1 (MQ=255)
                               ccACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGt                                        >  1:5574060/1‑71 (MQ=255)
                                 acaGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGttt                                      <  1:5498938/71‑1 (MQ=255)
                                  caGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGtttt                                     <  1:3007695/71‑1 (MQ=255)
                                    gAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAg                                                               <  1:5652737/43‑1 (MQ=255)
                                    gAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGttttt                                    >  1:324917/1‑70 (MQ=255)
                                     aaTGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGtttt                                     <  1:4267762/68‑1 (MQ=255)
                                     aaTGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGtttt                                     <  1:333433/68‑1 (MQ=255)
                                       tGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACa                                                                >  1:1275916/1‑39 (MQ=255)
                                          tCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGgtg                                                           >  1:1903412/1‑41 (MQ=255)
                                          tCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACggg                                              >  1:4442518/1‑53 (MQ=255)
                                          tCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACgg                                               >  1:3008823/1‑53 (MQ=255)
                                            aTGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTa                                >  1:3061706/1‑66 (MQ=255)
                                            aTGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAg                           >  1:1055736/1‑71 (MQ=255)
                                             tGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTc                                                        <  1:4275667/41‑1 (MQ=255)
                                             tGTTCGGTTCCGGAGACGGACCCATGATGTACAGGGTGTTc                                                        <  1:4775169/41‑1 (MQ=255)
                                              gTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAggg                         >  1:1284140/1‑71 (MQ=255)
                                               ttCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAg                                            <  1:3889965/51‑1 (MQ=255)
                                               ttCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGtt                                       >  1:3038035/1‑56 (MQ=255)
                                               ttCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTaa                               >  1:260483/1‑64 (MQ=255)
                                               ttCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAggg                         >  1:675934/1‑70 (MQ=255)
                                                tCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTa                       >  1:251769/1‑71 (MQ=255)
                                                tCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTa                       >  1:359129/1‑71 (MQ=255)
                                                 cGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTAc                      >  1:5628799/1‑71 (MQ=255)
                                                  ggTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGa                                              <  1:1371647/46‑1 (MQ=255)
                                                  ggTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGa                                              <  1:663643/46‑1 (MQ=255)
                                                  ggTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGa                                              <  1:3070475/46‑1 (MQ=255)
                                                      ccGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCg                 >  1:3917869/1‑71 (MQ=255)
                                                       cGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCg                 <  1:4329675/70‑1 (MQ=255)
                                                        ggAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTc               <  1:5267283/71‑1 (MQ=255)
                                                        ggAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTc               <  1:4430501/71‑1 (MQ=255)
                                                           gaCGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAg            <  1:1688737/71‑1 (MQ=255)
                                                              ggccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc         >  1:88760/4‑71 (MQ=255)
                                                              ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAgg                          >  1:3280802/1‑54 (MQ=255)
                                                              ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc         <  1:4683451/71‑1 (MQ=255)
                                                              ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc         >  1:4277363/1‑71 (MQ=255)
                                                              ggACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGcc         <  1:831753/71‑1 (MQ=255)
                                                               gccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc         >  1:1969608/3‑70 (MQ=255)
                                                               gccccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGccc         >  1:5585024/3‑70 (MQ=255)
                                                                 cccATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCAt      >  1:2876034/1‑71 (MQ=255)
                                                                  ccATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGa             >  1:1727957/1‑63 (MQ=255)
                                                                  ccATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGa             >  1:3538263/1‑63 (MQ=255)
                                                                   cATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac    <  1:1498445/71‑1 (MQ=255)
                                                                   cATAGTGTGCAGGGTGTGCAAATAGCGGAAGGAGTTTTTGGTCACCAGCGTACGACCGTCCAGCCCCAtcc    <  1:28585/71‑3 (MQ=255)
                                                                    aTGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATac    >  1:3725049/1‑70 (MQ=255)
                                                                     tGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATacca  >  1:1355781/1‑71 (MQ=255)
                                                                      |                                                                     
ATTTCTTGAAGTTCAGCGGCAGTTTTTCAGACCACAGAATGGTCATGTTCGGTTCCGGAGACGGACCCATGGTGTACAGGGTGTTCAGGAAACGGAAGCTGTTTTTGGTAACCAGGGTACGACCGTCGAGGCCCATACCA  >  Exported/609095‑609234

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: