Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,688,183 Δ1 bp 7.4% coding (783/2061 nt) pilJ ← twitching motility protein PilJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,688,1830A. 7.4% 166.2 / 7.0 54coding (783/2061 nt)pilJtwitching motility protein PilJ
Reads supporting (aligned to +/- strand):  ref base A (24/26);  new base . (2/2);  total (26/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 = 9.54e-01

GCCAGCCAGGTTCTCGAAACCGTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGG  >  NC_002947/5688042‑5688322
                                                                                                                                             |                                                                                                                                           
gccagccaGGTTCTCGAAACCGTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACt                                                                                                                                      >  8:269071/1‑149 (MQ=255)
    gccagGTTCTCGAAACCGTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGaa                                                                                                                                  <  5:144870/149‑1 (MQ=255)
     ccagGTTCTCGAAACCGTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGaa                                                                                                                                  >  6:297544/1‑148 (MQ=255)
            ctcGAAACCGTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCg                                                                                                                           <  1:19258/148‑1 (MQ=255)
             tcGAAACCGTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGc                                                                                                                         >  4:255417/1‑149 (MQ=255)
               gAAACCGTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGa                                                                                                                       >  8:287888/1‑149 (MQ=255)
                    cGTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACAGAACCGGCAACGAACTGGAACAACTCGGCGATTTCg                                                                                                                  >  2:300038/1‑149 (MQ=255)
                     gTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCgg                                                                                                                 <  5:290321/149‑1 (MQ=255)
                         gCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGCTTTCGGCCa                                                                                                              <  7:287888/148‑1 (MQ=255)
                          ccAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTgg                                                                                                                                    >  2:97489/1‑125 (MQ=255)
                          ccAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTgg                                                                                                                                    <  1:97489/125‑1 (MQ=255)
                             ggTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGC‑ACGAACTGGAACAACTCGGCGATTTCGGCCAGCcgcg                                                                                                        >  3:287978/1‑149 (MQ=255)
                             ggTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGC‑ACGAACTGGAACAACTCGGCGATTTCGGCCAGCcgcg                                                                                                        >  3:292815/1‑149 (MQ=255)
                                        tcgtcgAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTc                                                                                              >  3:143985/1‑149 (MQ=255)
                                          gtcgAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCg                                                                                             >  6:65994/1‑148 (MQ=255)
                                          gtcgAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCg                                                                                             <  5:65994/148‑1 (MQ=255)
                                          gtcgAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCg                                                                                             <  5:57002/148‑1 (MQ=255)
                                              aGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTc                                                                                        <  5:11247/149‑1 (MQ=255)
                                                  gggTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGaa                                                                                                                                        >  1:190083/1‑97 (MQ=255)
                                                  gggTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGaa                                                                                                                                        <  2:190083/97‑1 (MQ=255)
                                                   ggTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCt                                                                                      >  5:31057/1‑146 (MQ=255)
                                                   ggTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCt                                                                                      <  6:31057/146‑1 (MQ=255)
                                                   ggTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGa                                                                                   >  5:291404/1‑149 (MQ=255)
                                                      ttGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGAcgc                                                                                >  1:86367/1‑149 (MQ=255)
                                                        gCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGAcgcg                                                                               <  2:154447/148‑1 (MQ=255)
                                                         cGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCgg                                                                              <  5:185550/148‑1 (MQ=255)
                                                                cTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTg                                                                      >  6:24855/1‑149 (MQ=255)
                                                                   aaGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGa                                                                    <  4:163192/148‑1 (MQ=255)
                                                                   aaGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGAt                                                                   >  6:162748/1‑149 (MQ=255)
                                                                    aGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGAtt                                                                  >  7:159208/1‑149 (MQ=255)
                                                                     gATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCTAACCGGC‑ACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGt                                                                <  4:287978/149‑1 (MQ=255)
                                                                     gATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACGGGC‑ACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGt                                                                <  4:292815/149‑1 (MQ=255)
                                                                        gTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGa                                                                    >  6:327117/1‑143 (MQ=255)
                                                                        gTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGa                                                                    <  5:327117/143‑1 (MQ=255)
                                                                           gCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCa                                                            <  6:291404/148‑1 (MQ=255)
                                                                            cccGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGAcgcg                                                                               <  4:59405/128‑1 (MQ=255)
                                                                            cccGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGAcgcg                                                                               >  3:59405/1‑128 (MQ=255)
                                                                            cccGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCAtt                                                          >  4:93943/1‑149 (MQ=255)
                                                                              cgcgGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGc                                                        <  7:269071/149‑1 (MQ=255)
                                                                                         gAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATc                                                <  4:143985/146‑1 (MQ=255)
                                                                                         gAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATg                                             >  3:150683/1‑149 (MQ=255)
                                                                                          aaCACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGc                                            >  4:52473/1‑149 (MQ=255)
                                                                                           acacGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGcc                                           <  3:255417/149‑1 (MQ=255)
                                                                                                      ttCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGt                                                                <  8:183385/117‑1 (MQ=255)
                                                                                                      ttCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGt                                                                >  7:183385/1‑117 (MQ=255)
                                                                                                                  ttCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTGCGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCa                                      >  1:229193/1‑131 (MQ=255)
                                                                                                                  ttCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCa                                      <  2:229193/131‑1 (MQ=255)
                                                                                                                              aTCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGc        >  3:50217/1‑149 (MQ=255)
                                                                                                                               tCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGcc       <  5:24855/149‑1 (MQ=255)
                                                                                                                                ccACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCaaa                         >  7:220141/1‑130 (MQ=255)
                                                                                                                                ccACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCaaa                         <  8:220141/130‑1 (MQ=255)
                                                                                                                                  aCCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGaa    <  2:86367/149‑1 (MQ=255)
                                                                                                                                   ccGAACAGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAg   <  4:150683/149‑1 (MQ=255)
                                                                                                                                    cGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAgg  <  5:162748/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
GCCAGCCAGGTTCTCGAAACCGTTCGCCAGGTGCGACGCTTCGTCGAGCAGGGTTTGCGAAAGGCTGAAGATGTTGCCCGCGGCCTCTCGAACACGGAACAGTTCGGGCGAGGTTTCGAGGATTTCATCCACCGAACCGGCAACGAACTGGAACAACTCGGCGATTTCGGCCAGCCGCGCACGCGCGTCGGCGTCCTCGACGCGGGTGACCTGGATTGTCGCATTGCCATTGATCATGCCCTCCAGCACCTGGCCAAAGCGCCCGGCGTCACGGCCGAAGG  >  NC_002947/5688042‑5688322

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