Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 561,696 (G)6→7 coding (2000/2604 nt) sfmD → putative outer membrane export usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,6901.G86.2% 195.2 / 28.5 65coding (1994/2604 nt)sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base . (5/4);  new base G (24/32);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑G‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGC  >  NC_000913/561595‑561784
                                                                                                 |                                                                                                
ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑g‑g‑gg                                                                                            >  1:1515329/1‑100 (MQ=255)
ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑g‑g‑gg                                                                                            >  2:1263291/1‑100 (MQ=255)
ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑g‑g‑gg                                                                                            >  1:93053/1‑100 (MQ=255)
ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑g‑g                                                                                               >  1:1346815/1‑98 (MQ=255)
    aaCCTGAGTTATCGCGTAAGCCAGGGTGATACCTCGAAGAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑g‑g‑gggggc                                                                                        >  1:1121873/1‑98 (MQ=255)
    aaCCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCAAGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑g‑g‑gggggc                                                                                        >  1:361625/1‑98 (MQ=255)
      ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑G‑G‑GGGGCggg                                                                                      <  1:113869/100‑1 (MQ=255)
       cTGAGTTATCACGTAAGCCAGGGTGATCCCTCGAATAATGGGTACACGGGCAGCGCCACGGCCCACTGGCAGGCCGCTTACGGTACGC‑t                                                                                                   >  1:1020193/1‑89 (MQ=255)
       cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑g‑g‑gggggcggg                                                                                     <  2:213123/100‑6 (MQ=255)
        tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGt                                                                                     <  1:247476/100‑1 (MQ=255)
        tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGt                                                                                     <  2:1510627/100‑1 (MQ=255)
         gAGTTATCACGTAAGACAGGGTGATACCTCGAACAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑g‑g‑gggggcggg                                                                                     <  2:178814/98‑6 (MQ=255)
           gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGtata                                                                                  >  1:1668079/1‑100 (MQ=255)
            ttATCACGTAAGCCTGGGTTTTACCTCGAATAATGTATACACGGGCAGGGCCTCGGCAAACTGTCAGGCCGCTTTCGGTACGC‑T‑g‑g‑gggggggggtata                                                                                 <  1:413084/99‑6 (MQ=255)
            ttATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑G‑G‑GGGGCGGGTATAAc                                                                                <  2:1142356/100‑1 (MQ=255)
            ttATCACGTAAGCCAGGGTGATACCTCGAATAAGGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATaa                                                                                 <  1:1717696/100‑1 (MQ=255)
             tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAAc                                                                                <  2:1507560/100‑1 (MQ=255)
              aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACt                                                                               <  1:547907/100‑1 (MQ=255)
              aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACt                                                                               <  2:505443/100‑1 (MQ=255)
              aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCAACCGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑G‑G‑GGGGCGGGTATAACTa                                                                              >  1:1095582/1‑100 (MQ=255)
                  cGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGCCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGa                                                                           >  1:1122152/1‑100 (MQ=255)
                     aaGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACcg                                                                        <  1:1761523/100‑1 (MQ=255)
                      aGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACcgc                                                                       >  1:1382076/1‑100 (MQ=255)
                       gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACcgcg                                                                      >  2:983874/1‑100 (MQ=255)
                        ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACcgc                                                                       >  1:696388/1‑98 (MQ=255)
                         cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGAt                                                                    >  1:1350650/1‑100 (MQ=255)
                           gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCa                                                                  >  2:928121/1‑100 (MQ=255)
                           gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCa                                                                  >  1:1042779/1‑100 (MQ=255)
                              tGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACa                                                               <  2:1668090/100‑1 (MQ=255)
                              tGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACa                                                               <  2:1658604/100‑1 (MQ=255)
                                  aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                           <  2:716709/100‑1 (MQ=255)
                                  aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                           <  1:506499/100‑1 (MQ=255)
                                   ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                          >  1:1210974/1‑100 (MQ=255)
                                   ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGAATAACTACGACCGCGATCAACATGACg                                                          >  2:561176/1‑100 (MQ=255)
                                    ctgtaataatGGTTACACGGGCAGCGGCGCGGCAAACTGGCAGGCCGCTTACGGGACGC‑TGG‑G‑GGGGCGGGTATAAGTAGGACCGCGAGCAACATGag                                                           <  1:766804/94‑2 (MQ=255)
                                      cGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACTC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGtt                                                        <  1:1656222/99‑1 (MQ=255)
                                      cGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTa                                                       >  1:1519969/1‑100 (MQ=255)
                                         ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                    <  1:531948/100‑1 (MQ=255)
                                         ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑GAG‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                    >  2:1001120/1‑100 (MQ=255)
                                              gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                               <  2:361627/100‑1 (MQ=255)
                                                 tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                            <  2:1306903/100‑1 (MQ=255)
                                                   cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTCACGGTACGC‑T‑GAG‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGGCAACTGGCAGCTGTc                                          >  1:1804167/1‑100 (MQ=255)
                                                     ngggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                        >  1:1411316/2‑100 (MQ=255)
                                                      gggCAGCGCCACGGCAAACTGGCAGGACGCTTACGGTACGCAG‑G‑G‑GGGGCGGGTATAACTACGACCGCGATAAACATGACGTTAACTGGCAGCTGTCTgg                                       >  1:1029825/1‑100 (MQ=255)
                                                       ggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGc                                      >  1:325130/1‑100 (MQ=255)
                                                            cgcCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtgt                                 <  2:443495/100‑1 (MQ=255)
                                                            cgcCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtgt                                 <  1:983866/100‑1 (MQ=255)
                                                              ccACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTgg                               <  2:1121884/100‑1 (MQ=255)
                                                              ccACGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTgg                               <  2:1270485/100‑1 (MQ=255)
                                                                aCGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                             >  1:1511136/1‑100 (MQ=255)
                                                                 cGGCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCg                            <  1:483625/100‑1 (MQ=255)
                                                                   gCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                          <  1:1115577/100‑1 (MQ=255)
                                                                   gCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCgg                           >  1:1414765/1‑99 (MQ=255)
                                                                   gCAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTAGTCgg                           >  2:389375/1‑99 (MQ=255)
                                                                    cAAACTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGc                         <  2:250875/100‑1 (MQ=255)
                                                                     aaaCTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCa                        >  1:348031/1‑100 (MQ=255)
                                                                       aCTGGCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTTTCTGGCGGTGTGGTCGGGCa                        <  1:844550/98‑1 (MQ=255)
                                                                          ggCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaa                   >  2:212720/1‑100 (MQ=255)
                                                                          ggCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaa                   >  2:428191/1‑100 (MQ=255)
                                                                          ggCAGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaa                    <  1:1655771/99‑1 (MQ=255)
                                                                           tgaGGCCGCTTATGGTACGCTT‑G‑T‑TAGAAGGGTATGTCGAGTACCGCGATCAACATGAGGTTAACTGGCAGCTGTCTGGCGTTGTGGTCGGGCATGaaaa                  <  1:505523/98‑1 (MQ=255)
                                                                            cAGGCCGCTTACAGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAAt                 >  2:1194285/1‑100 (MQ=255)
                                                                             aGGCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATg                >  2:813714/1‑100 (MQ=255)
                                                                              ggCCGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg               <  1:235027/100‑1 (MQ=255)
                                                                                ccGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa             <  2:1788971/100‑1 (MQ=255)
                                                                                 cGCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCAt            <  2:486181/100‑1 (MQ=255)
                                                                                  gCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATa           <  2:118640/100‑1 (MQ=255)
                                                                                  gCTTACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATa           <  2:865656/100‑1 (MQ=255)
                                                                                  gCTTACGGTACGC‑T‑G‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATaa          >  2:1149620/1‑100 (MQ=255)
                                                                                     tACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGTCAGCTGTCTGGCGGTGTTGTCGGGCATGAAAATGGCATAACg        >  1:603249/1‑100 (MQ=255)
                                                                                     tACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGTTGTCTGGCGGTGTGGTCGGGCATGGAAATGGCATAACg        <  1:1120117/100‑1 (MQ=255)
                                                                                     tACGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAAATGGCAGCTGACTAACGGTGTGGTCAAGCATGAAAATGGCATAACg        >  2:553696/1‑100 (MQ=255)
                                                                                       cGGTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGc       <  1:942641/99‑1 (MQ=255)
                                                                                         gTACGC‑TGG‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGa    >  2:245502/1‑100 (MQ=255)
                                                                                           aCGC‑T‑G‑GAGGGGCAGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTTTCCGGCGGTGTGGTCGGTCATGAAAATGGCATAACGCTGagc  <  1:1462320/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGC‑T‑G‑G‑GGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGC  >  NC_000913/561595‑561784

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