Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,002,443 G→A 7.4% intergenic (+90/+33) yghU → / ← hybG predicted S‑transferase/hydrogenase 2 accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported2,002,4430GA 7.4% 292.7 / 8.0 94intergenic (+90/+33)yghU/hybGpredicted S‑transferase/hydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base G (46/40);  new base A (3/4);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

CAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGC  >  Exported/2002373‑2002509
                                                                        |                                                                  
caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                                    >  1:3460609/1‑71 (MQ=255)
caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                                                    >  1:549561/1‑71 (MQ=255)
  aGCCATCCAGGCCGGATAAGGCGTTCCCGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGcc                                                                  >  1:300041/1‑71 (MQ=255)
  aGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGc                                                                   >  1:5213317/1‑70 (MQ=255)
     cATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCcgcg                                                               >  1:3412365/1‑71 (MQ=255)
         taGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCt                                                            >  1:1747449/2‑71 (MQ=11)
         cAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCt                                                            >  1:547717/1‑70 (MQ=255)
          aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTa                                                          <  1:2424158/71‑1 (MQ=255)
           ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCtt                                                           <  1:2369869/70‑1 (MQ=1)
           ggCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTa                                                          >  1:5575257/1‑71 (MQ=1)
           ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                                         <  1:398235/71‑1 (MQ=255)
           ggCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTAt                                                         <  1:1976584/71‑1 (MQ=255)
            gCCGGATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTAt                                                         <  1:366548/71‑1 (MQ=1)
             ccGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCa                                                       <  1:2557997/71‑1 (MQ=255)
                gATAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAgg                                                     >  1:103501/1‑71 (MQ=1)
                 aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAgac                                                    <  1:2140569/71‑3 (MQ=1)
                 aTAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGcc                                                   >  1:299414/1‑71 (MQ=255)
                   aaGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTa                                                 >  1:502717/1‑71 (MQ=255)
                    aGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCt                                                  <  1:143691/70‑1 (MQ=11)
                     ggCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTAc                                                <  1:4056090/71‑1 (MQ=1)
                      gCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGa                                              <  1:5425052/71‑1 (MQ=255)
                       cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                                             <  1:4972006/71‑1 (MQ=255)
                       cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                                             <  1:2040050/71‑1 (MQ=255)
                          tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                          >  1:1811793/1‑71 (MQ=255)
                          tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                          >  1:2818161/1‑71 (MQ=255)
                          tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                          >  1:2566937/1‑71 (MQ=255)
                          tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                          >  1:445070/1‑71 (MQ=255)
                          tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                          >  1:5073097/1‑71 (MQ=255)
                           cACGCCGCATCCGGCAATCGTTCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                          <  1:2863666/70‑1 (MQ=255)
                           cACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGa                                              >  1:448828/1‑66 (MQ=255)
                           cACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                                          >  1:1488499/1‑70 (MQ=255)
                            aCGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCt                                         <  1:270106/70‑1 (MQ=255)
                            aCGCCGCATCCGGCAATCGAGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTa                                        <  1:3990030/71‑1 (MQ=255)
                               ccgcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                       >  1:3214259/1‑69 (MQ=255)
                               ccgcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                                     <  1:5619053/71‑1 (MQ=255)
                                 gcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                   <  1:1283253/71‑1 (MQ=255)
                                 gcATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                   <  1:1862813/71‑1 (MQ=255)
                                     ccGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcg                               >  1:853808/1‑71 (MQ=255)
                                      cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                       >  1:1066588/1‑62 (MQ=255)
                                      cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc                              >  1:3899269/1‑71 (MQ=255)
                                      cGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc                              >  1:2264226/1‑71 (MQ=255)
                                        gCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc                            >  1:313665/1‑71 (MQ=255)
                                         cAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt                           >  1:3221371/1‑71 (MQ=255)
                                         cAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCt                           >  1:4051214/1‑71 (MQ=255)
                                          aaTCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                                             <  1:4997144/52‑1 (MQ=255)
                                          aaTCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTg                          >  1:4606671/1‑71 (MQ=255)
                                           aTCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                         <  1:1021144/71‑1 (MQ=255)
                                            tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                                     >  1:5196056/1‑58 (MQ=255)
                                            tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                         >  1:4487556/1‑70 (MQ=255)
                                            tCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt                        <  1:5235391/71‑1 (MQ=255)
                                             cGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                   >  1:4082724/1‑59 (MQ=255)
                                             cGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgcg                             >  1:3916511/1‑65 (MQ=255)
                                             cGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGt                        <  1:5005017/70‑1 (MQ=255)
                                              gTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatca                                    <  1:307170/57‑1 (MQ=255)
                                              gTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                      >  1:2209611/1‑71 (MQ=255)
                                               tGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAAt                     <  1:3547729/71‑1 (MQ=255)
                                                 cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg                    <  1:2037908/70‑1 (MQ=255)
                                                 cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg                    <  1:3538600/70‑1 (MQ=255)
                                                 cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATg                    <  1:3390203/70‑1 (MQ=255)
                                                 cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                   <  1:2162253/71‑1 (MQ=255)
                                                 cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                   <  1:887715/71‑1 (MQ=255)
                                                 cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                   >  1:4828888/1‑71 (MQ=255)
                                                 cacaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                   <  1:3471460/71‑1 (MQ=255)
                                                  acaATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGt                   <  1:5081828/70‑1 (MQ=255)
                                                   caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACgcgc                            <  1:1787944/60‑1 (MQ=255)
                                                   caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTg                          >  1:4726473/1‑62 (MQ=255)
                                                   caATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg                 <  1:615478/71‑1 (MQ=255)
                                                    aaTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCg                 >  1:3729234/1‑70 (MQ=255)
                                                      aTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                         <  1:1094859/61‑1 (MQ=255)
                                                      aTGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTgg                         <  1:3056775/61‑1 (MQ=255)
                                                        gCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                   <  1:1928746/49‑1 (MQ=255)
                                                        gCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                      >  1:3842880/1‑62 (MQ=255)
                                                        gCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAAt             >  1:3180678/1‑71 (MQ=255)
                                                         ccTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAAt             <  1:4275079/70‑1 (MQ=255)
                                                         ccTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATc            >  1:499980/1‑71 (MQ=255)
                                                         ccTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATc            >  1:506836/1‑71 (MQ=255)
                                                             aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAg                                           >  1:5421369/1‑36 (MQ=255)
                                                             aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTaa                                       <  1:16249/40‑1 (MQ=255)
                                                             aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCa          <  1:522617/69‑1 (MQ=255)
                                                             aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCa          >  1:4843897/1‑69 (MQ=255)
                                                             aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCa          <  1:5367826/69‑1 (MQ=255)
                                                             aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAtt        >  1:179636/1‑71 (MQ=255)
                                                             aTGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAt         <  1:3104497/70‑1 (MQ=255)
                                                              tGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatcat                                   <  1:2751855/43‑1 (MQ=255)
                                                              tGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCAttt       >  1:1786443/1‑71 (MQ=255)
                                                               gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCAcgc                              >  1:5508818/1‑47 (MQ=255)
                                                               gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                      <  1:5073057/55‑1 (MQ=255)
                                                               gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                      <  1:369893/55‑1 (MQ=255)
                                                               gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTg      >  1:2480610/1‑71 (MQ=255)
                                                               gCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTg      >  1:1637329/1‑71 (MQ=255)
                                                                cGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAAtcatca                                 >  1:1548345/1‑43 (MQ=255)
                                                                   cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                                     >  1:1670158/1‑36 (MQ=255)
                                                                   cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAatc                                     >  1:4582596/1‑36 (MQ=255)
                                                                   cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTaa                      >  1:4973019/1‑51 (MQ=255)
                                                                   cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGgcgc  <  1:4883610/71‑1 (MQ=255)
                                                                   cGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGgcg   >  1:3603018/1‑70 (MQ=255)
                                                                        |                                                                  
CAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGGTAATGTCGTAATCCATTTGGCGC  >  Exported/2002373‑2002509

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: