Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,250,324 T→A 5.8% D123V (GAC→GTC)  dhaH ← fused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,250,3240TA 5.8% 293.7 / 8.2 104D123V (GAC→GTC) dhaHfused predicted dihydroxyacetone‑specific PTS enzyme HPr component and EI component
Reads supporting (aligned to +/- strand):  ref base T (50/48);  new base A (3/3);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

Only 100 of 106 total aligned reads displayed.

AGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAGTGTACCTTCGAC  >  AP009048/1250257‑1250391
                                                                   |                                                                   
agagATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCa                                                                   >  1:3847043/1‑70 (MQ=255)
agagATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCa                                                                   <  1:1501998/70‑1 (MQ=255)
 gagaTTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCa                                                                   >  1:2378054/1‑69 (MQ=255)
  agaTTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAg                                                                <  1:436550/71‑1 (MQ=255)
  agaTTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAg                                                                <  1:6837382/71‑1 (MQ=255)
   gaTTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCa                                                                   >  1:3890039/1‑67 (MQ=255)
   gaTTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCa                                                                   >  1:2898535/1‑67 (MQ=255)
      ttCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATaa                                                            >  1:7766216/1‑71 (MQ=255)
       tCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATaa                                                            <  1:179899/70‑1 (MQ=255)
       tCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATaa                                                            <  1:1660837/70‑1 (MQ=255)
       tCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATaa                                                            <  1:2764037/70‑1 (MQ=255)
         aGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACtt                                                         >  1:2275459/1‑71 (MQ=255)
         aGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACtt                                                         >  1:3358083/1‑71 (MQ=255)
         aGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACt                                                          <  1:7008566/70‑1 (MQ=255)
         aGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCAAGCGCATGCATGGCGTCAAAGATAACt                                                          <  1:7015524/70‑1 (MQ=255)
         aGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTACCAGCGCATGCATGGCGTCaaa                                                                 <  1:7751085/63‑1 (MQ=255)
          gtggCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACttt                                                        <  1:4290107/71‑1 (MQ=255)
           tgtCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGa                                                               >  1:1062619/1‑63 (MQ=255)
           tgtCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATa                                                             >  1:5311285/1‑65 (MQ=255)
           tgtCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATa                                                             >  1:6036363/1‑65 (MQ=255)
           tgtCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTg                                                       <  1:2952593/71‑1 (MQ=255)
           tgtCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTg                                                       <  1:4635861/71‑1 (MQ=255)
            gtCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTg                                                       >  1:447149/1‑70 (MQ=255)
            gtCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGt                                                      >  1:2292102/1‑71 (MQ=255)
            gtCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGt                                                      >  1:210586/1‑71 (MQ=255)
             tCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGGGCATGCATGGCGTCAAAGATAACTTTGTc                                                     >  1:1084740/1‑71 (MQ=255)
             tCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATATCTTTGTc                                                     >  1:5482042/1‑71 (MQ=255)
             tCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGt                                                      <  1:7320531/70‑1 (MQ=255)
             tCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGt                                                      <  1:6583428/70‑1 (MQ=255)
             tCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGt                                                      <  1:7604461/70‑1 (MQ=255)
              cGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTc                                                     >  1:1706385/1‑70 (MQ=255)
                gaggaCGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGat                                                  <  1:2749459/71‑1 (MQ=255)
                  ggaCGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGatat                                                <  1:4889147/71‑1 (MQ=255)
                  ggaCGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGatat                                                <  1:4503497/71‑1 (MQ=255)
                     cGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCg                                             >  1:1394334/1‑71 (MQ=255)
                     cGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCg                                             >  1:2608125/1‑71 (MQ=255)
                       gTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATaa                                                            <  1:2313013/54‑1 (MQ=255)
                       gTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGcc                                           >  1:1909805/1‑71 (MQ=255)
                        tAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGccc                                          >  1:1754914/1‑71 (MQ=255)
                        tAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGccc                                          >  1:3968818/1‑71 (MQ=255)
                        tAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGccc                                          >  1:5520827/1‑71 (MQ=255)
                        tAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGccc                                          >  1:3890056/1‑71 (MQ=255)
                         aaaCCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATc                                               >  1:5868368/1‑65 (MQ=255)
                           aCCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATaa                                                            >  1:2868928/1‑50 (MQ=255)
                            cccAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGa                                      <  1:1394056/71‑1 (MQ=255)
                             ccAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGa                                      <  1:6178527/70‑1 (MQ=255)
                             ccAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAg                                     <  1:791887/71‑1 (MQ=255)
                             ccAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAg                                     <  1:1263219/71‑1 (MQ=255)
                               aGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAgg                                    >  1:5693256/1‑70 (MQ=255)
                               aGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGc                                   <  1:2220712/71‑1 (MQ=255)
                               aGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGc                                   <  1:4985820/71‑1 (MQ=255)
                                 ttgttCACGTTTGGCTTCCAGCGCATGCATGGCGTc                                                                    <  1:2388295/36‑1 (MQ=255)
                                 ttgttCACGTTTGGCTTCCAGCGCATGCATGGCGTc                                                                    <  1:6530815/36‑1 (MQ=255)
                                 ttgttCACGTTTGGCTTCCAGCGCATGCATGGCGTc                                                                    <  1:1535484/36‑1 (MQ=255)
                                 ttgttCACGTTTGGCTTCCAGCGCATGCATGGCGTCa                                                                   <  1:6873049/37‑1 (MQ=255)
                                 ttgttCACGTTTGGCTTCCAGCGCATGCATGGCGTCa                                                                   <  1:2407970/37‑1 (MQ=255)
                                   gttCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACt                                                          <  1:5590474/44‑1 (MQ=255)
                                   gttCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACt                                                          <  1:5453100/44‑1 (MQ=255)
                                   gttCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACt                                                          <  1:3652069/44‑1 (MQ=255)
                                   gttCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACt                                                          <  1:2409381/44‑1 (MQ=255)
                                    ttCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCgcgc                              <  1:3113135/71‑1 (MQ=255)
                                     tCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCt                             >  1:6840308/1‑71 (MQ=255)
                                     tCACGTTTGGCTTCCAGCGCATGCATGGCGACAAAGTTATCt                                                          <  1:27713/42‑1 (MQ=37)
                                       aCGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGat                                                  <  1:4275096/48‑1 (MQ=255)
                                       aCGTTTGGCTTCCAGCGCATGCATGGCGACAAAGTTATCttt                                                        >  1:1883557/1‑42 (MQ=25)
                                       aCGTTTGGCTTCCAGCGCATGCATGGCGACAAAGTTATCttt                                                        >  1:1488091/1‑42 (MQ=25)
                                       aCGTTTGGCTTCCAGCGCATGCATGGCGACAAAGTTATCttt                                                        >  1:2289224/1‑42 (MQ=25)
                                        cGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGAc                          >  1:7824322/1‑71 (MQ=255)
                                          tttGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTg                                                       <  1:5930286/40‑1 (MQ=255)
                                          tttGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCg                        <  1:2977470/71‑1 (MQ=255)
                                            tGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGt                                                      >  1:193006/1‑39 (MQ=255)
                                            tGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTc                                                     >  1:1973741/1‑40 (MQ=255)
                                            tGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGAcc                         >  1:5112412/1‑68 (MQ=255)
                                               cTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGccccc                                        <  1:1280106/50‑1 (MQ=255)
                                                ttCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGa                                                   >  1:33852/1‑38 (MQ=255)
                                                ttCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGa                                                   >  1:228234/1‑38 (MQ=255)
                                                 tCCAGCGCATGCATGGCGTCAAAGATAACTTTGTc                                                     >  1:4183617/1‑35 (MQ=255)
                                                 tCCAGCGCATGCATGGCGTCAAAGATAACTTTGTc                                                     >  1:2180010/1‑35 (MQ=255)
                                                 tCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGa                                      >  1:547179/1‑50 (MQ=255)
                                                 tCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGAc                          >  1:1608813/1‑62 (MQ=255)
                                                  ccAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTtgctg                  <  1:1436459/69‑1 (MQ=255)
                                                    aGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAg              >  1:4524779/1‑71 (MQ=255)
                                                       gcatgcatGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCg                        <  1:1627581/58‑1 (MQ=255)
                                                         atgcatgGCGTCAAAGATAACTTTGTCGATATCCGCCCCCg                                       >  1:6161633/1‑41 (MQ=255)
                                                          tgcatgGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAGTGTAcc        <  1:3966548/71‑1 (MQ=255)
                                                           gcatgGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAGTGTAcc        <  1:4830609/70‑1 (MQ=255)
                                                           gcatgGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAGTGTACCt       >  1:1910279/1‑71 (MQ=255)
                                                            catgGCGTCAAAGATAACTTTGTCGATATCCGccccc                                        <  1:6364699/37‑1 (MQ=255)
                                                            catgGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAGTGTACCt       <  1:6617587/70‑1 (MQ=255)
                                                             atgGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAg              >  1:5931822/1‑62 (MQ=255)
                                                             atgGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAGTGTACCTTc     >  1:1075832/1‑71 (MQ=255)
                                                               ggCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAgg                                    >  1:2899143/1‑38 (MQ=255)
                                                               ggCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGt                       >  1:6260086/1‑51 (MQ=255)
                                                               ggCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGt                       >  1:1569647/1‑51 (MQ=255)
                                                               ggCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAgtgt           >  1:4295170/1‑63 (MQ=255)
                                                               ggCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAGTGTACCTTCGa   >  1:4990305/1‑71 (MQ=255)
                                                                gCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAGTGTACCTTCGAc  >  1:1911695/1‑71 (MQ=255)
                                                                  gTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAgt             <  1:1891971/58‑1 (MQ=255)
                                                                   tCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCctcg                               <  1:5990347/39‑1 (MQ=255)
                                                                   tCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCcgcg                               <  1:1668208/39‑1 (MQ=255)
                                                                   |                                                                   
AGAGATTTCAGTGTCGGAGGACGGTAAACCCAGTTGTTCACGTTTGGCTTCCAGCGCATGCATGGCGTCAAAGATAACTTTGTCGATATCCGCCCCCGAGGCCGCGCTGACCGTTGCTGCCAGTGTACCTTCGAC  >  AP009048/1250257‑1250391

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: