Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 2,002,416 A→G 11.7% intergenic (+63/+60) 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,4160AG11.7% 180.9 / 12.4 68intergenic (+63/+60)yghU/hybGpredicted S‑transferase/hydrogenase 2 accessory protein
Reads supporting (aligned to +/- strand):  ref base A (36/24);  new base G (3/5);  total (39/29)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.68e-01

AGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGG  >  Exported/2002346‑2002486
                                                                      |                                                                        
aGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                          >  1:628797/1‑71 (MQ=255)
aGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                          <  1:5217679/71‑1 (MQ=255)
          gTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACaa                                                               >  1:990705/1‑71 (MQ=255)
          gTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACaa                                                               <  1:4667310/71‑1 (MQ=255)
              gtgtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGc                                                            <  1:914092/70‑1 (MQ=255)
              gtgtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGc                                                            <  1:3137013/70‑1 (MQ=255)
               tgtTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAAtaga                                                           >  1:4184107/1‑67 (MQ=255)
                 ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTg                                                         <  1:1579360/70‑1 (MQ=255)
                 ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGa                                                        <  1:3091002/71‑1 (MQ=255)
                 ttCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGa                                                        <  1:2235435/71‑1 (MQ=255)
                  tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCATGCCGCATCCGGCAATCGTGCACAATGCCTGAt                                                       >  1:2414298/1‑71 (MQ=255)
                  tcgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGAt                                                       >  1:3347479/1‑71 (MQ=255)
                   cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATg                                                      >  1:3395797/1‑71 (MQ=255)
                   cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATg                                                      >  1:3300464/1‑71 (MQ=255)
                   cgtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATg                                                      >  1:800646/1‑71 (MQ=255)
                    gtcgCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGc                                                     <  1:2251105/71‑1 (MQ=255)
                       gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGAc                                                  <  1:3554549/71‑1 (MQ=255)
                       gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGAc                                                  >  1:3892702/1‑71 (MQ=255)
                       gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGAc                                                  >  1:4786456/1‑71 (MQ=255)
                       gcagcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGAc                                                  >  1:2601440/1‑71 (MQ=255)
                          gcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCt                                              >  1:2718058/1‑71 (MQ=255)
                          gcaAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCt                                              >  1:1706087/1‑71 (MQ=255)
                           caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                             >  1:549561/1‑71 (MQ=255)
                           caAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTg                                             >  1:3460609/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)
                                           gATAAGGCGTTCACGCCGCATCCGGCAATCGTGcac                                                                  <  1:1027287/36‑1 (MQ=255)
                                            aTAAGGCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAgac                             <  1:2140569/71‑3 (MQ=1)
                                            aTAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGcc                            >  1:299414/1‑71 (MQ=255)
                                              aaGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGc                                                >  1:5086911/1‑50 (MQ=255)
                                              aaGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTa                          >  1:502717/1‑71 (MQ=255)
                                               aGGCGTTCACGCCGCATCCGGCAGTCGTGCGCCATGCCTGATGCGACGCTTTGCGCGTCTTATCAGGCCt                           <  1:143691/70‑1 (MQ=11)
                                               aGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCt                                                          <  1:4704908/39‑1 (MQ=255)
                                                ggCGTTCACGCCGCATCCGGCAGTCGTGCGCAGATGCCTGATGCGACGCTAACGCGTCTTATCAGGCCTAc                         <  1:4056090/71‑1 (MQ=1)
                                                ggCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGc                                                     <  1:4087008/43‑1 (MQ=255)
                                                 gCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGa                       <  1:5425052/71‑1 (MQ=255)
                                                  cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                      <  1:2040050/71‑1 (MQ=255)
                                                  cGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGaa                      <  1:4972006/71‑1 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:2818161/1‑71 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:1811793/1‑71 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:2566937/1‑71 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:5073097/1‑71 (MQ=255)
                                                     tCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGc                   >  1:445070/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)
                                                                      |                                                                        
AGGGGTAAGGGTTGGTGTTCGTCGCAGCAAGCCATCCAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAATCGTGCACAATGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGAAAGCTAATCATCACGCGCTGG  >  Exported/2002346‑2002486

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: