Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 4,523,210 C→T T87I (ACC→ATC)  fimE → tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6064,523,2100CT99.4% 690.7 / ‑4.2 174T87I (ACC→ATC) fimEtyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (75/98);  total (76/98)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

Only 100 of 182 total aligned reads displayed.

TCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGC  >  REL606/4523136‑4523281
                                                                          |                                                                       
 cGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGAtcc                                                                       >  2:178923/1‑76 (MQ=255)
   ccGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAg                                                                     <  2:1485179/76‑1 (MQ=255)
    cGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAgg                                                                    <  1:1101054/76‑1 (MQ=255)
    cGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAgg                                                                    >  1:932383/1‑76 (MQ=255)
     gACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGa                                                                   <  1:4069202/76‑1 (MQ=255)
     gACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGa                                                                   <  1:4825301/76‑1 (MQ=255)
      aCTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGaa                                                                  >  2:3983309/1‑76 (MQ=255)
       cTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAAc                                                                 >  1:623662/1‑76 (MQ=255)
        tGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACg                                                                >  1:2708494/1‑76 (MQ=255)
         gaagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGt                                                               >  1:1562721/1‑76 (MQ=255)
         gaagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGt                                                               <  2:170495/76‑1 (MQ=255)
         gaagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGt                                                               <  1:2074508/76‑1 (MQ=255)
         gaagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGt                                                               <  1:4929628/76‑1 (MQ=255)
         gaagaaCGGATTTTCAACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGt                                                               <  1:1672444/76‑1 (MQ=255)
          aagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTg                                                              <  2:771141/76‑1 (MQ=255)
          aagaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTg                                                              >  2:2923987/1‑76 (MQ=255)
           agaaCGGATTTTCTACCGTTCACCCGTTACGTTTTGCTGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGc                                                             >  1:4513174/1‑76 (MQ=255)
             aaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTa                                                           <  2:1065207/76‑1 (MQ=255)
             aaCGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTa                                                           >  1:4683917/1‑76 (MQ=255)
                ggATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACt                                                        <  2:2181234/76‑1 (MQ=255)
                 gATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTg                                                       <  2:1397518/76‑1 (MQ=255)
                 gATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTg                                                       >  1:6113413/1‑76 (MQ=255)
                  aTTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTgg                                                      >  2:3038680/1‑76 (MQ=255)
                  aTTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTgg                                                      >  2:2610069/1‑76 (MQ=255)
                  aTTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTgg                                                      <  2:5658246/76‑1 (MQ=255)
                  aTTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTgg                                                      <  1:1375589/76‑1 (MQ=255)
                  aTTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTgg                                                      <  1:1937785/76‑1 (MQ=255)
                   ttttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGa                                                     <  2:4946644/76‑1 (MQ=255)
                    tttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGaa                                                    <  1:1533625/76‑1 (MQ=255)
                    tttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGaa                                                    <  1:2327175/76‑1 (MQ=255)
                     ttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGaaa                                                   <  2:29582/76‑1 (MQ=255)
                     ttCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGaaa                                                   <  2:5001568/76‑1 (MQ=255)
                       cTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAgg                                                 >  1:1430646/1‑76 (MQ=255)
                       cTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAgg                                                 >  1:3314600/1‑76 (MQ=255)
                       cTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAgg                                                 >  1:774376/1‑76 (MQ=255)
                           cGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCt                                             <  1:4422971/76‑1 (MQ=255)
                           cGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCt                                             >  2:4745821/1‑76 (MQ=255)
                           cGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCt                                             <  1:3148309/76‑1 (MQ=255)
                           cGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCt                                             <  1:5733054/76‑1 (MQ=255)
                            gTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTg                                            <  1:3975190/76‑1 (MQ=255)
                            gTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTg                                            >  2:198772/1‑76 (MQ=255)
                            gTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTg                                            <  2:589232/76‑1 (MQ=255)
                            gTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTg                                            <  2:847676/76‑1 (MQ=255)
                             ttCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGa                                           >  1:1932645/1‑76 (MQ=255)
                             ttCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGa                                           <  2:1227210/76‑1 (MQ=255)
                              tCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGAc                                          >  2:2293809/1‑76 (MQ=255)
                              tCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGAc                                          >  2:5270872/1‑76 (MQ=255)
                               cACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGAcc                                         >  1:1203745/1‑76 (MQ=255)
                                aCCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCg                                        <  2:3398080/76‑1 (MQ=255)
                                aCCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCg                                        <  1:5857829/76‑1 (MQ=255)
                                 cccGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCgg                                       <  1:1786600/76‑1 (MQ=255)
                                  ccGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCggc                                      >  2:4698678/1‑75 (MQ=255)
                                  ccGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGa                                      >  2:236978/1‑76 (MQ=255)
                                   cGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGAc                                     <  2:2169598/76‑1 (MQ=255)
                                   cGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACATGCTAACTGGAAAGGCGCTGACCGGAc                                     >  2:1512834/1‑76 (MQ=255)
                                    gTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACt                                    <  2:652292/76‑1 (MQ=255)
                                     ttACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTg                                   <  2:3213584/76‑1 (MQ=255)
                                      tACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGa                                  >  2:119542/1‑76 (MQ=255)
                                      tACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGa                                  <  2:1908973/76‑1 (MQ=255)
                                      tACGTGTTGAGGAGCGTGAAGCCGTGTAGCGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCAGACTGa                                  <  2:2464168/76‑1 (MQ=255)
                                        cGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACg                                <  1:2616303/76‑1 (MQ=255)
                                         gTTTTGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGc                               >  2:6033731/1‑76 (MQ=255)
                                          ttttGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCt                              >  1:4325755/1‑76 (MQ=255)
                                          tggtgacgaGCGGGAAGCCGTGGAGCGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCt                              <  2:983303/73‑1 (MQ=255)
                                           tttGATGATCTTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCta                             <  2:2692166/76‑1 (MQ=255)
                                           tttGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCta                             <  2:5161245/76‑1 (MQ=255)
                                           tttGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCta                             <  1:3232971/76‑1 (MQ=255)
                                           tttGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCta                             >  1:5283388/1‑76 (MQ=255)
                                            ttGATGAGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCtat                            >  2:868255/1‑76 (MQ=255)
                                             tgatgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCtata                           >  1:1289915/1‑76 (MQ=255)
                                             tgatgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCtata                           <  1:852771/76‑1 (MQ=255)
                                             tgatgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGCGCTAACTGGAAAGGCGCTGACCGGACTGACGCtata                           <  2:5754285/76‑1 (MQ=255)
                                              gatgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTatat                          <  2:2663659/76‑1 (MQ=255)
                                              gatgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTatat                          >  2:1314492/1‑76 (MQ=255)
                                              gatgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTatat                          <  1:5267667/76‑1 (MQ=255)
                                               atgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAACGGCGCTGACCGGACTGACGTGAtatt                         >  2:3656548/1‑76 (MQ=255)
                                               atgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtatt                         <  2:679969/76‑1 (MQ=255)
                                               atgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtatt                         >  2:238983/1‑76 (MQ=255)
                                                tgaGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtattt                        >  2:2838315/1‑76 (MQ=255)
                                                  aGCGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTAtatttat                      >  1:3521247/1‑76 (MQ=255)
                                                    cGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATatttattt                    >  2:42927/1‑76 (MQ=255)
                                                    cGTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATatttattt                    <  1:6125754/76‑1 (MQ=255)
                                                     gTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtc                   >  1:4882150/1‑76 (MQ=255)
                                                     gTGAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtc                   >  2:6278655/1‑76 (MQ=255)
                                                       gAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtctc                 <  2:4380911/76‑1 (MQ=255)
                                                       gAAGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTtctc                 <  1:5986091/76‑1 (MQ=255)
                                                        aaGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcg                <  2:6247333/76‑1 (MQ=255)
                                                        aaGCCGTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTcg                <  2:5443040/76‑1 (MQ=255)
                                                             gTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCcgcg           >  1:1059976/1‑76 (MQ=255)
                                                             gTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCcgcg           >  2:3559159/1‑76 (MQ=255)
                                                             gTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCcgcg           <  1:3663440/76‑1 (MQ=255)
                                                             gTGGAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTAGCcgcg           >  2:1655495/1‑76 (MQ=255)
                                                                gAACGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGa        <  1:104892/76‑1 (MQ=255)
                                                                 aaCGCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAg       >  2:4856233/1‑76 (MQ=255)
                                                                    gCTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCg    >  2:4133773/1‑76 (MQ=255)
                                                                     cTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCgg   >  1:684349/1‑76 (MQ=255)
                                                                     cTGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCgg   >  1:3579770/1‑76 (MQ=255)
                                                                      tGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGc  <  2:1371970/76‑1 (MQ=255)
                                                                      tGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGc  >  1:5392646/1‑76 (MQ=255)
                                                                      tGGATCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGc  >  1:2213882/1‑76 (MQ=255)
                                                                          |                                                                       
TCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGGAAAGGCGCTGACCGGACTGACGCTATATTTATTTCTCGCCGCGGGAGTCGGC  >  REL606/4523136‑4523281

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