Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,595,518 C→T 25.2% A360A (GCG→GCA arsB ← arsenite/antimonite transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,595,5180CT25.2% 223.8 / 78.5 127A360A (GCG→GCAarsBarsenite/antimonite transporter
Reads supporting (aligned to +/- strand):  ref base C (51/44);  new base T (20/12);  total (71/56)
Fisher's exact test for biased strand distribution p-value = 4.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 134 total aligned reads displayed.

GCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATA  >  minE/2595448‑2595587
                                                                      |                                                                     
gcgtagcCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                                                       >  1:2131489/1‑71 (MQ=255)
           cTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGa                                                             >  1:3749195/1‑70 (MQ=255)
              ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                                                          <  1:496167/70‑1 (MQ=255)
              ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATaa                                                          <  1:3592426/70‑1 (MQ=255)
              ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAAc                                                         <  1:686166/71‑1 (MQ=255)
              ccAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAAc                                                         >  1:2744289/1‑71 (MQ=255)
              ccAATCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGacga                                                          <  1:1759982/70‑4 (MQ=255)
               cAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACg                                                        >  1:3597233/1‑71 (MQ=255)
                aaTCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                                                       >  1:1142962/1‑71 (MQ=255)
                aaTCGGGGTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGc                                                       <  1:2056505/71‑1 (MQ=255)
                 aTTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGcc                                                      <  1:448498/71‑1 (MQ=255)
                 aTTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGc                                                       >  1:1826369/1‑70 (MQ=255)
                  ttGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCa                                                     <  1:966097/71‑1 (MQ=255)
                     gggTGATTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgtc                                                  <  1:1039489/71‑3 (MQ=255)
                       gTAATTTTCTGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGc                                                >  1:3629754/1‑71 (MQ=255)
                       gTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                                                       <  1:45937/48‑1 (MQ=255)
                        tAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                               <  1:1716718/71‑1 (MQ=255)
                         gaTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCg                                              >  1:2310222/2‑71 (MQ=255)
                          aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                                                       >  1:3061579/1‑45 (MQ=255)
                          aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATc                                                                       >  1:3104766/1‑45 (MQ=255)
                          aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAAGGCCAGATGCCGt                                             >  1:3657593/1‑71 (MQ=255)
                          aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                                              <  1:3632583/70‑1 (MQ=255)
                          aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                             >  1:2282745/1‑71 (MQ=255)
                          aTTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                             >  1:3045419/1‑71 (MQ=255)
                          aTTTTCGGTCCCAAATCGCACCCAATCACATTGGCATAAACCATc                                                                       >  1:3530414/1‑45 (MQ=255)
                          aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCttt                                                               >  1:3659092/1‑53 (MQ=255)
                          aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCttt                                                               >  1:1575010/1‑53 (MQ=255)
                          aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                                                    <  1:2776004/64‑1 (MQ=255)
                          aTTTTCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                                                    <  1:373726/64‑1 (MQ=255)
                           ttttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGctt                                                                   <  1:1153299/48‑1 (MQ=255)
                             ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                          >  1:1318215/1‑71 (MQ=255)
                             ttCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgct                                          >  1:1834918/1‑71 (MQ=255)
                              tCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctg                                         >  1:1079007/1‑71 (MQ=255)
                              tCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAg                                                    >  1:1243003/1‑60 (MQ=255)
                              tCGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtactc                                         >  1:1335928/1‑70 (MQ=255)
                               cGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGcct                                                                   >  1:3654912/1‑42 (MQ=255)
                               cGGGCCTAAATCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGtactc                                         <  1:1014757/70‑2 (MQ=255)
                                ggTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctgc                                        >  1:1098596/1‑70 (MQ=255)
                                ggTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGcc                                       <  1:2138816/71‑1 (MQ=255)
                                 gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCa                                      <  1:3081292/71‑1 (MQ=255)
                                 gTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCa                                      <  1:3379826/71‑1 (MQ=255)
                                  tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTg                                                              <  1:803192/46‑1 (MQ=255)
                                  tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTg                                                              <  1:2643761/46‑1 (MQ=255)
                                  tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGa                                                             >  1:884215/1‑47 (MQ=255)
                                  tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtg                                            <  1:1221579/64‑1 (MQ=255)
                                  tCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtg                                            <  1:3699881/64‑1 (MQ=255)
                                   cccAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgctg                                         >  1:3733789/1‑66 (MQ=255)
                                   cccAAATCGAAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCAt                                     >  1:556/1‑70 (MQ=255)
                                     cAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGAt                                                  >  1:1126735/1‑55 (MQ=255)
                                     cAAATCGCAGCCAATCACATTGGCATAAAACATCGCTTCTTTGATaa                                                          >  1:1108861/1‑47 (MQ=255)
                                      aaaTCGCAGCCAATCACATTGGCATAAACCATCGctt                                                                   >  1:3549646/1‑37 (MQ=255)
                                      aaaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAAt                                 >  1:2797251/1‑71 (MQ=255)
                                      aaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgt                                                   >  1:609924/1‑52 (MQ=255)
                                      aaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAgt                                                   >  1:2970104/1‑52 (MQ=255)
                                      aaaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAt                                 >  1:3493069/1‑71 (MQ=255)
                                       aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                           >  1:96736/1‑60 (MQ=255)
                                       aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGtgc                                           >  1:3524820/1‑60 (MQ=255)
                                       aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                >  1:2405114/1‑71 (MQ=255)
                                       aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                >  1:2105044/1‑71 (MQ=255)
                                       aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                >  1:899769/1‑71 (MQ=255)
                                       aaTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                >  1:1132875/1‑71 (MQ=255)
                                       aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGta                                            >  1:743546/1‑58 (MQ=255)
                                       aaTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAAtc                                >  1:1814704/1‑70 (MQ=255)
                                        aTCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCa                                      >  1:3412466/1‑64 (MQ=255)
                                        aTCACAGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCaa                                  <  1:3923911/68‑1 (MQ=255)
                                          cGCAGCCAATCACATTGGCATAAACCATCGGTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGa                              >  1:810901/1‑70 (MQ=255)
                                            cAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACaa                           <  1:744751/71‑1 (MQ=255)
                                             aGCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACagc                          >  1:2383520/1‑69 (MQ=255)
                                             aGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAc                          <  1:1588070/71‑1 (MQ=255)
                                              gCCAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAGcg                         <  1:925202/71‑1 (MQ=255)
                                               ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAg                                                    >  1:200458/1‑43 (MQ=255)
                                               ccAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcg                         >  1:2386721/1‑70 (MQ=255)
                                                cAATCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                       >  1:3797257/1‑70 (MQ=255)
                                                cAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgc                        <  1:1894529/70‑1 (MQ=255)
                                                 aaTCACATTGGCATAAATCATTGCCTCTTTGACGACGCCAGTCGCCGTACTCCCGTCAATCGACAgcgca                       <  1:2834251/70‑2 (MQ=255)
                                                 aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATc                                    <  1:3800795/57‑1 (MQ=255)
                                                 aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                      >  1:3837660/1‑71 (MQ=255)
                                                 aaTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACgcgc                      <  1:301922/71‑1 (MQ=255)
                                                  aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGcc                                               <  1:794708/45‑1 (MQ=255)
                                                  aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCg                                              >  1:3140959/1‑46 (MQ=255)
                                                  aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGt                                             >  1:2633605/1‑47 (MQ=255)
                                                  aTCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCa                                      <  1:1290411/54‑1 (MQ=255)
                                                   tCACATTGGCATAAACCATCGCTTCTTTGATAACGCTAGATGCCGTGCTGCCATCAATGGACAACGCGcc                     >  1:3386963/1‑70 (MQ=255)
                                                   tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATg                                >  1:2580474/1‑59 (MQ=255)
                                                   tCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGcc                     >  1:2944604/1‑70 (MQ=255)
                                                     acaTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTTCTGCCATCAATGGACAACGCGCCAAc                  <  1:3891851/71‑1 (MQ=255)
                                                      caTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAAcc                 <  1:1891077/71‑1 (MQ=255)
                                                       aTTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCa                <  1:2081315/71‑1 (MQ=255)
                                                        ttGGCATAAACCATCGCTTCTTTGATATCGCCAGTTGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg               <  1:1018982/71‑1 (MQ=255)
                                                         tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg               >  1:2753474/1‑70 (MQ=255)
                                                         tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg               >  1:283908/1‑70 (MQ=255)
                                                         tGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAg               >  1:331902/1‑70 (MQ=255)
                                                             aTAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCaa                                  <  1:3172012/47‑1 (MQ=255)
                                                             aTAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgc                        <  1:778934/57‑1 (MQ=255)
                                                              tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgc                        >  1:2829819/1‑56 (MQ=255)
                                                              tAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCg          <  1:1832053/70‑1 (MQ=255)
                                                               aaaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc        >  1:3243634/1‑71 (MQ=255)
                                                                aaCCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTc        <  1:62314/70‑1 (MQ=255)
                                                                     tCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCata  >  1:1873291/1‑71 (MQ=255)
                                                                      cGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAAcgc                        <  1:1589313/48‑1 (MQ=255)
                                                                      |                                                                     
GCGTAGCCAGGCTACCAATTGGGGTAATTTTCGGTCCCAAATCGCAGCCAATCACATTGGCATAAACCATCGCTTCTTTGATAACGCCAGATGCCGTGCTGCCATCAATGGACAACGCGCCAACCAGTACCGTCGGCATA  >  minE/2595448‑2595587

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: