Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 1,894,072 C→T 5.6% G293D (GGT→GAT)  yoaA ← conserved hypothetical protein with nucleoside triphosphate hydrolase domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090481,894,0720CT 5.6% 344.6 / 7.0 107G293D (GGT→GAT) yoaAconserved hypothetical protein with nucleoside triphosphate hydrolase domain
Reads supporting (aligned to +/- strand):  ref base C (47/53);  new base T (4/2);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 4.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

Only 100 of 109 total aligned reads displayed.

TTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCAC  >  AP009048/1894002‑1894140
                                                                      |                                                                    
ttCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTAc                                                                      >  1:3551967/1‑71 (MQ=255)
       gtgtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATa                                                               <  1:6562956/71‑1 (MQ=255)
        tgtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcg                                                                    >  1:5674757/1‑65 (MQ=255)
          tCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGAt                                                                >  1:2034727/1‑67 (MQ=255)
          tCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAAc                                                             <  1:226884/70‑1 (MQ=255)
            aTCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCt                                                           <  1:6840635/70‑1 (MQ=255)
            aTCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                          >  1:4168154/1‑71 (MQ=255)
            aTCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                          >  1:1594494/1‑71 (MQ=255)
             tCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTgg                                                         >  1:3214245/1‑71 (MQ=255)
              cGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGc                                                        >  1:5925818/1‑71 (MQ=255)
                    gTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTAc                                                                      <  1:2341815/51‑1 (MQ=255)
                     tAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgcg                                                                  <  1:4232065/54‑1 (MQ=255)
                     tAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGAt                                                                <  1:1047389/56‑1 (MQ=255)
                     tAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGAt                                                                <  1:1413541/56‑1 (MQ=255)
                     tAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                          <  1:707406/62‑1 (MQ=255)
                      aaaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAg                                                >  1:3531643/1‑71 (MQ=255)
                       aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                          >  1:3355634/1‑60 (MQ=255)
                       aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTg                                                          >  1:2014350/1‑60 (MQ=255)
                       aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGc                                                        >  1:7337298/1‑62 (MQ=255)
                       aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGt                                               >  1:2573817/1‑71 (MQ=255)
                       aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGt                                               >  1:849128/1‑71 (MQ=255)
                       aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGt                                               >  1:64905/1‑71 (MQ=255)
                       aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGt                                               >  1:2274406/1‑71 (MQ=255)
                        aaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGtt                                              <  1:7404179/71‑1 (MQ=255)
                         aaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAc                                                    <  1:7164441/64‑1 (MQ=255)
                          aTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTg                                             >  1:5723860/1‑70 (MQ=255)
                          aTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTg                                             >  1:5525724/1‑70 (MQ=255)
                          aTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGc                                            <  1:5793724/71‑1 (MQ=255)
                           tGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCt                                                           <  1:4734079/55‑1 (MQ=255)
                             cccGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGa                                         <  1:5738999/71‑1 (MQ=255)
                             cccGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGa                                         <  1:3602014/71‑1 (MQ=255)
                                gCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAcc                                                   <  1:1972638/58‑1 (MQ=255)
                                gCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACg                                       >  1:543348/1‑70 (MQ=255)
                                gCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACg                                       >  1:918541/1‑70 (MQ=255)
                                  tGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATcgcg                                                                  <  1:2864746/41‑1 (MQ=255)
                                   gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgcg                                                                  <  1:4858716/40‑1 (MQ=255)
                                   gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgcg                                                                  <  1:3799555/40‑1 (MQ=255)
                                   gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAAcc                                                            <  1:1172613/46‑1 (MQ=255)
                                   gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAAcc                                                            <  1:1290441/46‑1 (MQ=255)
                                   gAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGaa                                     <  1:1835964/69‑1 (MQ=255)
                                    aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTATcgcggt                                                                >  1:439278/1‑39 (MQ=39)
                                    aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTATcgcggt                                                                >  1:6463714/1‑39 (MQ=39)
                                    aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTATcgcggt                                                                >  1:7637168/1‑39 (MQ=39)
                                    aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTATcgcggt                                                                >  1:3978366/1‑39 (MQ=39)
                                    aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGaaaa                                   <  1:4708944/70‑1 (MQ=255)
                                     aTTTGCGGATTAGCTAACAGCTCACGCAGGTTACcg                                                                    >  1:4040470/1‑36 (MQ=255)
                                     aTTTGCGGATTAGCTAACAGCTCACGCAGGTTACcg                                                                    >  1:1226960/1‑36 (MQ=255)
                                     aTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATc                                 <  1:244987/71‑1 (MQ=255)
                                      tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGAc                                        >  1:2792172/1‑63 (MQ=255)
                                      tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGaa                                     >  1:7476306/1‑66 (MQ=255)
                                      tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATcc                                >  1:7922475/1‑71 (MQ=255)
                                      tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATcc                                >  1:120683/1‑71 (MQ=255)
                                        tGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCt                               <  1:5990634/70‑1 (MQ=255)
                                         gCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                              <  1:1761756/70‑1 (MQ=255)
                                            gATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCa                          >  1:2312491/1‑71 (MQ=255)
                                             aTTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCAc                         >  1:6120007/1‑71 (MQ=255)
                                              ttAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCAct                        >  1:3604734/1‑71 (MQ=255)
                                               tAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCa                                           >  1:5940023/1‑51 (MQ=255)
                                               tAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCActc                       >  1:2943314/1‑71 (MQ=255)
                                               tAGCTAACAGCTCACGAAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCActc                       >  1:4248/1‑71 (MQ=255)
                                                aGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCActc                       <  1:2014086/70‑1 (MQ=255)
                                                 gCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCActct                      <  1:1716188/70‑1 (MQ=255)
                                                 gCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCActct                      <  1:5337205/70‑1 (MQ=255)
                                                 gCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTg                     <  1:3424302/71‑1 (MQ=255)
                                                  cTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAg                                                <  1:7182468/43‑1 (MQ=255)
                                                  cTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATc                                 <  1:5638870/58‑1 (MQ=255)
                                                    aaCAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCActc                       <  1:3025799/66‑1 (MQ=255)
                                                     aCAGCTCACGCAGGTTAGCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTg                     <  1:5054149/67‑1 (MQ=255)
                                                     aCAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGc                  <  1:3759803/70‑1 (MQ=255)
                                                      cAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCaa                <  1:1191487/71‑1 (MQ=255)
                                                       aGCTTACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                              <  1:5517036/56‑1 (MQ=255)
                                                       aGCTCACGCAGGTTACCGCGATAACCTGGCTCACCg                                                  >  1:614970/1‑36 (MQ=255)
                                                       aGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                              <  1:1245278/56‑1 (MQ=255)
                                                       aGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTg                              <  1:7887156/56‑1 (MQ=255)
                                                        gCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGa              <  1:3010725/71‑1 (MQ=255)
                                                          tCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGa                                      >  1:316608/1‑45 (MQ=255)
                                                          tCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCaa                <  1:3208127/67‑1 (MQ=255)
                                                          tCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGAc             >  1:4304933/1‑70 (MQ=255)
                                                           cACGCAGGTTACCGCGATAACCTGGCTCACGGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGa           >  1:7773749/1‑71 (MQ=255)
                                                            aCGCAGGTTGCCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACg            <  1:2050321/69‑1 (MQ=255)
                                                            aCGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACg            <  1:414878/69‑1 (MQ=255)
                                                            aCGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACg            <  1:4567176/69‑1 (MQ=255)
                                                            aCGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGAt          <  1:351499/71‑1 (MQ=255)
                                                            aCGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGAt          <  1:342174/71‑1 (MQ=255)
                                                            aCGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGAt          <  1:1808657/71‑1 (MQ=255)
                                                             cGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGAt          >  1:3872531/1‑70 (MQ=255)
                                                              gCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCa                          <  1:4790747/53‑1 (MQ=255)
                                                              gCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCa                          <  1:2525572/53‑1 (MQ=255)
                                                               cAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAg       >  1:7023318/1‑71 (MQ=255)
                                                                aGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGa                                         >  1:3129682/1‑36 (MQ=255)
                                                                aGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGa                                         >  1:7654250/1‑36 (MQ=255)
                                                                aGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGa                                         >  1:3502387/1‑36 (MQ=255)
                                                                 ggTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGa                                      <  1:6823287/38‑1 (MQ=255)
                                                                  gTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCAct                        <  1:5896421/51‑1 (MQ=255)
                                                                  gTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGAt          <  1:746515/65‑1 (MQ=255)
                                                                    tACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa        >  1:3030759/1‑65 (MQ=255)
                                                                    tACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCAc  >  1:1579375/1‑71 (MQ=255)
                                                                    tACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCAc  >  1:511649/1‑71 (MQ=255)
                                                                     aCCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa        >  1:6969927/1‑64 (MQ=255)
                                                                     aCCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCAc  <  1:6363165/70‑1 (MQ=255)
                                                                      |                                                                    
TTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCAC  >  AP009048/1894002‑1894140

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: