Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,462,844 (G)6→7 coding (501/588 nt) gspJ → Type II secretion system protein GspJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,462,8381.G90.0% 160.5 / 9.7 50coding (495/588 nt)gspJType II secretion system protein GspJ
Reads supporting (aligned to +/- strand):  ref base . (2/3);  new base G (16/29);  total (18/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTCTGTGGATGGCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGT  >  NC_000913/3462701‑3462969
                                                                                                                                           |                                                                                                                                      
aCTTCTGTGGATGGCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATgg                                                                                                                                      >  1:589537/1‑140 (MQ=255)
aCTTCTGTGGATGGCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATgg                                                                                                                                      <  1:370026/140‑1 (MQ=255)
 cTTCTGTGGATGGCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATggg                                                                                                                                     <  2:589537/140‑1 (MQ=255)
          aTGGCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACa                                                                                                                          >  2:473775/1‑140 (MQ=255)
            ggCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAgc                                                                                                                        >  2:362742/1‑140 (MQ=255)
            ggCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAgc                                                                                                                        <  1:618372/140‑1 (MQ=255)
            ggCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAgc                                                                                                                        >  2:523770/1‑140 (MQ=255)
             gCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAgcg                                                                                                                       <  1:473775/140‑1 (MQ=255)
              cAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAgcgc                                                                                                                      >  2:411941/1‑140 (MQ=255)
              cAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAgcgc                                                                                                                      >  2:581723/1‑140 (MQ=255)
              cAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAgcgc                                                                                                                      >  1:535237/1‑140 (MQ=255)
              cAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGGTACGTCTGCAGACTCACCAATATGGGGGGGTTACAgcgc                                                                                                                      >  2:294430/1‑140 (MQ=255)
                 gaaggaTCAACCCGCCCCGGCGTTGCTGGAACATGCCGATAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGt                                                                                                                   <  2:203222/140‑1 (MQ=255)
                 gaaggaTCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGt                                                                                                                   <  1:581723/140‑1 (MQ=255)
                  aaggaTCAACCCGCCAAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGtt                                                                                                                  <  2:381747/140‑1 (MQ=255)
                        cAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCAt                                                                                                            <  1:294430/140‑1 (MQ=255)
                        cAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCAt                                                                                                            <  1:601514/140‑1 (MQ=255)
                        cAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCGTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACACCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCAt                                                                                                            <  2:134417/140‑1 (MQ=255)
                        cAACCCGCCCAGGCGATCCTGGTACATGTCGCGAGCTTTTTTCTGGAACGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCAt                                                                                                            <  2:304118/140‑1 (MQ=255)
                             cGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCt                                                                                                       <  2:447349/140‑1 (MQ=255)
                                ccAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAg                                                                                                    <  2:158433/140‑1 (MQ=255)
                                   ggCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAg                                                                                                 <  1:562785/140‑1 (MQ=255)
                                        tGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGc                                                                                            <  2:535237/140‑1 (MQ=255)
                                          cTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCga                                                                                          <  2:637529/140‑1 (MQ=255)
                                            ggAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCgaga                                                                                        <  1:411941/140‑1 (MQ=255)
                                                  tGTCGAGAGCTTTCTTCTGTAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGAAATATGGGGGTGTTACAGCGCCGTTTTGCATTACCTGAGCAGATAGCGAGAGAagag                                                                                  >  2:179981/1‑140 (MQ=255)
                                                  tGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGt                                                                                 >  2:152078/1‑140 (MQ=255)
                                                     cGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGATTACAGCGCCGTGTTGGATTACCTGAGCAGTTAGCGAGAGAAGAGTCg                                                                               >  2:245721/1‑140 (MQ=255)
                                                       agagCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGcc                                                                             >  2:638547/1‑140 (MQ=255)
                                                                    ggAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGc                                                                              >  1:426332/1‑126 (MQ=255)
                                                                    ggAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGc                                                                              <  2:426332/126‑1 (MQ=255)
                                                                      aaaGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGc                                                              <  1:362742/140‑1 (MQ=255)
                                                                       aaGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGcc                                                             <  1:485571/140‑1 (MQ=255)
                                                                       aaGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGcc                                                             <  2:656179/140‑1 (MQ=255)
                                                                        aGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCg                                                            <  2:65469/140‑1 (MQ=255)
                                                                                    gAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAAcc                                                <  2:619036/140‑1 (MQ=255)
                                                                                    gAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAAcc                                                >  2:72954/1‑140 (MQ=255)
                                                                                     aaaGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCa                                               >  2:701652/1‑140 (MQ=255)
                                                                                      aaGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCAt                                              <  1:547398/140‑1 (MQ=255)
                                                                                                         tcAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACgcggc                           >  1:389241/1‑140 (MQ=255)
                                                                                                          cAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACgcggcg                          <  2:385297/140‑1 (MQ=255)
                                                                                                          cAGTCACGTTACATCTGAAGACTCAGCAATATGGGAGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACgcggcg                          <  2:217644/140‑1 (MQ=255)
                                                                                                           aGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACgcggcgc                         >  1:720741/1‑139 (MQ=255)
                                                                                                           aGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGt                         >  1:143415/1‑140 (MQ=255)
                                                                                                               aCGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTcgcg                     <  1:106784/140‑1 (MQ=255)
                                                                                                               aCGTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTcgcg                     <  2:622261/140‑1 (MQ=255)
                                                                                                                 gTTACATCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCgct                   <  1:523770/140‑1 (MQ=255)
                                                                                                                   tACATCTGCAGACTCAGCAATACGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCgctgc                 <  1:245721/140‑1 (MQ=255)
                                                                                                                   tACATCTGCAGACTCAGCAATATGGGGGGCTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCgctgc                 <  1:270357/140‑1 (MQ=255)
                                                                                                                      aTCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCa              <  1:179981/140‑1 (MQ=255)
                                                                                                                      aTCTGCAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCa              <  1:122516/140‑1 (MQ=255)
                                                                                                                           cAGACTCAGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACgcggcg                          >  1:249933/1‑123 (MQ=255)
                                                                                                                           cAGACACAGCAATATGGGGGGGTTACAGCGCAGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACgcggcg                          <  2:249933/123‑1 (MQ=255)
                                                                                                                                  aGCAATATGGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGt  >  1:138004/1‑140 (MQ=255)
                                                                                                                                           |                                                                                                                                      
ACTTCTGTGGATGGCAGGAAGGATCAACCCGCCCAGGCGATGCTGGAACATGTCGAGAGCTTTCTTCTGGAAAGTAACAGCGGGGAAAGCCAGGAACTTCCTCTCTCAGTCACGTTACATCTGCAGACTCAGCAATATGGGGGGTTACAGCGCCGTTTTGCATTACCTGAGCAGTTAGCGAGAGAAGAGTCGCCAGCTCAGACGCAGGCCGGGAATAATAACCATGAATAATGAACAACGCGGCGTCGCGCTGCTCATTGTACTCATGT  >  NC_000913/3462701‑3462969

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: