Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 750,624 C→T 7.5% A516T (GCA→ACA)  ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913750,6240CT 7.5% 118.9 / 2.5 53A516T (GCA→ACA) ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (29/20);  new base T (2/2);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.54e-01

ATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGT  >  NC_000913/750478‑750764
                                                                                                                                                  |                                                                                                                                            
aTAGCCAACCCGCCCCCCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGATGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAAAGACGGCGAAATCGAGATATCTCTCCAGTCAACAAAATAACAATTAAAACctgctg                                                                                                                                            >  2:655500/1‑149 (MQ=255)
    ccAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGAt                                                                                                                                        >  2:1010/1‑149 (MQ=255)
           gACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCg                                                                                                                                 >  1:819988/1‑149 (MQ=255)
           gACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCg                                                                                                                                 >  2:366310/1‑149 (MQ=255)
            cccaccGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCgg                                                                                                                                <  2:242307/148‑1 (MQ=255)
                ccGTTACCGAAGGGGAGCGAAAGCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGCCGGCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGt                                                                                                                            <  1:178802/149‑1 (MQ=255)
                 cGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTg                                                                                                                           >  2:626650/1‑149 (MQ=255)
                     aCCGCAGGGGAGCTAAATCGACAGATAATGCTTGTTGTCGCGAGCTTTATCTTCGTAGTGGGTCGTTTTGAATGACGTCGAAATCGATATATCGCTCCAGTCTCCAATATCAACATTAACACCTGCTGTGATGTTGGCGGTCGTGGAgg                                                                                                                       <  2:768234/149‑1 (MQ=255)
                       cGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGTGATATCTATCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGTCGGTCGTGGTGGCg                                                                                                                     >  2:257963/1‑149 (MQ=255)
                       cGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCCCCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCg                                                                                                                     >  2:382453/1‑149 (MQ=255)
                        gAAGGGGAGCGAAATCGACAGAAAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCACCAGTCCCCAATATCAACCTTAAAACCTGCTGTGATGTTGGCGGTCGTGGCGGCGt                                                                                                                    >  2:932174/1‑149 (MQ=255)
                                 cGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTcc                                                                                                           >  2:494901/1‑149 (MQ=255)
                                 cGAAATCGACAGATAAATCTCGTTGTCGCGATCATTATCTTCGTAGTGGGTCGTATTGAATGACCTCGAAATCGAGATATCTCACCAGTCACCAATATAAACATTAACACCTGCTTTGAAGAATGCGTTCGTGGATGCGTCTCCATTcc                                                                                                           >  2:665364/1‑149 (MQ=255)
                                  gAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGTTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTcca                                                                                                          >  1:1018067/1‑149 (MQ=255)
                                  gAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGTTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTcca                                                                                                          >  1:518573/1‑149 (MQ=255)
                                          aaGATAAAGCTGGTTGTAGCGATCTGTATCTTCGTAGGGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTTCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGttt                                                                                                  <  2:562573/148‑1 (MQ=255)
                                            gATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATTTCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGa                                                                                                <  2:185009/149‑1 (MQ=255)
                                             aTAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGGCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGAt                                                                                               <  1:298371/149‑1 (MQ=255)
                                             aTAAATCTGGATGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGAt                                                                                               <  2:413312/149‑1 (MQ=255)
                                               aaaTCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTTCATTCCACCAGGTTTGATGt                                                                                             >  1:992044/1‑149 (MQ=255)
                                               aaaTCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGt                                                                                             >  1:13041/1‑149 (MQ=255)
                                                 aTCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAg                                                                                           >  2:799920/1‑149 (MQ=255)
                                                    tGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAg                                                                                        <  1:754988/149‑1 (MQ=255)
                                                     ggTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGTTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAgg                                                                                       <  2:1018067/149‑1 (MQ=255)
                                                     ggTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGTTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAgg                                                                                       <  2:518573/149‑1 (MQ=255)
                                                     ggTTGTAGCTATCTTTATCTTCGTTGTGGGTCGGATTGAATGACGTCGAAATCGAGATCTCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGTCGTCTGCATTCCACCATGTTTGAGGTAGCAgg                                                                                       <  2:188212/149‑1 (MQ=255)
                                                        tctCGCGATCTTTATCTTCGTTGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTg                                                                                    <  1:958186/147‑1 (MQ=255)
                                                          tCGCGATCTTTATCTTCGTAGTGGGTCGTTTTGAATGACGGCGATATCGAGATATCTCTCCAGTCACCTATCTCAACTTTTAACCCTGCTGTGATGCTGGCGGTCGGGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGc                                                                                  <  2:993451/149‑1 (MQ=255)
                                                          tCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGc                                                                                  >  2:623830/1‑149 (MQ=255)
                                                            gcgATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGTTCCTGGAGCCGTCTGCATTAAACCAGGTTTGAAGTAGCAGGTTggaga                                                                                >  2:958186/1‑146 (MQ=255)
                                                                 cTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGGTTGATGTAGCAGGTTGGCGTAAAGa                                                                           >  2:526887/1‑149 (MQ=255)
                                                                        ttCGTTGTTGGTCGTATTGAATGACGTCGAATTCGAGATATCTCTCCAGTAACCTATATCATCATTACCAGCTGCTGTGATGTTGGCTGTCGTGGAGGCGGCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGt                                                                    <  1:655500/149‑1 (MQ=255)
                                                                         tCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGtt                                                                   <  2:379961/149‑1 (MQ=255)
                                                                         tCGTAGGGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGGCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGtt                                                                   <  1:228494/149‑1 (MQ=255)
                                                                            tAGTGGGGCGTATTGACTGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAACCCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAg                                                                <  2:534715/149‑1 (MQ=255)
                                                                                gggTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCAGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACAAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTggg                                                            <  1:435331/149‑1 (MQ=255)
                                                                                             tGACGTCGAAATCGAGATATCTCTCCAGTCCCCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTCATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGc                                               >  1:498221/1‑149 (MQ=255)
                                                                                                cGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGATGGccga                                            >  2:832742/1‑147 (MQ=255)
                                                                                                 gTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGTATGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCAc                                           >  1:695538/1‑149 (MQ=255)
                                                                                                 gTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATCTACCAGGTTGGCGTAAAGATTGAGGTTTAGTGCGCTAATCGGTTGGCCCAc                                           >  2:56098/1‑149 (MQ=255)
                                                                                                    gAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGa                                        >  2:809479/1‑149 (MQ=255)
                                                                                                    gAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCAGCATTCCACCAGGTTTGATGTAGCAGGTTTGCGTTAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGa                                        >  1:963052/1‑149 (MQ=255)
                                                                                                     aaaTCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGGGGGGTATTCGGTTGGCCCACGGAt                                       >  1:659376/1‑149 (MQ=255)
                                                                                                         cGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCCGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGc                                   >  2:214151/1‑149 (MQ=255)
                                                                                                                tctctcCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCCCCAGGTTTGATGTAGCAGGTGGGCGTAAAGATTGAGGGTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCgg                            >  1:364495/1‑148 (MQ=255)
                                                                                                                 ctctcCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTc                           >  1:194183/1‑149 (MQ=255)
                                                                                                                     ccAGTCACCAATCTCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGGTGGCGTAAAGATTGAGGTTTTGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGtt                       <  2:819988/149‑1 (MQ=255)
                                                                                                                         tCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGTTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTc                   >  2:266333/1‑149 (MQ=255)
                                                                                                                              aaTATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCt              <  1:799920/149‑1 (MQ=255)
                                                                                                                                aatCAACATTAAAACCTGCTGTGATGGTGGAGGTCGTGGAGGCGTCTGCATTCCACCAGGATAGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTgc            <  2:104979/148‑1 (MQ=255)
                                                                                                                                   cAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCa         <  2:194183/149‑1 (MQ=255)
                                                                                                                                     aCATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCAGATCGTCGGTTTTTCGTCCTGCGCATc       >  2:104457/1‑149 (MQ=255)
                                                                                                                                        ttAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCt    >  1:847549/1‑149 (MQ=255)
                                                                                                                                          aaaaCCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATTTAGCAGGTTGGCGTAAAGATTGAGGTTTATAGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGt  >  2:588143/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
ATAGCCAACCCGACCACCGTTACCGAAGGGGAGCGAAATCGACAGATAAATCTGGTTGTCGCGATCTTTATCTTCGTAGTGGGTCGTATTGAATGACGTCGAAATCGAGATATCTCTCCAGTCACCAATATCAACATTAAAACCTGCTGTGATGTTGGCGGTCGTGGAGGCGTCTGCATTCCACCAGGTTTGATGTAGCAGGTTGGCGTAAAGATTGAGGTTTAGTGGGGTAATCGGTTGGCCCACGGATAAGCTGATCGTCTGTTTTTCGTCCTGCGCATCGCTGT  >  NC_000913/750478‑750764

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