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

CCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGC  >  NC_000913/561595‑561784
                                                                                                 |                                                                                                
ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg                                                                                            >  1:1515329/1‑100 (MQ=255)
ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg                                                                                            >  2:1263291/1‑100 (MQ=255)
ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg                                                                                            >  1:93053/1‑100 (MQ=255)
ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgg                                                                                               >  1:1346815/1‑98 (MQ=255)
    aaCCTGAGTTATCGCGTAAGCCAGGGTGATACCTCGAAGAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggc                                                                                        >  1:1121873/1‑98 (MQ=255)
    aaCCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCAAGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggc                                                                                        >  1:361625/1‑98 (MQ=255)
      ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCggg                                                                                      <  1:113869/100‑1 (MQ=255)
       cTGAGTTATCACGTAAGCCAGGGTGATCCCTCGAATAATGGGTACACGGGCAGCGCCACGGCCCACTGGCAGGCCGCTTACGGTACGCt                                                                                                   >  1:1020193/1‑89 (MQ=255)
       cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcggg                                                                                     <  2:213123/100‑6 (MQ=255)
        tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGt                                                                                     <  1:247476/100‑1 (MQ=255)
        tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGt                                                                                     <  2:1510627/100‑1 (MQ=255)
         gAGTTATCACGTAAGACAGGGTGATACCTCGAACAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcggg                                                                                     <  2:178814/98‑6 (MQ=255)
           gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtata                                                                                  >  1:1668079/1‑100 (MQ=255)
            ttATCACGTAAGCCTGGGTTTTACCTCGAATAATGTATACACGGGCAGGGCCTCGGCAAACTGTCAGGCCGCTTTCGGTACGCTgggggggggggtata                                                                                 <  1:413084/99‑6 (MQ=255)
            ttATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAAc                                                                                <  2:1142356/100‑1 (MQ=255)
            ttATCACGTAAGCCAGGGTGATACCTCGAATAAGGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATaa                                                                                 <  1:1717696/100‑1 (MQ=255)
             tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAAc                                                                                <  2:1507560/100‑1 (MQ=255)
              aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACt                                                                               <  1:547907/100‑1 (MQ=255)
              aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACt                                                                               <  2:505443/100‑1 (MQ=255)
              aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCAACCGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTa                                                                              >  1:1095582/1‑100 (MQ=255)
                  cGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGCCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGa                                                                           >  1:1122152/1‑100 (MQ=255)
                     aaGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcg                                                                        <  1:1761523/100‑1 (MQ=255)
                      aGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcgc                                                                       >  1:1382076/1‑100 (MQ=255)
                       gCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcgcg                                                                      >  2:983874/1‑100 (MQ=255)
                        ccAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcgc                                                                       >  1:696388/1‑98 (MQ=255)
                         cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGAt                                                                    >  1:1350650/1‑100 (MQ=255)
                           gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCa                                                                  >  2:928121/1‑100 (MQ=255)
                           gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCa                                                                  >  1:1042779/1‑100 (MQ=255)
                              tGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACa                                                               <  2:1668090/100‑1 (MQ=255)
                              tGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACa                                                               <  2:1658604/100‑1 (MQ=255)
                                  aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                           <  2:716709/100‑1 (MQ=255)
                                  aCCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGAc                                                           <  1:506499/100‑1 (MQ=255)
                                   ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                          >  1:1210974/1‑100 (MQ=255)
                                   ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGAATAACTACGACCGCGATCAACATGACg                                                          >  2:561176/1‑100 (MQ=255)
                                    ctgtaataatGGTTACACGGGCAGCGGCGCGGCAAACTGGCAGGCCGCTTACGGGACGCTGGGGGGGCGGGTATAAGTAGGACCGCGAGCAACATGag                                                           <  1:766804/94‑2 (MQ=255)
                                      cGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACTCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGtt                                                        <  1:1656222/99‑1 (MQ=255)
                                      cGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTa                                                       >  1:1519969/1‑100 (MQ=255)
                                         ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                    <  1:531948/100‑1 (MQ=255)
                                         ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGAGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                    >  2:1001120/1‑100 (MQ=255)
                                              gggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                               <  2:361627/100‑1 (MQ=255)
                                                 tACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTg                                            <  2:1306903/100‑1 (MQ=255)
                                                   cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTCACGGTACGCTGAGGGGGCGGGTATAACTACGACCGCGATCAACATGACGGCAACTGGCAGCTGTc                                          >  1:1804167/1‑100 (MQ=255)
                                                     ngggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                        >  1:1411316/2‑100 (MQ=255)
                                                      gggCAGCGCCACGGCAAACTGGCAGGACGCTTACGGTACGCAGGGGGGGCGGGTATAACTACGACCGCGATAAACATGACGTTAACTGGCAGCTGTCTgg                                       >  1:1029825/1‑100 (MQ=255)
                                                       ggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGc                                      >  1:325130/1‑100 (MQ=255)
                                                            cgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtgt                                 <  2:443495/100‑1 (MQ=255)
                                                            cgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtgt                                 <  1:983866/100‑1 (MQ=255)
                                                              ccACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTgg                               <  2:1121884/100‑1 (MQ=255)
                                                              ccACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTgg                               <  2:1270485/100‑1 (MQ=255)
                                                                aCGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTc                             >  1:1511136/1‑100 (MQ=255)
                                                                 cGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCg                            <  1:483625/100‑1 (MQ=255)
                                                                   gCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCggg                          <  1:1115577/100‑1 (MQ=255)
                                                                   gCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCgg                           >  1:1414765/1‑99 (MQ=255)
                                                                   gCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTAGTCgg                           >  2:389375/1‑99 (MQ=255)
                                                                    cAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGc                         <  2:250875/100‑1 (MQ=255)
                                                                     aaaCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCa                        >  1:348031/1‑100 (MQ=255)
                                                                       aCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTTTCTGGCGGTGTGGTCGGGCa                        <  1:844550/98‑1 (MQ=255)
                                                                          ggCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaa                   >  2:212720/1‑100 (MQ=255)
                                                                          ggCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaa                   >  2:428191/1‑100 (MQ=255)
                                                                          ggCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaa                    <  1:1655771/99‑1 (MQ=255)
                                                                           tgaGGCCGCTTATGGTACGCTTGTTAGAAGGGTATGTCGAGTACCGCGATCAACATGAGGTTAACTGGCAGCTGTCTGGCGTTGTGGTCGGGCATGaaaa                  <  1:505523/98‑1 (MQ=255)
                                                                            cAGGCCGCTTACAGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAAt                 >  2:1194285/1‑100 (MQ=255)
                                                                             aGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATg                >  2:813714/1‑100 (MQ=255)
                                                                              ggCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATgg               <  1:235027/100‑1 (MQ=255)
                                                                                ccGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa             <  2:1788971/100‑1 (MQ=255)
                                                                                 cGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCAt            <  2:486181/100‑1 (MQ=255)
                                                                                  gCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATa           <  2:118640/100‑1 (MQ=255)
                                                                                  gCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATa           <  2:865656/100‑1 (MQ=255)
                                                                                  gCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATaa          >  2:1149620/1‑100 (MQ=255)
                                                                                     tACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGTCAGCTGTCTGGCGGTGTTGTCGGGCATGAAAATGGCATAACg        >  1:603249/1‑100 (MQ=255)
                                                                                     tACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGTTGTCTGGCGGTGTGGTCGGGCATGGAAATGGCATAACg        <  1:1120117/100‑1 (MQ=255)
                                                                                     tACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAAATGGCAGCTGACTAACGGTGTGGTCAAGCATGAAAATGGCATAACg        >  2:553696/1‑100 (MQ=255)
                                                                                       cGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGc       <  1:942641/99‑1 (MQ=255)
                                                                                         gTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGa    >  2:245502/1‑100 (MQ=255)
                                                                                           aCGCTGGAGGGGCAGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTTTCCGGCGGTGTGGTCGGTCATGAAAATGGCATAACGCTGagc  <  1:1462320/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTGAGC  >  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: