Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,048,285 A→C 8.4% V136G (GTC→GGC)  gfcC ← putative O‑antigen capsule production periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,048,2850AC 8.4% 144.3 / 4.9 57V136G (GTC→GGC) gfcCputative O‑antigen capsule production periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (32/20);  new base C (3/2);  total (35/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 = 9.66e-01

CGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCC  >  NC_000913/1048145‑1048429
                                                                                                                                            |                                                                                                                                                
cgcCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACgg                                                                                                                                           <  3:68164/148‑1 (MQ=255)
     gCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGAt                                                                                                                                      >  3:395313/1‑148 (MQ=255)
     gCAAGGCGGGGGTGATCTTGCAGATAAACGGTAACGCTACTGCCTGCCTGCCAGTGCAATTCCCCGGCCCCAGAAAACGCACCCAGAAGATTAATACTAAAAGGACGCTTTACGGTATAGAGCGTATAGTCCCCGACCAAcctggtat                                                                                                                                      >  8:154867/1‑141 (MQ=255)
         ggCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCt                                                                                                                                  >  7:309272/1‑148 (MQ=255)
                        gCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATcc                                                                                                                        <  4:395313/143‑1 (MQ=255)
                        gCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATcc                                                                                                                        <  6:45777/143‑1 (MQ=255)
                        gCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCAc                                                                                                                      <  8:309272/145‑1 (MQ=255)
                              aaTCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACgg                                                                                                                                           <  1:30826/118‑1 (MQ=255)
                              aaTCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACgg                                                                                                                                           >  2:30826/1‑118 (MQ=255)
                                   gTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAg                                                                                                        <  1:210193/148‑1 (MQ=255)
                                         cTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTc                                                                                                    >  1:118481/1‑146 (MQ=255)
                                            cGGCCTGCCGGCCCTGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCCGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAg                                                                                                 <  6:123672/146‑1 (MQ=255)
                                             ggCCTGCCTGCCAGGGCAAATGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCGGACCAACGGGGGATTGCTGTTCTCCTCCACTCGTACAAAATCAGGGTCCAGttt                                                                                              <  6:37605/148‑1 (MQ=255)
                                                  gcctgccAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATcc                                                                                                                        >  6:205782/1‑117 (MQ=255)
                                                  gcctgccAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATcc                                                                                                                        <  5:205782/117‑1 (MQ=255)
                                                  gcctgccAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCg                                                                                         >  3:70364/1‑148 (MQ=255)
                                                          gggCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGa                                                                                 <  2:118481/148‑1 (MQ=255)
                                                           ggCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGAc                                                                                >  2:79264/1‑148 (MQ=255)
                                                               aTTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGt                                                                            >  4:400583/1‑148 (MQ=255)
                                                                 tGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGa                                                                          >  8:413770/1‑148 (MQ=255)
                                                                    ccGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGatat                                                                       >  3:98480/1‑148 (MQ=255)
                                                                         cccAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGAcc                                                                               >  6:237388/1‑135 (MQ=255)
                                                                         cccAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGAcc                                                                               <  5:237388/135‑1 (MQ=255)
                                                                            aGAAACCGCACCCAGCAGAGTAATATTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATattgagattg                                                               <  5:124368/148‑1 (MQ=255)
                                                                            aGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACTGTATAGAGCGTATAGTCCCCTACCAACGGGGGAGTGCTGTTGTCATCCACTCCTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATattgagatt                                                                >  6:352105/1‑147 (MQ=255)
                                                                             gAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCCGTGATattgagaat                                                                >  6:124368/1‑144 (MQ=255)
                                                                             gAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATattgagattga                                                              <  3:400583/148‑1 (MQ=255)
                                                                                        cagcagAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCGCCGACCACCGGGGGATTGCTGTTCTCATCCAGTCGTACAAAATCAGGGTCCAGTATGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTgg                                                   >  6:309290/1‑148 (MQ=255)
                                                                                        cagcagAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTgg                                                   >  8:368354/1‑148 (MQ=255)
                                                                                             gagTAATAGTCACCGGACTCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGCTTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCAc                                              >  1:342848/1‑148 (MQ=255)
                                                                                                       cACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATc                                                                                                          <  5:285361/78‑1 (MQ=255)
                                                                                                       cACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATc                                                                                                          >  6:285361/1‑78 (MQ=255)
                                                                                                       cACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATattgaga                                                                  <  6:395166/118‑1 (MQ=255)
                                                                                                       cACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATattgaga                                                                  >  5:395166/1‑118 (MQ=255)
                                                                                                       cACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATc                                    >  1:222135/1‑148 (MQ=255)
                                                                                                             aCGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTgcc                              >  3:210849/1‑148 (MQ=255)
                                                                                                             aCGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTgcc                              >  7:60272/1‑148 (MQ=255)
                                                                                                                cTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACAGATTTAATCGTTgccgcc                           <  4:98480/148‑1 (MQ=255)
                                                                                                                  gTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTgccgccc                          >  2:232928/1‑146 (MQ=255)
                                                                                                                   tACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCcaca                        >  4:385379/1‑148 (MQ=255)
                                                                                                                     cGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACAtc                      >  3:197138/1‑148 (MQ=255)
                                                                                                                      ggTATAGAGCGTATAGTCCCCGCCCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGa                                           >  6:194363/1‑126 (MQ=255)
                                                                                                                      ggTATAGAGCGTATAGTCCCCGCCCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGa                                           >  6:187943/1‑126 (MQ=255)
                                                                                                                      ggTATAGAGCGTATAGTCCCCGCCCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGa                                           <  5:194363/126‑1 (MQ=255)
                                                                                                                      ggTATAGAGCGTATAGTCCCCGCCCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGa                                           <  5:187943/126‑1 (MQ=255)
                                                                                                                             agCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAg              >  6:226223/1‑148 (MQ=255)
                                                                                                                               cGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCt            >  6:346672/1‑148 (MQ=255)
                                                                                                                                   taGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTg        <  8:60272/148‑1 (MQ=255)
                                                                                                                                    aGTCCCCGCCCAACGGGGGATTGCTGTTCTCCTCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGAGATTGAGCTTGAGTAACTGCTGGCGCACGGATGAAATCGTTGCCGCCCCAGCGTCGTCAGCTTCTg        >  2:372419/1‑147 (MQ=255)
                                                                                                                                    aGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTg                                >  3:186290/1‑123 (MQ=255)
                                                                                                                                    aGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTg                                <  4:186290/123‑1 (MQ=255)
                                                                                                                                     gTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCt      >  2:118645/1‑148 (MQ=255)
                                                                                                                                       ccccGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATcgtcgt                 >  3:68722/1‑135 (MQ=255)
                                                                                                                                       ccccGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATcgtcgt                 <  4:68722/135‑1 (MQ=255)
                                                                                                                                         ccGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTccc  <  1:79264/148‑1 (MQ=255)
                                                                                                                                         ccGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGATTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCTCGGAGTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTccc  >  5:399881/1‑148 (MQ=255)
                                                                                                                                          cGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTcc   >  3:148951/1‑146 (MQ=255)
                                                                                                                                            |                                                                                                                                                
CGCCTGCAAGGCGGGGGTGATCTTGCAGATAATCGGTAACGCTACGGCCTGCCTGCCAGGGCAATTGCCCGGCCCCAGAAACCGCACCCAGCAGAGTAATAGTCACCGGACGCTGTACGGTATAGAGCGTATAGTCCCCGACCAACGGGGGATTGCTGTTCTCATCCACTCGTACAAAATCAGGGTCCAGTTTGACCGGTAAGCGACCAGTGATATTGAGATTGAGTAACTGCTGGCGCACGGATTTAATCGTTGCCGCCACATCGTCGTCAGCTTCTGCTTCCC  >  NC_000913/1048145‑1048429

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