Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 4,971,387 Δ1 bp 7.7% coding (328/1110 nt) flgI ← Flagellar P‑ring protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24404,971,3870T. 7.7% 136.5 / 7.1 52coding (328/1110 nt)flgIFlagellar P‑ring protein
Reads supporting (aligned to +/- strand):  ref base T (20/25);  new base . (2/2);  total (24/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CGTCAAAGCCACCCACGACCAGGTTGCCCTGGGCAATGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGAATTGCGACAGC  >  KT2440/4971261‑4971511
                                                                                                                              |                                                                                                                            
cGTCAAAGCCACCCACGACCAGGTTGCCCTGGTCAATGGCGTAAACGTTGCCGTAGATACCCTTGAGCGGGGTCATCAACACGCTGACGCCGCGCAGGCTCTTCGAGTTTACGAAAGAAGACACCGGGAtttaaa                                                                                                                      >  6:47676/1‑130 (MQ=255)
cGTCAAAGCCACCCACGACCAGGTTGCCCTGGGCAATGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGCCACGGTGATGTcca                                                                                                                      >  2:26389/1‑135 (MQ=255)
  tCAAAGCCACCCACGACCAGGTTGCCCTGGGCAATGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTccacc                                                                                                                    <  4:304764/135‑1 (MQ=255)
        ccacccacGACCAGGTTGCCCTGGGCAATGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCCTCAACCGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGGCGACACGGGGATGTCCACCACCTgg                                                                                                              >  3:21294/1‑135 (MQ=255)
               cGACCAGGTTGCCCTGGGCAATGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACG‑‑GATGTCCACCACCTGGCCCGGCTTg                                                                                                     >  1:306873/1‑135 (MQ=255)
               cGACCAGGTTGCCCTGGGCAATGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACG‑‑GATGTCCACCACCTGGCCCGGCTTg                                                                                                     >  1:310967/1‑135 (MQ=255)
                 aCCAGGTTGCCCTGGGCAATGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTg                                                                                                     <  2:201314/135‑1 (MQ=255)
                             tGGGCAATGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGGCATCAACCGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGc                                                                                         <  1:304481/135‑1 (MQ=255)
                                  aaTGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTc                                                                                    >  7:262905/1‑135 (MQ=255)
                                   aTGGCGTAACCGTTGCCGTCGATACCCTTGAGCGGGGGCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCg                                                                                   <  5:315287/135‑1 (MQ=255)
                                   aTGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCg                                                                                   <  3:88529/135‑1 (MQ=255)
                                   aTGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGGCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCg                                                                                   <  5:145/135‑1 (MQ=255)
                                       cGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGCCACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCAt                                                                               >  5:318140/1‑135 (MQ=255)
                                       cGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCAt                                                                               >  5:111614/1‑135 (MQ=255)
                                            aCGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCAc                                                                           <  8:60792/134‑1 (MQ=255)
                                                   cGTCGATACCCTTGAGCGGGGTCATCCACAGGCTCCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACCGTGGCAGGTCGGCATGCACCGACACTg                                                                   >  5:305388/1‑135 (MQ=255)
                                                   cGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGCCACTg                                                                   >  4:258270/1‑135 (MQ=255)
                                                    gTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCCCCACCTGGCCCGGCTTGGCGAACGGTGGCCGGTCGGCCTGCCCCGACCCTGc                                                                  >  6:212688/1‑135 (MQ=255)
                                                      cGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCgg                                                                <  6:300992/135‑1 (MQ=255)
                                                      cGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCgg                                                                <  1:76001/135‑1 (MQ=255)
                                                      cGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCgg                                                                <  1:32661/135‑1 (MQ=255)
                                                        aTACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCa                                                              <  6:264391/135‑1 (MQ=255)
                                                           cccTTGAGCGGGGTCATCAACAGGCCGCCGCCGCGCAGGCGCTGCTACTTCCCGATCCACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCATGTCGGCATGCACCGACACTGCGGCAACg                                                           <  5:221084/135‑1 (MQ=255)
                                                               tGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTAGAGTTTCCGATCGACGACCCGGTGATGTCCACCACCAGGCCCGGCTTGGCGCACGGTGGCAGGTCGGCCTGCACCCACCCTGCGGCAACGttct                                                       >  2:15207/1‑135 (MQ=255)
                                                                      ggTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTg                                                <  3:234027/135‑1 (MQ=255)
                                                                         catcaCCAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCCCCACCTGGCCCGGCTTGGCGACCGGTGGCAGGTCGGCATGCACCGACACTGCGGCCACGTTCTTCAGCTGCAc                                             >  2:258359/1‑135 (MQ=255)
                                                                         catcaACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCAc                                             >  8:73907/1‑135 (MQ=255)
                                                                         catcaACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCAc                                             >  3:194792/1‑135 (MQ=255)
                                                                         catcaACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCAc                                             >  3:190773/1‑135 (MQ=255)
                                                                         catcaACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCAc                                             >  2:53385/1‑135 (MQ=255)
                                                                                cggcggcCGCCGCGCAGGCGCTTCGAGTTACCGAGCGACGACACGGGGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCg                                      <  6:32189/134‑1 (MQ=255)
                                                                                     gccgccGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAAcc                                 <  5:28571/135‑1 (MQ=255)
                                                                                     gccgccGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAAcc                                 >  6:28571/1‑135 (MQ=255)
                                                                                     gccgccGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAAAGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGATGCACGTTGCCGGAAcc                                 >  6:178693/1‑135 (MQ=255)
                                                                                       cgccgcGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAAccgg                               >  5:298284/1‑135 (MQ=255)
                                                                                        gccgcGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAAccggc                              <  1:231874/135‑1 (MQ=255)
                                                                                        gccgcGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAAccggc                              <  1:330299/135‑1 (MQ=255)
                                                                                        gccgcGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAAccggc                              <  5:89733/135‑1 (MQ=255)
                                                                                                   tctTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGa                   <  1:269596/135‑1 (MQ=255)
                                                                                                         aGTTACCGATCGACGACACG‑‑GATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGaa            <  2:310967/134‑1 (MQ=255)
                                                                                                         aGTTACCGATCGACGACACG‑‑GATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGaa            <  2:306873/134‑1 (MQ=255)
                                                                                                              ccGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGAATTg         <  3:172468/134‑1 (MQ=255)
                                                                                                               ctagcgacgaCACGGTGATGTCCACCACCTGGCCCGGCTGGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGAATTGCg       <  6:208822/131‑1 (MQ=255)
                                                                                                                gATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATg                                                                              >  7:233705/1‑63 (MQ=255)
                                                                                                                gATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATg                                                                              <  8:233705/63‑1 (MQ=255)
                                                                                                                gATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCCCCTTGATGCCGAATTGCg       >  7:287994/1‑134 (MQ=255)
                                                                                                                gATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGAATTg         <  7:59151/132‑1 (MQ=255)
                                                                                                                gATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGAATTg         >  8:59151/1‑132 (MQ=255)
                                                                                                                gATCGACGACACGGTGATGTCCACCACCTGGACCCGCTTGGCGAACCGTGGCAGGTCGGCATGCACCGACACAGCGGCAACGATCTTCAGCTGCACGTTGCCGGAACCGGCCTGCACCTTGATGCCGAATTGCg       >  7:95046/1‑134 (MQ=255)
                                                                                                                   cgacgaCACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGAATTGCGACAg   <  2:301413/135‑1 (MQ=255)
                                                                                                                     acgacACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGAATTGCGACAGc  <  8:262905/134‑1 (MQ=255)
                                                                                                                     acgacACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGAATTGCGACAGc  <  4:71174/134‑1 (MQ=255)
                                                                                                                              |                                                                                                                            
CGTCAAAGCCACCCACGACCAGGTTGCCCTGGGCAATGGCGTAAACGTTGCCGTCGATACCCTTGAGCGGGGTCATCAACAGGCTGCCGCCGCGCAGGCTCTTCGAGTTACCGATCGACGACACGGTGATGTCCACCACCTGGCCCGGCTTGGCGAACGGTGGCAGGTCGGCATGCACCGACACTGCGGCAACGTTCTTCAGCTGCACGTTGCCGGAACCGGCCGGCACCTTGATGCCGAATTGCGACAGC  >  KT2440/4971261‑4971511

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