Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,181,048 T→C 8.3% P282P (CCA→CCG)  yehA ← putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,181,0480TC 8.3% 129.1 / 2.5 48P282P (CCA→CCG) yehAputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base T (16/28);  new base C (2/2);  total (18/30)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

GCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTCCGTGATGTTTGATTTTGTTGACCTAATAAATTTT  >  CP009273/2180901‑2181193
                                                                                                                                                   |                                                                                                                                                 
gcCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTg                                                                                                                                                  <  1:523087/149‑1 (MQ=255)
 tcagtgttAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTAGCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTgg                                                                                                                                                 <  1:302808/143‑1 (MQ=255)
 ccAGATTTCACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCCGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTgg                                                                                                                                                 <  2:809420/149‑1 (MQ=255)
         aaCTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTGCCTTTAAATTTGTATATCATTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCAt                                                                                                                                         >  1:630102/1‑149 (MQ=255)
         aaCTTTGCTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCAt                                                                                                                                         >  1:593340/1‑149 (MQ=255)
                  gggaTCGTGGACGTAGTAACTGGCTGTTATCGGAATTGAATCCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTATATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCaa                                                                                                                                <  1:392149/146‑1 (MQ=255)
                    tATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGa                                                                                                                              <  1:50149/149‑1 (MQ=255)
                      tCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAAt                                                                                                                            <  2:700349/149‑1 (MQ=255)
                       cGTGGACGTAGTAACTGGCTGTTATCGGAATTGAATCTGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATa                                                                                                                           <  2:539766/149‑1 (MQ=255)
                       cGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATa                                                                                                                           >  2:89318/1‑149 (MQ=255)
                        gTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATaa                                                                                                                          >  2:306079/1‑149 (MQ=255)
                            aCGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATcc                                                                                                                      >  2:149022/1‑149 (MQ=255)
                              gtagtaACTGGCTGTTATCGGTATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATAGCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTa                                                                                                                    <  1:501139/149‑1 (MQ=255)
                                  taACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATAGCCTTGACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGcc                                                                                                                <  2:195232/149‑1 (MQ=255)
                                        gCTTTTATCGGAATTGACACCTTTGTAGCGGTAATTTTTGATGTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATc                                                                                                          <  1:248543/149‑1 (MQ=255)
                                        gCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTCGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATc                                                                                                          <  2:341630/149‑1 (MQ=255)
                                         cTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCt                                                                                                         <  1:757191/149‑1 (MQ=255)
                                          tGTTATCGGAAGTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGGGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCtt                                                                                                        <  2:599263/149‑1 (MQ=255)
                                               tccgAATTCAACCCGTTGTATCTGTAAGTTTTGACTTCGCCGCACCTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCa                                                                                                   <  2:459341/146‑1 (MQ=255)
                                                cgtaaTTGACATCGGTGTATCTGTAAGTTTTGACTTCGCCGAAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCGACGCCATTATCTTCCGCAc                                                                                                  <  2:660515/146‑1 (MQ=255)
                                                 ggAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACCGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCAcc                                                                                                 >  1:821559/1‑149 (MQ=255)
                                                 ggAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACCGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTAACTTCCGCAcc                                                                                                 >  1:451513/1‑149 (MQ=255)
                                                   aaTTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCtt                                                                                               <  1:167403/149‑1 (MQ=255)
                                                    aTTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCttt                                                                                              <  2:40294/149‑1 (MQ=255)
                                                     ttGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCGGTTAGTGATATTAAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCGTACGCCCTTAGCTTCCGCCCTTTTc                                                                                             >  1:260967/1‑149 (MQ=255)
                                                        aaaCCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGt                                                                                          >  1:121373/1‑149 (MQ=255)
                                                        aaaCCGTCGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTTCTATTAGTGATATTTAAAAGTGTGAGTGGAACACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAga                                                                                          >  2:201787/1‑148 (MQ=255)
                                                          aCCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTAt                                                                                        >  2:263394/1‑149 (MQ=255)
                                                          aCCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAACTGGTGTAGCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATGTTCCCCACCTTTAAGtct                                                                                        >  1:361764/1‑147 (MQ=255)
                                                             gTTGTATTTGTAAGTTTTGACTTCGCCGCAACTTCCGTTTAATTTGTATATCCTTTACTGCTGTTAGTGCTATTTAAAAGTGTGACTGGCGATCCATTACTTTCAAGAATAAATCCTACGCCATTTTCTTCCGCACCTTTCAGAATtgt                                                                                     <  1:371030/149‑1 (MQ=255)
                                                                       tAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTCGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAATAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTg                                                                           >  2:752796/1‑149 (MQ=255)
                                                                          gTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGtt                                                                        <  1:715081/149‑1 (MQ=255)
                                                                          gTTTTGACTCCGCCCTCACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGGGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGtt                                                                        <  1:664896/149‑1 (MQ=255)
                                                                               gACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACaa                                                                   <  1:763112/149‑1 (MQ=255)
                                                                                            cTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGAt                                                      >  2:133418/1‑149 (MQ=255)
                                                                                                  ttAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGGATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATac                                                >  2:398028/1‑149 (MQ=255)
                                                                                                             tatCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGACCTTGATAAATACACCTGCAtgtt                                     >  1:86954/1‑149 (MQ=255)
                                                                                                              atCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTc                                    <  2:841901/149‑1 (MQ=255)
                                                                                                              atCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTc                                    <  2:535044/149‑1 (MQ=255)
                                                                                                                  tttACTGCTGTTAGTGCTATAAAAGATGGTGACTGGCGAACCATTACTTTCAAGAATAAATCATACGCCATAATCTTCCCCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAATCAGAACTTGATAAATACACCTGCATTTTCCGTg                                >  1:442637/1‑149 (MQ=255)
                                                                                                                        gctgTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTGCTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAGGACAGGACTTGATAAATACACCTGCATTTTCCGTGATGttt                          <  1:149022/149‑1 (MQ=255)
                                                                                                                         ctgTTAGTGATATTTAAGAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCGGTTTTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTCCGTGATGTTTg                         <  2:319995/149‑1 (MQ=255)
                                                                                                                               gggATATTTACAAGTTGTACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTAGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTCCGTGATGTTTGATTttg                   <  1:564463/147‑1 (MQ=255)
                                                                                                                                       tAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATGAATAACCCTACATGTTCCGTGATGTTTGATTTTGTTGACCta           >  1:828806/1‑149 (MQ=255)
                                                                                                                                               tgACCGGCGATCCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTCCGTGATGTTTGATTTTGTTGACCTAATAAAttt   <  2:451513/149‑1 (MQ=255)
                                                                                                                                               tgACCGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTCCGTGATGTTTGATTTTGTTGACCTAATAAAttt   <  2:821559/149‑1 (MQ=255)
                                                                                                                                                gcgtGGCGTACCATTACTTTCAAGAATAATTCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTCCGTGATGTTTGATTTTGTTGACCTAATAAAtttt  <  1:201787/146‑1 (MQ=255)
                                                                                                                                                gACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTCCGTGATGTTTGATTTTGTTGACCTAATAAAtttt  <  2:121373/149‑1 (MQ=255)
                                                                                                                                                gACTGGCGAACCATTACTTTCAAGAATAAATCCTAAGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTCCGTGATGTTTGATTTTGTTGACCTAATAAAtttt  <  1:398028/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCCAGATTTAACTTTGTTTGTATCGTAGACGTAGTAACTGGCTGTTATCGGAATTGAAACCGTTGTATCTGTAAGTTTTGACTTCGCCGCAACTTCCTTTAAATTTGTATATCCTTTACTGCTGTTAGTGATATTTAAAAGTGTGACTGGCGAACCATTACTTTCAAGAATAAATCCTACGCCATTATCTTCCGCACCTTTCAGTATTGTGTTGCTGTTGGTTAACAAGTCAGAACTTGATAAATACACCTGCATTTTCCGTGATGTTTGATTTTGTTGACCTAATAAATTTT  >  CP009273/2180901‑2181193

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑