Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AP009048 4,119,900 C→T 6.2% T63I (ACA→ATA) ‡ hofQ → predicted fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AP0090484,119,9000CT 6.2% 513.8 / 15.7 160T63I (ACA→ATA) ‡hofQpredicted fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base C (68/81);  new base T (4/6);  total (73/87)
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.19e-01

Only 100 of 160 total aligned reads displayed.

GCTGAACAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAGGAAGG  >  AP009048/4119833‑4119969
                                                                   |                                                                     
gctgAACAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACa                                                                      <  1:5901246/69‑1 (MQ=255)
gctgAACAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAg                                                                     <  1:1161132/70‑1 (MQ=255)
   gAACAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGt                                                                 >  1:1497167/1‑71 (MQ=255)
    aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGt                                                                 <  1:5784882/70‑1 (MQ=255)
    aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATg                                                                  <  1:3941294/69‑1 (MQ=255)
    aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGtt                                                                >  1:1874807/1‑71 (MQ=255)
    aaCAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGtt                                                                >  1:1864862/1‑71 (MQ=255)
    aaCAGGAGAAGTTGAACCTGGGCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGt                                                                 <  1:4221674/70‑1 (MQ=255)
      cAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTTTCGTTACATCTAACAGATGTTcc                                                              >  1:3838553/1‑71 (MQ=255)
        ggAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            >  1:1028221/1‑71 (MQ=255)
        ggAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            <  1:4253403/71‑1 (MQ=255)
        ggAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            >  1:3831095/1‑71 (MQ=255)
         gagaAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTg                                                           <  1:2604982/71‑1 (MQ=255)
          agaAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGtt                                                                >  1:2891174/1‑65 (MQ=255)
             aGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGtaa                                                               >  1:6110285/1‑61 (MQ=38)
             aGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGtaa                                                               >  1:397929/1‑61 (MQ=38)
             aGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGtaa                                                               >  1:2260818/1‑61 (MQ=38)
              gTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGt                                                                 >  1:5680683/1‑60 (MQ=255)
              gTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGc                                                      <  1:3659028/71‑1 (MQ=255)
               ttGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            <  1:884058/64‑1 (MQ=255)
               ttGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCt                                                            <  1:6230181/64‑1 (MQ=255)
               ttGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCa                                                     <  1:4918012/71‑1 (MQ=255)
                tGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGc                                                      >  1:2770888/1‑69 (MQ=255)
                tGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAg                                                    >  1:740456/1‑71 (MQ=255)
                tGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAg                                                    >  1:7923536/1‑71 (MQ=255)
                 gAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTcc                                                              <  1:833801/60‑1 (MQ=255)
                    ccTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTgg                                                          <  1:269478/61‑1 (MQ=255)
                    ccTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAg                                                    <  1:7860940/67‑1 (MQ=255)
                     cTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGa                                                                    <  1:1276064/50‑1 (MQ=255)
                     cTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCAc                                                >  1:5715029/1‑70 (MQ=255)
                      tGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGt                                                                 <  1:7813468/52‑1 (MQ=37)
                      tGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGt                                                                 <  1:4596050/52‑1 (MQ=37)
                      tGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTGTTAGATGt                                                                 <  1:1226900/52‑1 (MQ=37)
                      tGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAgg                                                   >  1:1157202/1‑66 (MQ=255)
                      tGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCa                                                 >  1:1923321/1‑68 (MQ=255)
                         tCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACa                                            >  1:5075195/1‑70 (MQ=255)
                         tCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACa                                            >  1:7491781/1‑70 (MQ=255)
                              tCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGc                                                      <  1:4826251/55‑1 (MQ=255)
                              tCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACa                                            >  1:4606978/1‑65 (MQ=255)
                              tCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACa                                            >  1:2759606/1‑65 (MQ=255)
                              tCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       >  1:5161784/1‑70 (MQ=255)
                              tCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       >  1:7203071/1‑70 (MQ=255)
                               cGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGt                                      <  1:2603239/70‑1 (MQ=255)
                                gCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCAc                                                <  1:6789828/59‑1 (MQ=255)
                                gCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCAc                                                <  1:5299571/59‑1 (MQ=255)
                                  cAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGt                                   >  1:7635062/1‑70 (MQ=255)
                                  cAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTg                                  <  1:5089195/71‑1 (MQ=255)
                                    gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:830633/71‑1 (MQ=255)
                                    gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:4063159/71‑1 (MQ=255)
                                    gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:4401420/71‑1 (MQ=255)
                                    gACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:6614369/71‑1 (MQ=255)
                                     aCGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaa                               <  1:6460168/71‑1 (MQ=255)
                                      cGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              <  1:2362923/71‑1 (MQ=255)
                                        tCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAgc                            <  1:6537090/71‑1 (MQ=255)
                                        tCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAgc                            <  1:6168190/71‑1 (MQ=255)
                                          aGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGt                                      <  1:7919428/59‑1 (MQ=255)
                                          aGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGt                                      <  1:4065147/59‑1 (MQ=255)
                                          aGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAgcgc                          >  1:6279705/1‑71 (MQ=255)
                                              gTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAGTTTGTAGTgc                                 >  1:2045013/1‑59 (MQ=38)
                                               tACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCa                                                 >  1:3522496/1‑43 (MQ=255)
                                               tACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCa                                                 >  1:6546734/1‑43 (MQ=255)
                                               tACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGAc                     >  1:5566265/1‑71 (MQ=255)
                                               tACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGAc                     >  1:1677194/1‑71 (MQ=255)
                                                aCGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAg                                    >  1:7811810/1‑55 (MQ=255)
                                                aCGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACt                    >  1:2204757/1‑71 (MQ=255)
                                                aCGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACt                    >  1:5560762/1‑71 (MQ=255)
                                                 cGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCAc                                                >  1:6502240/1‑42 (MQ=255)
                                                 cGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCAc                                                >  1:5125792/1‑42 (MQ=255)
                                                 cGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:3122338/58‑1 (MQ=255)
                                                 cGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:2058481/58‑1 (MQ=255)
                                                 cGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                <  1:3955342/58‑1 (MQ=255)
                                                  ggTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGa                  >  1:1720993/1‑71 (MQ=255)
                                                   gtgtCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGa                  >  1:6314149/1‑70 (MQ=255)
                                                    tgtCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       <  1:5326709/48‑1 (MQ=255)
                                                      tCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACt                                        <  1:7280521/45‑1 (MQ=255)
                                                      tCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACt                                        <  1:5644004/45‑1 (MQ=255)
                                                      tCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              >  1:6985714/1‑55 (MQ=255)
                                                      tCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              >  1:5397204/1‑55 (MQ=255)
                                                       cGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGa                  >  1:6986587/1‑66 (MQ=255)
                                                        gTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAAcgc            >  1:7293633/1‑71 (MQ=255)
                                                        gTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAAcgc            <  1:2874736/71‑1 (MQ=255)
                                                        gTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAAcg             <  1:2598065/70‑1 (MQ=255)
                                                        gTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAAcg             <  1:3702694/70‑1 (MQ=255)
                                                         ttACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       <  1:7683829/43‑1 (MQ=255)
                                                         ttACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       <  1:2563357/43‑1 (MQ=255)
                                                         ttACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTg                                       <  1:4173257/43‑1 (MQ=255)
                                                          tACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaa                                >  1:6783838/1‑49 (MQ=255)
                                                           aCATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCgg          >  1:1615460/1‑70 (MQ=255)
                                                           aCATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCgg          >  1:7548555/1‑70 (MQ=255)
                                                            cATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              >  1:1062207/1‑49 (MQ=255)
                                                            cATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              >  1:1901965/1‑49 (MQ=255)
                                                             aTCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAg       <  1:2610923/71‑1 (MQ=255)
                                                               cTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGaaaa                              <  1:471168/46‑1 (MQ=255)
                                                               cTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAg                             <  1:4261020/47‑1 (MQ=255)
                                                               cTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAgg      >  1:1394960/1‑70 (MQ=255)
                                                               cTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAgg      >  1:4145598/1‑70 (MQ=255)
                                                                  aCAGATGTTCCCTGGAAGCGGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAGGaa    <  1:5840770/69‑1 (MQ=255)
                                                                  aCAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAGGaa    <  1:31473/69‑1 (MQ=255)
                                                                  aCAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAGGa     <  1:1393200/68‑1 (MQ=255)
                                                                  aCAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAGGAAgg  <  1:1352109/71‑1 (MQ=255)
                                                                   |                                                                     
GCTGAACAGGAGAAGTTGAACCTGGTCGTGTCGCCAGACGTCAGCGGTACGGTGTCGTTACATCTAACAGATGTTCCCTGGAAGCAGGCACTACAAACTGTAGTGAAAAGCGCCGGACTGATAACGCGGCAGGAAGG  >  AP009048/4119833‑4119969

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: