Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 1,318,494 C→A 8.1% R29L (CGT→CTT)  cld ← regulator of length of O‑antigen component of lipopolysaccharide chains

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported1,318,4940CA 8.1% 257.0 / 5.3 74R29L (CGT→CTT) cldregulator of length of O‑antigen component of lipopolysaccharide chains
Reads supporting (aligned to +/- strand):  ref base C (44/24);  new base A (3/3);  total (47/27)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTAT  >  Exported/1318424‑1318563
                                                                      |                                                                     
cGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccac                                                                       >  1:4304901/1‑71 (MQ=255)
cGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccac                                                                       >  1:1560906/1‑71 (MQ=255)
    aCTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccactcca                                                                   >  1:545137/1‑71 (MQ=255)
    aCTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccacgcca                                                                   >  1:123203/1‑71 (MQ=255)
    aCTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccacgcc                                                                    >  1:1039777/1‑70 (MQ=255)
     cTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccacgccac                                                                  <  1:3891506/71‑1 (MQ=255)
     cTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccacgcc                                                                    <  1:1088641/69‑1 (MQ=255)
     cTGCCAAATATCCAATAGCCAGGGCAATACCCACAATGACGGAAATGATGATTGTCATCTTgccacgcca                                                                   >  1:1171043/1‑70 (MQ=255)
       gCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:103777/1‑71 (MQ=255)
        ccAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTgccacgccac                                                                  >  1:625005/1‑68 (MQ=255)
        ccAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCcaca                                                                 <  1:1990516/69‑1 (MQ=255)
          aaaTATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAAc                                                               >  1:581977/1‑69 (MQ=255)
          aaaTATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTg                                                             <  1:1222752/71‑1 (MQ=255)
          aaaTATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTg                                                             <  1:2886021/71‑1 (MQ=255)
          aaaTATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTg                                                             <  1:3513195/71‑1 (MQ=255)
           aaTATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGc                                                            <  1:2302326/71‑1 (MQ=255)
           aaTATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGc                                                            <  1:2882755/71‑1 (MQ=255)
             tatCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCAc                                                          >  1:4043087/1‑71 (MQ=255)
                ccAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAg                                                       >  1:1359094/1‑71 (MQ=255)
                ccAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAg                                                       >  1:3869985/1‑71 (MQ=255)
                ccAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAg                                                       >  1:351636/1‑71 (MQ=255)
                   aTAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTa                                                     >  1:2987852/1‑70 (MQ=255)
                   aTAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTa                                                     >  1:7924/1‑70 (MQ=255)
                   aTAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTa                                                     >  1:584109/1‑70 (MQ=255)
                         aGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaat                                              <  1:547538/71‑1 (MQ=255)
                         aGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaat                                              <  1:1000444/71‑1 (MQ=255)
                         aGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaat                                              <  1:3504386/71‑1 (MQ=255)
                         aGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaat                                              <  1:2500913/71‑1 (MQ=255)
                           ggCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatca                                            <  1:3870646/71‑1 (MQ=255)
                           ggCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatca                                            <  1:519675/71‑1 (MQ=255)
                           ggCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatca                                            <  1:335507/71‑1 (MQ=255)
                            gCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatca                                            >  1:1918423/1‑70 (MQ=255)
                            gCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatca                                            >  1:143619/1‑70 (MQ=255)
                            gCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaat                                              <  1:3691540/68‑1 (MQ=255)
                             cAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatcaaa                                          <  1:2599096/71‑1 (MQ=255)
                              aaTAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:2686598/1‑48 (MQ=255)
                              aaTAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:3524413/1‑48 (MQ=255)
                              aaTAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:1106164/1‑48 (MQ=255)
                              aaTAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:3000824/1‑48 (MQ=255)
                              aaTAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:1512760/1‑48 (MQ=255)
                              aaTAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:1156656/1‑48 (MQ=255)
                               aTAGCCACAATGACGGAAATTATGATTGTCATCTTGCCACGCCACAAc                                                               <  1:3760073/48‑1 (MQ=255)
                               aTAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAAc                                                               <  1:329881/48‑1 (MQ=255)
                                 aGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAAc                                                               >  1:2242182/1‑46 (MQ=255)
                                  gCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatcaaatcaat                                     >  1:4090435/1‑71 (MQ=255)
                                  gCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatcaaatcaat                                     >  1:2198823/1‑71 (MQ=255)
                                  gCCACAATGACGGAAATGATGATTGTCATCTTGCCAAGCCACAACTGCACTAGTaaa                                                   <  1:3999543/57‑1 (MQ=255)
                                  gCCACAATGACGGAAATGATGATTGTCATCTTGCCAAGCCACAACTGCACTAGTaaa                                                   <  1:3808384/57‑1 (MQ=255)
                                  gCCACAATGACGGAAATGATGATTGTCATCTTGCCAAGCCACAACTGCACTAGTaaa                                                   <  1:183652/57‑1 (MQ=255)
                                    cacaATGACGGAAATGATGATTGTCATCTTGCCAAGCCACAACTGCACTAGTaaatc                                                 >  1:4399820/1‑57 (MQ=255)
                                    cacaATGACGGAAATGATGATTGTCATCTTGCCAAGCCACAACTGCACTAGTaaatc                                                 >  1:2060300/1‑57 (MQ=255)
                                    cacaATGACGGAAATGATGATTGTCATCTTGCCAAGCCACAACTGCACTAGTaaatc                                                 >  1:2084680/1‑57 (MQ=255)
                                     acaATGACGGAAATGATGATTGTCATCTTgccacgacac                                                                  <  1:3388688/39‑1 (MQ=255)
                                       aaTGACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:2990357/1‑39 (MQ=255)
                                       aaTGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatc                                                 >  1:1301430/1‑54 (MQ=255)
                                       aaTGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatc                                                 >  1:1125125/1‑54 (MQ=255)
                                       aaTGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatc                                                 >  1:3790625/1‑54 (MQ=255)
                                        aTGACGGAAATGATGATTGTCATCTTgccacgccac                                                                  <  1:399146/36‑1 (MQ=255)
                                        aTGACGGAAATGATGATTGTCATCTTgccacgccac                                                                  <  1:3184920/36‑1 (MQ=255)
                                        aTGACGGAAATGATGATTGTCATCTTgccacgccac                                                                  <  1:2190900/36‑1 (MQ=255)
                                        aTGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTc                               <  1:3612028/71‑1 (MQ=255)
                                         tGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatca                                            <  1:647478/57‑1 (MQ=255)
                                         tGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTcc                              >  1:1647223/1‑71 (MQ=255)
                                          gACGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:3905818/1‑36 (MQ=255)
                                           aCGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:3617335/1‑35 (MQ=255)
                                           aCGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:2622471/1‑35 (MQ=255)
                                           aCGGAAATGATGATTGTCATCTTGCCACGCCACaa                                                                >  1:4038462/1‑35 (MQ=255)
                                           aCGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatcaaa                                          >  1:4016544/1‑57 (MQ=255)
                                           aCGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTaaatcaatcaaa                                          >  1:1187040/1‑57 (MQ=255)
                                           aCGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCg                             <  1:855359/70‑1 (MQ=255)
                                           aCGGAAATGATGATTGTCATCTTCCCACGCCACaa                                                                >  1:4091766/1‑35 (MQ=255)
                                                aaTGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCAt                       >  1:3613917/1‑71 (MQ=255)
                                                aaTGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCAt                       >  1:4211749/1‑71 (MQ=255)
                                                   gatgatTGTCATCTTGCCACGCCACAACTGCACTAGTaaa                                                   >  1:2891519/1‑40 (MQ=255)
                                                   gatgatTGTCATCTTGCCACGCCACAACTGCACTAGTaaa                                                   >  1:35966/1‑40 (MQ=255)
                                                     tgatTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGttt                  >  1:11463/1‑71 (MQ=255)
                                                     tgatTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGttt                  >  1:3607107/1‑71 (MQ=255)
                                                           tCATCTTGCCACGCCACAACTGCACTAGTaaatcaatca                                            <  1:1036499/39‑1 (MQ=255)
                                                            cATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGcc              >  1:4292424/1‑68 (MQ=255)
                                                             aTCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGtttt                 >  1:429321/1‑64 (MQ=255)
                                                             aTCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGtttt                 >  1:2475233/1‑64 (MQ=255)
                                                             aTCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGtttt                 >  1:1221077/1‑64 (MQ=255)
                                                             aTCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAg           <  1:4134254/70‑1 (MQ=255)
                                                              tCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTc                         <  1:2510289/55‑1 (MQ=255)
                                                              tCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTc                         <  1:294987/55‑1 (MQ=255)
                                                                ttGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGaaa        >  1:2522999/1‑70 (MQ=255)
                                                                 tGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACa      >  1:2552135/1‑71 (MQ=255)
                                                                  gccacgccacAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGaaa        >  1:4372316/1‑68 (MQ=255)
                                                                  gccacgccacAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGaaa        <  1:1480291/68‑1 (MQ=255)
                                                                  gccacgccacAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGaaa        >  1:577873/1‑68 (MQ=255)
                                                                  gccacgccacAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGaaa        >  1:605695/1‑68 (MQ=255)
                                                                  gccacgccacAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACa      <  1:994618/70‑1 (MQ=255)
                                                                     acgccacAACTGCACTAGTaaatcaatcaaatcaatc                                    <  1:1430100/37‑1 (MQ=255)
                                                                     acgccacAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACattat  <  1:32069/71‑1 (MQ=255)
                                                                      |                                                                     
CGCTACTGCCAAATATCCAATAGCCAGGGCAATAGCCACAATGACGGAAATGATGATTGTCATCTTGCCACGCCACAACTGCACTAGTAAATCAATCAAATCAATCTGTTCCGGGTCATGGTTTTGCCCAGAAACATTAT  >  Exported/1318424‑1318563

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: