Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,595,782 A→T 13.6% G272G (GGT→GGA arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,595,7820AT13.6% 219.3 / 22.8 81G272G (GGT→GGAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base A (40/30);  new base T (6/5);  total (46/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTG  >  Exported/2595717‑2595843
                                                                 |                                                             
tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCTAGGGGGCACCGCGCAGGACTTTACCCg                                                           >  1:1994388/1‑70 (MQ=255)
tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTAccc                                                            >  1:3169769/1‑69 (MQ=255)
tAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCg                                                           >  1:302012/1‑70 (MQ=255)
  aaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTa                                                         <  1:489454/70‑1 (MQ=255)
  aaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAt                                                        >  1:3996750/1‑71 (MQ=255)
  aaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAt                                                        >  1:1811029/1‑71 (MQ=255)
  aaCCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTAt                                                        >  1:824228/1‑71 (MQ=255)
         aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                  >  1:471589/1‑70 (MQ=255)
         aGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGc                                                 <  1:4850556/71‑1 (MQ=255)
          gATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                <  1:5047097/71‑1 (MQ=255)
            tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATg                                              <  1:4321383/71‑1 (MQ=255)
            tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATg                                              <  1:3809251/71‑1 (MQ=255)
            tACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATg                                              <  1:3476054/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)
                tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTccg                                                            >  1:4369955/1‑52 (MQ=255)
                tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcct                                          >  1:3760988/1‑70 (MQ=255)
                tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcct                                          >  1:2363060/1‑70 (MQ=255)
                tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcct                                          >  1:4590037/1‑70 (MQ=255)
                tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcct                                          >  1:2643713/1‑70 (MQ=255)
                tGCCCAGTGAGAAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGAcc                                           >  1:2919275/1‑70 (MQ=255)
                      gtgagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCttttc                                    <  1:3026550/69‑2 (MQ=255)
                      gtgagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCttttc                                    <  1:2800717/69‑2 (MQ=255)
                      gtgagaAAATAACGATCTGCCATGGCGCACCGCGCAGGACTTTTCCGGTGTTGATCGAATGACCTCttttc                                    <  1:4383066/69‑2 (MQ=255)
                        gagaAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                  <  1:4328686/55‑1 (MQ=255)
                        gagaAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                  <  1:1917654/55‑1 (MQ=255)
                        gagaAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGc                                  <  1:1119414/71‑1 (MQ=255)
                         agaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                                      <  1:786577/50‑1 (MQ=255)
                         agaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTa                                                      <  1:2224188/50‑1 (MQ=255)
                          gaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                >  1:1095870/1‑55 (MQ=255)
                          gaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                >  1:3299263/1‑55 (MQ=255)
                          gaagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCa                                                >  1:2173623/1‑55 (MQ=255)
                           aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                        <  1:3610010/62‑1 (MQ=255)
                           aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                        <  1:3255475/62‑1 (MQ=255)
                           aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                        <  1:414333/62‑1 (MQ=255)
                           aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                        <  1:5078090/62‑1 (MQ=255)
                           aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACcgcg                                        <  1:2933304/62‑1 (MQ=255)
                           aagaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCCTATTAATCGCATGACcgcg                                        <  1:728229/62‑1 (MQ=255)
                            agaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTaa                                                     >  1:2548385/1‑48 (MQ=255)
                             gaTGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGtt                                      >  1:4425836/1‑62 (MQ=255)
                               tGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                           <  1:3021015/71‑1 (MQ=255)
                                gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                           >  1:4341645/1‑70 (MQ=255)
                                gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                           >  1:1406680/1‑70 (MQ=255)
                                gACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACa                           >  1:3648213/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)
                                         ccAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCg                                                  <  1:2629714/38‑1 (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:2505166/1‑56 (MQ=255)
                                                 cACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaat                        >  1:1514605/1‑56 (MQ=255)
                                                 cACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaat                        >  1:1948552/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:37659/1‑58 (MQ=255)
                                                     gcgcAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATc                  >  1:351181/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:862131/1‑71 (MQ=255)
                                                        cAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTg  <  1:4945267/71‑1 (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:1987051/48‑1 (MQ=255)
                                                            aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                     <  1:1945603/48‑1 (MQ=255)
                                                            aCTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAaataat                     <  1:1092064/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)
                                                                 |                                                             
TAAACCACCAGATACATGCCGAGCGAGAAGATGACAATCTGCCAGGGGGCACCGCGCAGGACTTTACCCGTATTAATCGCATGACCGCGTTTAGCGACGACAAATAATATCAGCGCGCCCACAGCTG  >  Exported/2595717‑2595843

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: