Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,316,160 G→A *891* (TAG→TAA rcsD → phosphotransfer intermediate protein in two‑component regulatory system with RcsBC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,316,1600GA98.4% 203.6 / NA 64*891* (TAG→TAArcsDphosphotransfer intermediate protein in two‑component regulatory system with RcsBC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (29/34);  minor base C (0/1);  total (29/35)
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

GTATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTA  >  NC_000913/2316026‑2316291
                                                                                                                                      |                                                                                                                                   
gTATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAg                                                                                                                                 >  1:528598/1‑139 (MQ=255)
       ccATGCTAAATCTGGTACCCGTCAAGCAGTTATGTAAAACACTGGAACATCTTATTCGTGAGAAGGATGTTCTAGGAATAGAAGAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACACGGTAGCCt                                                                                                                          <  2:649384/139‑1 (MQ=255)
        cATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTa                                                                                                                         <  1:667084/139‑1 (MQ=255)
            cTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTAc                                                                                                                     <  2:609925/139‑1 (MQ=255)
             tAAATCTGGTACCCGGCAAGGCGTTATGTGTAACGCTGGAACATCTGATTCGTGAAATGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTACCAAGGTAGCCTATTACa                                                                                                                    <  2:66805/139‑1 (MQ=255)
               aaTCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATg                                                                                                                  >  2:99763/1‑139 (MQ=255)
                 tCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGaa                                                                                                                <  2:566595/139‑1 (MQ=255)
                     gTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAat                                                                                                            >  2:417587/1‑139 (MQ=255)
                     gTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAat                                                                                                            >  1:536623/1‑139 (MQ=255)
                       aCCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAatat                                                                                                          <  1:367195/139‑1 (MQ=255)
                       aCCCGGCAACCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAatat                                                                                                          <  2:528598/139‑1 (MQ=255)
                           ggCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAAc                                                                                                      <  1:355379/139‑1 (MQ=255)
                               aGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                >  2:433505/1‑109 (MQ=255)
                               aGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                <  1:433505/109‑1 (MQ=255)
                                      aTGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGc                                                                                           <  1:326384/139‑1 (MQ=255)
                                          gAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGAt                                                                                       <  1:597542/139‑1 (MQ=255)
                                           aaaCGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAg                                                                                                                                 >  1:392551/1‑96 (MQ=255)
                                           aaaCGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAg                                                                                                                                 <  2:392551/96‑1 (MQ=255)
                                             aCGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGAc                                                                                    <  1:432366/139‑1 (MQ=255)
                                              cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                          >  2:645577/1‑100 (MQ=255)
                                              cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCt                                                                                                                          <  1:645577/100‑1 (MQ=255)
                                              cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAAc                                                                                                      >  2:97311/1‑120 (MQ=255)
                                              cGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAAc                                                                                                      <  1:97311/120‑1 (MQ=255)
                                                   gAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCg                                                                              >  2:252232/1‑139 (MQ=255)
                                                    aaCATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGa                                                                             <  1:417587/139‑1 (MQ=255)
                                                         cTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTc                                                                        >  1:237280/1‑139 (MQ=255)
                                                           gATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCCAGAGCTTGCTGTAACAAGGTAGCCTATTCCATGAACCATATGCACGTAATTATTGCCGATGACCATCCGATAGTCtt                                                                      >  2:324639/1‑139 (MQ=255)
                                                                 tGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCgg                                                                >  2:426458/1‑139 (MQ=255)
                                                                      aGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTc                                                           >  1:598947/1‑139 (MQ=255)
                                                                           gTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCaa                                                       >  1:599337/1‑138 (MQ=255)
                                                                           gTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCaa                                                       <  2:599337/138‑1 (MQ=255)
                                                                            ttCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACa                                                                                                                    >  1:244318/1‑76 (MQ=255)
                                                                            ttCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACa                                                                                                                    <  2:244318/76‑1 (MQ=255)
                                                                            ttCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAAt                                                     >  2:478839/1‑139 (MQ=255)
                                                                            ttCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAAt                                                     >  2:218579/1‑139 (MQ=255)
                                                                             tCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                >  2:373139/1‑63 (MQ=255)
                                                                             tCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAgg                                                                                                                                <  1:373139/63‑1 (MQ=255)
                                                                                   aaTAGTAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATTAACAATATGAACGTAATTATTGCCGATGACCATCCGATATTCTTGTTCGGTATTCGCAAATCACTTGa                                              <  1:478839/139‑1 (MQ=255)
                                                                                   aaTAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGa                                              <  1:218579/139‑1 (MQ=255)
                                                                                     tAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGc                                            >  2:31970/1‑139 (MQ=255)
                                                                                      aGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGc                                                                                                                            >  2:251730/1‑58 (MQ=255)
                                                                                      aGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGc                                                                                                                            <  1:251730/58‑1 (MQ=255)
                                                                                       gAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCaa                                          >  1:16982/1‑139 (MQ=255)
                                                                                        aaaaaTACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCaaa                                         <  2:598947/139‑1 (MQ=255)
                                                                                          aaaTACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAAtt                                       <  2:536623/139‑1 (MQ=255)
                                                                                             tACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAg                                    <  1:129602/139‑1 (MQ=255)
                                                                                                  cAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTAc                                                                                                                     <  1:85877/53‑1 (MQ=255)
                                                                                                  cAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTAc                                                                                                                     >  2:85877/1‑53 (MQ=255)
                                                                                                  cAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGa                                                                                     <  1:443543/85‑1 (MQ=255)
                                                                                                  cAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGa                                                                                     >  2:443543/1‑85 (MQ=255)
                                                                                                  cAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGt                               >  2:322477/1‑139 (MQ=255)
                                                                                                       aCATTGACAGTTATGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAAtg                          >  1:362508/1‑139 (MQ=255)
                                                                                                                    tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAg                                    <  2:361454/116‑1 (MQ=255)
                                                                                                                    tGTCAAGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAg                                    >  1:361454/1‑116 (MQ=255)
                                                                                                                        aaGAGCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGaa         <  2:231770/139‑1 (MQ=255)
                                                                                                                         agagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAg        >  1:469296/1‑139 (MQ=255)
                                                                                                                          gagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa       <  2:723082/139‑1 (MQ=255)
                                                                                                                          gagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa       <  2:16982/139‑1 (MQ=255)
                                                                                                                          gagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa       <  1:324639/139‑1 (MQ=255)
                                                                                                                          gagCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCACAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa       <  2:388533/139‑1 (MQ=255)
                                                                                                                          gagCTTGCTGTAACAAGGTAGCATATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTGGTTCGGTATTCGCAAATCAATTGATCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGa       <  1:525181/139‑1 (MQ=255)
                                                                                                                           agCTTGCTGTAACAAGGTAGCCTATTACATGAACAATATGAAGGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGAc      >  1:503948/1‑139 (MQ=255)
                                                                                                                               tGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTa  <  1:245640/139‑1 (MQ=255)
                                                                                                                               tGCTGTAACAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTa  >  2:311094/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                   
GTATTTGCCATGCTAAATCTGGTACCCGGCAAGCAGTTATGTGAAACGCTGGAACATCTGATTCGTGAGAAGGATGTTCCAGGAATAGAAAAATACATCAGCGACATTGACAGTTATGTCAAGAGCTTGCTGTAGCAAGGTAGCCTATTACATGAACAATATGAACGTAATTATTGCCGATGACCATCCGATAGTCTTGTTCGGTATTCGCAAATCACTTGAGCAAATTGAGTGGGTGAATGTTGTCGGCGAATTTGAAGACTCTA  >  NC_000913/2316026‑2316291

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