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.G98.6% 248.7 / ‑3.3 71coding (1994/2604 nt)sfmDputative outer membrane export usher protein
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (35/35);  total (35/36)
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.44e-01

GCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTG  >  NC_000913/561594‑561781
                                                                                                 |                                                                                           
gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggg                                                                                          >  2:1226156/1‑100 (MQ=255)
gCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTggg                                                                                          >  2:2815393/1‑100 (MQ=255)
 ccACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggg                                                                                         <  1:715806/100‑1 (MQ=255)
      aCCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcg                                                                                    >  2:1566947/1‑97 (MQ=255)
       ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcgg                                                                                   >  1:2126079/1‑96 (MQ=255)
       ccTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcg                                                                                    <  2:1367090/99‑4 (MQ=255)
        cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcggg                                                                                  <  2:806867/100‑6 (MQ=255)
        cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcggg                                                                                  <  1:2354081/100‑6 (MQ=255)
        cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcggg                                                                                  >  1:1981739/1‑95 (MQ=255)
        cTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTgggggggcggg                                                                                  <  1:1207062/100‑6 (MQ=255)
         tGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGt                                                                                  <  1:1571461/100‑1 (MQ=255)
           aGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtat                                                                                <  2:2807766/100‑1 (MQ=255)
           aGTGATCACGTAAGCCAGGGTGATACCTCGAACAATGGGTACACGGGCAGCGCCACGGCACACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGtata                                                                               <  1:2748061/100‑1 (MQ=255)
            gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtata                                                                               >  2:1984359/1‑100 (MQ=255)
            gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtata                                                                               <  2:2865961/100‑1 (MQ=255)
            gTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGtata                                                                               <  2:2883850/100‑1 (MQ=255)
              tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAAc                                                                             <  2:735515/100‑1 (MQ=255)
              tATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAAc                                                                             <  2:1037176/100‑1 (MQ=255)
               aTCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACt                                                                            <  2:1981725/100‑1 (MQ=255)
                 cACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTAc                                                                          >  2:1936111/1‑100 (MQ=255)
                 cACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTAc                                                                          >  1:195018/1‑100 (MQ=255)
                 cACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTAc                                                                          >  1:956180/1‑100 (MQ=255)
                     tAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGAcc                                                                      <  2:2602414/100‑1 (MQ=255)
                     tAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGAcc                                                                      >  2:1021289/1‑100 (MQ=255)
                     tAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGAACGCTGGGGGGGCGGGTAGAACTACGAcc                                                                      >  1:1406344/1‑100 (MQ=255)
                      aaGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACcg                                                                     <  1:2357251/100‑1 (MQ=255)
                          cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGAt                                                                 <  2:940019/100‑1 (MQ=255)
                          cAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGACAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGAt                                                                 >  2:734864/1‑100 (MQ=255)
                            gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCa                                                               >  1:2136630/1‑100 (MQ=255)
                            gggTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCa                                                               <  2:2872026/100‑1 (MQ=255)
                              gTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAAc                                                             >  1:2886912/1‑100 (MQ=255)
                                gATACCTCGAATAATGGGTACACGGGCAGCGCCACGGGAAGCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATCACTACGACCGCGATCAACa                                                            <  2:780163/99‑1 (MQ=255)
                                    ccTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACg                                                       >  1:1810968/1‑100 (MQ=255)
                                     cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                      >  2:1943622/1‑100 (MQ=255)
                                     cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                      >  1:1294475/1‑100 (MQ=255)
                                     cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                      <  2:550909/100‑1 (MQ=255)
                                     cTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGt                                                      <  2:1618745/100‑1 (MQ=255)
                                        gAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTaa                                                   >  2:822024/1‑100 (MQ=255)
                                          ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                 >  1:43730/1‑100 (MQ=255)
                                          ataatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACt                                                 >  1:347119/1‑100 (MQ=255)
                                           taatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTg                                                <  1:2070790/100‑1 (MQ=255)
                                            aatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTgg                                               <  1:2138357/100‑1 (MQ=255)
                                            aatGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTgg                                               >  2:315597/1‑100 (MQ=255)
                                             atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGc                                              >  1:875671/1‑100 (MQ=255)
                                             atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGc                                              <  2:3348/100‑1 (MQ=255)
                                             atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGc                                              <  2:2809766/100‑1 (MQ=255)
                                             atGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGc                                              <  1:1733679/100‑1 (MQ=255)
                                              tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                             <  2:2868871/100‑1 (MQ=255)
                                              tGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCa                                             <  1:1726856/100‑1 (MQ=255)
                                               ggttACACGGGCAGCGCCACGGCAAGCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAg                                            <  2:1721267/97‑1 (MQ=255)
                                                ggTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGc                                           >  1:691079/1‑100 (MQ=255)
                                                    cacGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTc                                       >  1:376888/1‑100 (MQ=255)
                                                     acGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCt                                      <  1:1049631/100‑1 (MQ=255)
                                                     acGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCt                                      <  2:2064079/100‑1 (MQ=255)
                                                      cGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTg                                     >  1:1755567/1‑100 (MQ=255)
                                                        ggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGc                                   >  1:567081/1‑100 (MQ=255)
                                                        ggCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGc                                   <  2:711235/100‑1 (MQ=255)
                                                         gCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCg                                  >  1:1074108/1‑100 (MQ=255)
                                                           aGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgt                                >  1:786457/1‑100 (MQ=255)
                                                            gcgcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGgtg                               >  2:2834741/1‑100 (MQ=255)
                                                              gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             >  2:2199090/1‑100 (MQ=255)
                                                              gcCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             >  1:1422687/1‑100 (MQ=255)
                                                              gcCACGGCAAACTGGCAGGCCGATTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGtgtg                             >  2:1677222/1‑100 (MQ=255)
                                                                cACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGt                           <  2:379033/100‑1 (MQ=255)
                                                                   ggCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGCATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCgg                        >  2:1584704/1‑100 (MQ=255)
                                                                      aaaCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCa                     <  1:2822757/100‑1 (MQ=255)
                                                                       aaCTGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCAt                    <  1:734866/100‑1 (MQ=255)
                                                                          ttccgatctgCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATAAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaa                 <  1:423325/91‑1 (MQ=255)
                                                                          tGGCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaa                 <  1:822025/100‑1 (MQ=255)
                                                                            gCAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGaaaa               <  2:2057948/100‑1 (MQ=255)
                                                                             cAGGCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGCCGTTAACTGGCAGCTGTCTGGCGGCGTGGTCGGGCATGAAAAt              >  2:3923/1‑100 (MQ=255)
                                                                               ggCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATg             <  2:1787129/99‑1 (MQ=255)
                                                                                gCCGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGACGGTGTGGTCGGGCATGAAACTgg            >  1:844817/1‑99 (MQ=255)
                                                                                 ccGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCa          >  1:1471108/1‑100 (MQ=255)
                                                                                  cGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCAt         <  1:1226160/100‑1 (MQ=255)
                                                                                  cGCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCAt         <  2:2136617/100‑1 (MQ=255)
                                                                                   gCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATa        >  1:556367/1‑100 (MQ=255)
                                                                                   gCTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCGGTCTGGCGGTGTGGTCGGGCATGaaagatcggaa        >  2:423324/1‑92 (MQ=255)
                                                                                    cTTACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATaa       >  1:1238141/1‑100 (MQ=255)
                                                                                     ttACGGTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAAc      <  2:789449/100‑1 (MQ=255)
                                                                                         ggTACGCTGGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTg  >  2:247272/1‑100 (MQ=255)
                                                                                                 |                                                                                           
GCCACAACCTGAGTTATCACGTAAGCCAGGGTGATACCTCGAATAATGGGTACACGGGCAGCGCCACGGCAAACTGGCAGGCCGCTTACGGTACGCTGGGGGGCGGGTATAACTACGACCGCGATCAACATGACGTTAACTGGCAGCTGTCTGGCGGTGTGGTCGGGCATGAAAATGGCATAACGCTG  >  NC_000913/561594‑561781

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