Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 609,291 A→T 10.0% T321T (ACT→ACA pflB ← pyruvate formate lyase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported609,2910AT10.0% 264.8 / 15.7 98T321T (ACT→ACApflBpyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (55/33);  new base T (7/3);  total (62/36)
Fisher's exact test for biased strand distribution p-value = 7.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTT  >  Exported/609223‑609358
                                                                    |                                                                   
aGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAg                                                                    >  1:4514107/1‑70 (MQ=255)
    ccATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTGCGGAGTACg                                                                >  1:3484677/1‑70 (MQ=255)
      aTACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAg                                                             >  1:4629968/1‑71 (MQ=255)
       tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                            <  1:3820930/71‑1 (MQ=255)
       tACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAgg                                                            <  1:1712763/71‑1 (MQ=255)
         ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                                                           <  1:821195/70‑1 (MQ=255)
         ccaccGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGa                                                           <  1:3361070/70‑1 (MQ=255)
             cGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                                       <  1:1657633/70‑1 (MQ=255)
               aTAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                                                    >  1:4952182/1‑71 (MQ=255)
                tAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                                      <  1:2988351/68‑1 (MQ=255)
                 aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGc                                                        >  1:1460356/1‑65 (MQ=255)
                 aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                   >  1:5282517/1‑70 (MQ=255)
                 aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAc                                                    >  1:4756260/1‑69 (MQ=255)
                 aGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                                  >  1:399609/1‑71 (MQ=255)
                  gATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                                  <  1:4702818/70‑1 (MQ=255)
                  acttCCGTCGCCCAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCa                                                  >  1:5649674/3‑70 (MQ=255)
                      cGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACg                                              >  1:4007674/1‑70 (MQ=255)
                           gCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGtt                                         >  1:1803539/1‑70 (MQ=255)
                            cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                                       <  1:3545860/71‑1 (MQ=255)
                            cccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTc                                       <  1:5186223/71‑1 (MQ=255)
                             ccAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                      >  1:3512965/1‑71 (MQ=255)
                              cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                                     <  1:2266684/71‑1 (MQ=255)
                              cAGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                                     <  1:3130415/71‑1 (MQ=255)
                              cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                                       <  1:42153/53‑1 (MQ=255)
                              cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                                       <  1:3613448/53‑1 (MQ=255)
                              cAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                                       <  1:4897229/53‑1 (MQ=255)
                               aGATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCAt                                     >  1:3455217/1‑70 (MQ=255)
                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                     >  1:4891136/1‑70 (MQ=255)
                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                     >  1:3794498/1‑70 (MQ=255)
                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                     >  1:3589419/1‑70 (MQ=255)
                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                     >  1:123839/1‑70 (MQ=255)
                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                     >  1:2193704/1‑70 (MQ=255)
                               aGATCGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCGAACCATACGCAGTTTCAt                                     >  1:4897799/1‑70 (MQ=255)
                                gATTGGGTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                   >  1:3407185/1‑71 (MQ=255)
                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTTCGGAGTACGCAGGAAGCGa                                                      >  1:190192/1‑52 (MQ=255)
                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGc                                                               <  1:2008814/43‑1 (MQ=255)
                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                                     >  1:4145257/1‑53 (MQ=255)
                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                                     >  1:5165173/1‑53 (MQ=255)
                                gATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGaa                                                     >  1:2798119/1‑53 (MQ=255)
                                 aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                                                     >  1:440374/1‑36 (MQ=255)
                                 aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGa                                                                     >  1:5447180/1‑36 (MQ=255)
                                 aTCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                                  <  1:3607582/71‑1 (MQ=255)
                                  ttgggTCGCCGGAGAACAGCGAATCAAATTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                   >  1:1024587/3‑69 (MQ=255)
                                  tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCa                                                              >  1:352507/1‑42 (MQ=255)
                                  tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGa                                                      >  1:1816350/1‑50 (MQ=255)
                                  tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                                >  1:3601350/1‑56 (MQ=255)
                                  tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                                >  1:23699/1‑56 (MQ=255)
                                  tCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                 <  1:3888221/71‑1 (MQ=255)
                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                   >  1:3362506/1‑52 (MQ=255)
                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                   >  1:5637492/1‑52 (MQ=255)
                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                   >  1:2447193/1‑52 (MQ=255)
                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAAcc                                                   >  1:5177806/1‑52 (MQ=255)
                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                        >  1:722803/1‑63 (MQ=255)
                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                      >  1:5626028/1‑65 (MQ=255)
                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                 <  1:2468755/70‑1 (MQ=255)
                                   cGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAcc                                 <  1:5527717/70‑1 (MQ=255)
                                   cGGGTCGCCAGAGAACAGTTCATCGTATTACGGAGTACGCAGGAAGCg                                                       >  1:5453611/1‑48 (MQ=255)
                                     ggTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCa                                                  <  1:2814762/51‑1 (MQ=255)
                                        cGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtg                            >  1:4023666/1‑70 (MQ=255)
                                         gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCg                                                       >  1:3567643/1‑42 (MQ=255)
                                         gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGAc                                  >  1:4516847/1‑63 (MQ=255)
                                         gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                          >  1:3271539/1‑71 (MQ=255)
                                         gCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAggtggt                          >  1:5303577/1‑71 (MQ=255)
                                           cAGAGAACAGTTCATCGTATTCCTGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                        <  1:5031927/71‑1 (MQ=255)
                                           cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGt                                          <  1:228596/53‑1 (MQ=255)
                                           cAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                        <  1:278814/71‑1 (MQ=255)
                                            agagAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAg                               >  1:1972603/1‑63 (MQ=255)
                                             gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATa                                                <  1:4348556/45‑1 (MQ=255)
                                             gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                       >  1:878281/1‑69 (MQ=255)
                                             gagaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCACCTGGTCat                       >  1:2931471/1‑69 (MQ=255)
                                               gaACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCa                    <  1:1885619/71‑1 (MQ=255)
                                                 aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGCTTCATGACCAGGTGGTCAACCAtt                  >  1:1719359/1‑71 (MQ=255)
                                                 aCAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGAAGTTTCATGACCAGGTGGTCAACCAtt                  >  1:4089161/1‑71 (MQ=255)
                                                  cAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                  >  1:2999037/1‑70 (MQ=255)
                                                     ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                        <  1:4923374/45‑1 (MQ=255)
                                                     ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                        <  1:4537443/45‑1 (MQ=255)
                                                     ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                        <  1:64744/45‑1 (MQ=255)
                                                     ttCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGttt                                        <  1:2220146/45‑1 (MQ=255)
                                                      tCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAg                                           >  1:2543439/1‑41 (MQ=255)
                                                      tCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAg                                           >  1:4645613/1‑41 (MQ=255)
                                                       cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                      >  1:5633186/1‑45 (MQ=255)
                                                       cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCa                                      >  1:2829690/1‑45 (MQ=255)
                                                       cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCGACCATTTCCtac            >  1:3684893/1‑69 (MQ=255)
                                                       cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc            >  1:3693472/1‑71 (MQ=255)
                                                       cATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTgc            >  1:1981534/1‑71 (MQ=255)
                                                       cATCGTATTCCGGAGTAAGCAGGAAGCGAACCATACGCAGTTTCa                                      >  1:5174622/1‑45 (MQ=255)
                                                          caaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgc          >  1:1560493/4‑70 (MQ=255)
                                                           gTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGctt        >  1:4107251/1‑71 (MQ=255)
                                                           aaaTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGATGAAGTGATCGATCAGTTCCTgtgcc         <  1:939224/68‑2 (MQ=255)
                                                             aTTCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGa                                   >  1:4940931/1‑42 (MQ=255)
                                                              ttCCGGTGTACGCAGGAAGCGTACCATACGGATCTTCATGatg                                 >  1:4958715/1‑41 (MQ=255)
                                                              ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTcta                       >  1:958715/1‑51 (MQ=255)
                                                              ttCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGcttctt     >  1:5283728/1‑71 (MQ=255)
                                                               tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                  <  1:5318468/57‑1 (MQ=255)
                                                               tCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAtt                  <  1:4064226/57‑1 (MQ=255)
                                                                ccGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCaa                       >  1:1607554/1‑51 (MQ=255)
                                                                 cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt  <  1:4578514/71‑1 (MQ=255)
                                                                 cGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGtt  <  1:236672/71‑1 (MQ=255)
                                                                   gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCa                        <  1:3155037/47‑1 (MQ=255)
                                                                   gAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCAttt                 <  1:4530147/54‑1 (MQ=255)
                                                                    |                                                                   
AGGCCCATACCACCGATAGATTCGGTTGCCCAGATCGGGTCGCCAGAGAACAGTTCATCGTATTCCGGAGTACGCAGGAAGCGAACCATACGCAGTTTCATGACCAGGTGGTCAACCATTTCCTGCGCTTCTTGTT  >  Exported/609223‑609358

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: