Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,595,796 C→A 12.0% A268S (GCG→TCG)  arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,595,7960CA12.0% 294.6 / 15.8 91A268S (GCG→TCG) arsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (44/36);  new base A (5/6);  total (49/42)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGCC  >  Exported/2595726‑2595865
                                                                      |                                                                     
aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGc                                                                       <  1:4850556/71‑1 (MQ=255)
 gATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                                      <  1:5047097/71‑1 (MQ=255)
   tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATg                                                                    <  1:4321383/71‑1 (MQ=255)
   tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATg                                                                    <  1:3476054/71‑1 (MQ=255)
   tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATg                                                                    <  1:3809251/71‑1 (MQ=255)
   tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGAAt                                                                     <  1:196274/70‑1 (MQ=255)
    aCATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                                                                   >  1:1059677/1‑71 (MQ=255)
    aCATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                                                                   >  1:2745877/1‑71 (MQ=255)
       tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcct                                                                >  1:2643713/1‑70 (MQ=255)
       tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcct                                                                >  1:3760988/1‑70 (MQ=255)
       tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcct                                                                >  1:4590037/1‑70 (MQ=255)
       tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcct                                                                >  1:2363060/1‑70 (MQ=255)
       tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcc                                                                 >  1:2919275/1‑70 (MQ=255)
             gtgagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCttttc                                                          <  1:2800717/69‑2 (MQ=255)
             gtgagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCttttc                                                          <  1:4383066/69‑2 (MQ=255)
             gtgagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCttttc                                                          <  1:3026550/69‑2 (MQ=255)
               gagaAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGc                                                        <  1:1119414/71‑1 (MQ=255)
                 gaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                                      >  1:3299263/1‑55 (MQ=255)
                 gaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                                      >  1:1095870/1‑55 (MQ=255)
                 gaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                                      >  1:2173623/1‑55 (MQ=255)
                  aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                              <  1:2933304/62‑1 (MQ=255)
                  aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                              <  1:5078090/62‑1 (MQ=255)
                  aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                              <  1:3610010/62‑1 (MQ=255)
                  aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                              <  1:3255475/62‑1 (MQ=255)
                  aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                                              <  1:414333/62‑1 (MQ=255)
                  aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCCTATTAATCGCATGACcgcg                                                              <  1:728229/62‑1 (MQ=255)
                    gaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGtt                                                            >  1:4425836/1‑62 (MQ=255)
                      tGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                                 <  1:3021015/71‑1 (MQ=255)
                       gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                                 >  1:3648213/1‑70 (MQ=255)
                       gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                                 >  1:1406680/1‑70 (MQ=255)
                       gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                                                 >  1:4341645/1‑70 (MQ=255)
                        aCAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGa                                                                   >  1:1650575/1‑51 (MQ=255)
                        aCAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGAc                                                                  >  1:5147509/1‑52 (MQ=255)
                        aCAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACaa                                                <  1:1789208/70‑1 (MQ=255)
                              tGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGcgacga                                                   >  1:367128/1‑61 (MQ=255)
                                ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATc                                        <  1:2419454/70‑1 (MQ=255)
                                      cgCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCttttc                                                          <  1:4579620/45‑2 (MQ=255)
                                       tcACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCttttc                                                          <  1:2404317/44‑2 (MQ=255)
                                        cACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaat                                              >  1:1948552/1‑56 (MQ=255)
                                        cACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaat                                              >  1:2505166/1‑56 (MQ=255)
                                        cACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaat                                              >  1:1514605/1‑56 (MQ=255)
                                            gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATc                                        >  1:365335/1‑58 (MQ=255)
                                            gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATc                                        >  1:3270335/1‑58 (MQ=255)
                                            gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATc                                        >  1:816324/1‑58 (MQ=255)
                                            gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATc                                        >  1:351181/1‑58 (MQ=255)
                                            gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATc                                        >  1:37659/1‑58 (MQ=255)
                                            gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATc                                        >  1:3594951/1‑58 (MQ=255)
                                            gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTGTAGCGACGACAAATAATATc                                        >  1:3430696/1‑58 (MQ=255)
                                               cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCt                         >  1:3385120/1‑70 (MQ=255)
                                               cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg                        <  1:4945267/71‑1 (MQ=255)
                                               cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg                        >  1:862131/1‑71 (MQ=255)
                                                 ggACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaat                                              <  1:1812008/47‑1 (MQ=255)
                                                  gACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                           <  1:701102/49‑1 (MQ=255)
                                                  gACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                           <  1:963849/49‑1 (MQ=255)
                                                  gACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                           <  1:3044714/49‑1 (MQ=255)
                                                   aatttACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAg                                      >  1:3433485/3‑53 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                           <  1:1092064/48‑1 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                           <  1:1945603/48‑1 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                           <  1:1987051/48‑1 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                                            >  1:4658755/1‑47 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                                            >  1:2180914/1‑47 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                                            >  1:2027668/1‑47 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataa                                            >  1:1286755/1‑47 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAg                                      >  1:4102148/1‑53 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAg                                      >  1:3809427/1‑53 (MQ=255)
                                                   aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAg                                      >  1:906988/1‑53 (MQ=255)
                                                         cccGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCt              <  1:3293853/71‑1 (MQ=255)
                                                         cccGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCt              <  1:5062088/71‑1 (MQ=255)
                                                         cccGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCt              <  1:2250904/71‑1 (MQ=255)
                                                          ccGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCa                                       <  1:497072/45‑1 (MQ=255)
                                                           cGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                                           <  1:72146/40‑1 (MQ=255)
                                                            gTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGAAATGGCGCTCAc           >  1:397069/1‑71 (MQ=255)
                                                              aTTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCAcc          <  1:2219978/70‑1 (MQ=255)
                                                              aTTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCAcc          <  1:3593877/70‑1 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                           <  1:2406281/36‑1 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                           >  1:2499246/1‑36 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                           >  1:5032047/1‑36 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                           >  1:447674/1‑36 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                           >  1:5070175/1‑36 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                           >  1:3784264/1‑36 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                           >  1:1698826/1‑36 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                           >  1:4207923/1‑36 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAaataat                                           >  1:3468552/1‑36 (MQ=255)
                                                               ttAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCgg        <  1:3544329/71‑1 (MQ=255)
                                                                tAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCgg        <  1:4461401/70‑1 (MQ=255)
                                                                tAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCgg        >  1:4426727/1‑70 (MQ=255)
                                                                  aTCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGaa      >  1:4579846/1‑70 (MQ=255)
                                                                   tCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATg    <  1:3668930/71‑1 (MQ=255)
                                                                   tCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATg    <  1:2572019/71‑1 (MQ=255)
                                                                     gCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGcc  <  1:3168242/71‑1 (MQ=255)
                                                                     gCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGcc  <  1:2468804/71‑1 (MQ=255)
                                                                     gCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGcc  >  1:4622874/1‑71 (MQ=255)
                                                                      |                                                                     
AGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTGCAATGGCGCTCACCGGAATGCC  >  Exported/2595726‑2595865

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: