Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 136,262 G→T 8.0% A26A (GCC→GCA yadK ← predicted fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported136,2620GT 8.0% 247.4 / 5.6 75A26A (GCC→GCAyadKpredicted fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (39/30);  new base T (4/2);  total (43/32)
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

GTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAG  >  Exported/136193‑136329
                                                                     |                                                                   
gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGTAACAAGGTCAATATCCTGTCCAGCAATgg                                                                    <  1:2031585/71‑1 (MQ=255)
gTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATg                                                                     <  1:324417/70‑1 (MQ=255)
 ttACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGc                                                                   <  1:2589782/71‑1 (MQ=255)
         tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATgg                                                                    <  1:3234884/62‑1 (MQ=255)
         tCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATgg                                                                    <  1:813199/62‑1 (MQ=255)
            gCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGa                                                        >  1:1664984/1‑71 (MQ=255)
              ttttGCAGGTGCTGTTTTTTACATTGGCAATAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCg                                                         <  1:2293110/68‑1 (MQ=255)
              ttttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                       >  1:3072456/1‑70 (MQ=255)
              ttttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                      >  1:3146144/1‑71 (MQ=255)
              ttttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                      >  1:1878671/1‑71 (MQ=255)
                ttGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGc                                                     >  1:1298106/1‑70 (MQ=255)
                  gCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaa                                                   >  1:2832105/1‑70 (MQ=255)
                   cAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaaa                                                 <  1:1974341/71‑1 (MQ=255)
                    aGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGa                                                               <  1:816976/56‑1 (MQ=255)
                         cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAATCAGTaa                                           <  1:474191/71‑1 (MQ=255)
                         cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                           >  1:2257454/1‑71 (MQ=255)
                         cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                           <  1:1929242/71‑1 (MQ=255)
                         cTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTa                                            <  1:3296452/70‑1 (MQ=255)
                           gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                 <  1:3165427/47‑1 (MQ=255)
                           gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                 <  1:3060915/47‑1 (MQ=255)
                           gTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                 <  1:1872986/47‑1 (MQ=255)
                            ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                       >  1:2844503/1‑56 (MQ=255)
                            ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGAt                                                       >  1:2118535/1‑56 (MQ=255)
                            ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                      >  1:1980010/1‑57 (MQ=255)
                            ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATg                                         <  1:3647946/70‑1 (MQ=255)
                            ttttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGa                                        >  1:2227807/1‑71 (MQ=255)
                             tttttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGa                                        >  1:825858/1‑70 (MQ=255)
                                ttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAg                                                            >  1:236347/1‑47 (MQ=255)
                                ttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCaaa                                                  <  1:1809708/57‑1 (MQ=255)
                                ttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACa                                               >  1:2537597/1‑60 (MQ=255)
                                ttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacag                                     >  1:14904/1‑70 (MQ=255)
                                ttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacag                                     >  1:1337246/1‑70 (MQ=255)
                                ttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacag                                     >  1:1426440/1‑70 (MQ=255)
                                ttACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacag                                     >  1:563596/1‑70 (MQ=255)
                                  aCATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTa                                            <  1:3678544/61‑1 (MQ=255)
                                   cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTg                                                                 >  1:140765/1‑39 (MQ=255)
                                   cATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaa                                  >  1:1393415/1‑70 (MQ=255)
                                      tGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTaa                                           <  1:1224425/58‑1 (MQ=255)
                                        gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAat                            >  1:251632/1‑71 (MQ=255)
                                        gCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAat                            >  1:1760066/1‑71 (MQ=255)
                                          aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACCGAatta                          <  1:474922/71‑1 (MQ=255)
                                          aacaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAatta                          >  1:3522967/1‑71 (MQ=255)
                                           acaaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattat                         >  1:172814/1‑71 (MQ=255)
                                            caaGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatt                        <  1:390851/71‑1 (MQ=255)
                                              aGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAattatta                       >  1:2902860/1‑69 (MQ=255)
                                               ggTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaaca                                <  1:2448308/60‑1 (MQ=255)
                                               ggTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaaca                                <  1:2909496/60‑1 (MQ=255)
                                                 tCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacaga                              >  1:3414703/1‑60 (MQ=255)
                                                 tCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                   >  1:1595349/1‑71 (MQ=255)
                                                  cAATATCCTGTCCAGCAATGGCTGGACAGGCGATg                                                      >  1:644971/1‑35 (MQ=255)
                                                  cAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacag                               <  1:955396/58‑1 (MQ=255)
                                                  cAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacag                               <  1:1754258/58‑1 (MQ=255)
                                                  cAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                   >  1:2992987/1‑70 (MQ=255)
                                                  cAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCtt                   >  1:1881557/1‑70 (MQ=255)
                                                   aaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacagaa                             >  1:2301180/1‑59 (MQ=255)
                                                   aaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                  <  1:1364588/70‑1 (MQ=255)
                                                   aaTATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                  <  1:2213819/70‑1 (MQ=255)
                                                    atatCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcat                >  1:3084119/1‑71 (MQ=255)
                                                    atatCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcat                <  1:346318/71‑1 (MQ=255)
                                                          tGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGacagaaacag                               <  1:257604/50‑1 (MQ=255)
                                                           gTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcatc               >  1:1912543/1‑65 (MQ=255)
                                                           gTCCAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtcat                >  1:586292/1‑64 (MQ=255)
                                                           gTCCAGCAATGGCTGGCCAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCtt          >  1:2524142/1‑70 (MQ=255)
                                                           gTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTa         >  1:946270/1‑71 (MQ=255)
                                                             ccAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                    <  1:1076981/58‑1 (MQ=39)
                                                             ccAGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTctct                    <  1:3724718/58‑1 (MQ=39)
                                                              cAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGt      >  1:2047619/1‑71 (MQ=255)
                                                               aGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                  >  1:1224322/1‑58 (MQ=39)
                                                               aGCCATTGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtc                  >  1:1120112/1‑58 (MQ=39)
                                                               aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt     >  1:3478640/1‑71 (MQ=255)
                                                               aGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGt      >  1:1939872/1‑70 (MQ=255)
                                                               aGCAATGGCGGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTtca                 <  1:2496193/59‑1 (MQ=255)
                                                                 cAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGtt     >  1:3532899/1‑69 (MQ=255)
                                                                 cAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa   <  1:1624073/71‑1 (MQ=255)
                                                                 cAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa   <  1:1462950/71‑1 (MQ=255)
                                                                 cAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa   <  1:3051277/71‑1 (MQ=255)
                                                                  aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGa   <  1:3242786/70‑1 (MQ=255)
                                                                  aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAg  >  1:2261140/1‑71 (MQ=255)
                                                                  aaTGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAg  <  1:631974/71‑1 (MQ=255)
                                                                     |                                                                   
GTTACTGATTCCGCTTTTGCAGGTGCTGTTTTTTACATTGGCAACAAGGTCAATATCCTGTCCAGCAATGGCTGGACAGGCGATGCAAAACAGTAATGACAGAAACAGAATTATTCTCTTCATCAGCTTACGTTGAG  >  Exported/136193‑136329

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: