Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,118,164 G→A R81C (CGT→TGT)  wcaK ← colanic acid biosynthesis protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,118,1640GA97.6% 126.3 / NA 41R81C (CGT→TGT) wcaKcolanic acid biosynthesis protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (21/19);  minor base T (0/1);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACAT  >  NC_000913/2118029‑2118295
                                                                                                                                       |                                                                                                                                   
cgccgACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagag                                                                                                                                  <  1:945528/139‑1 (MQ=255)
 gccgACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagagg                                                                                                                                 <  2:1155113/139‑1 (MQ=255)
          gATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTa                                                                                                                        <  2:147609/139‑1 (MQ=255)
             aaTGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAAcg                                                                                                                     <  2:1156753/139‑1 (MQ=255)
                 gCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGcc                                                                                                                 >  2:415848/1‑139 (MQ=255)
                              aCAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagagga                                                                                                                                <  1:1046596/111‑1 (MQ=255)
                              aCAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagagga                                                                                                                                >  2:1046596/1‑111 (MQ=255)
                                  cagGCGCACGAAGTCGGTGAAGCCCTGGGCGATGGCGATATTGCGCAGCTGGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGAGACAGAGGACTTTTGTCACGCGCCCGACAACGCCCGCcgcg                                                                                                <  2:161613/139‑1 (MQ=255)
                                     gcgcACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtg                                                                                             >  1:317113/1‑139 (MQ=255)
                                        cACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgttg                                                                                          >  2:419998/1‑139 (MQ=255)
                                                         ccTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAg                                                                         >  1:414380/1‑139 (MQ=255)
                                                                   ggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCa                                                                                >  1:823212/1‑122 (MQ=255)
                                                                   ggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCa                                                                                <  2:550288/122‑1 (MQ=255)
                                                                   ggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCa                                                                                >  1:550288/1‑122 (MQ=255)
                                                                   ggcgatATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCa                                                                                <  2:823212/122‑1 (MQ=255)
                                                                             cgcAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTa                                                     >  2:663142/1‑139 (MQ=255)
                                                                              gcAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagag                                                                                                                                  >  1:7826/1‑61 (MQ=255)
                                                                              gcAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACagag                                                                                                                                  <  2:7826/61‑1 (MQ=255)
                                                                                  cTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGAc                                                                                                                               <  2:612179/60‑1 (MQ=255)
                                                                                  cTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGAc                                                                                                                               >  1:612179/1‑60 (MQ=255)
                                                                                    tGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATACTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGTCAACGCCCGCCGCGCGGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACTGGGCg                                              >  2:515188/1‑139 (MQ=255)
                                                                                     gCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                               >  1:1283120/1‑121 (MQ=255)
                                                                                     gCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                               <  2:1283120/121‑1 (MQ=255)
                                                                                     gCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCgg                                             >  1:851339/1‑139 (MQ=255)
                                                                                     gCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCgg                                             >  1:163702/1‑139 (MQ=255)
                                                                                          tgtCGGTGACTCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACTGAGGACTTTTTTAACACGCCCAACAACGCCCGCGGCGCTGTTGTTTTGTTTCATTTGCAGGAACAGCGCATCGCCCATTACCGGGCGGTTGag                                        <  2:322235/139‑1 (MQ=255)
                                                                                          tgtCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                        <  2:414380/139‑1 (MQ=255)
                                                                                          tgtCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                        <  2:317113/139‑1 (MQ=255)
                                                                                          tgtCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                        <  2:163702/139‑1 (MQ=255)
                                                                                          tgtCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                        <  2:851339/139‑1 (MQ=255)
                                                                                          tgtCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                        <  1:419998/139‑1 (MQ=255)
                                                                                                         agagCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgttgtg                                                                                        <  2:418128/76‑1 (MQ=255)
                                                                                                         agagCAATACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgttgtg                                                                                        >  1:418128/1‑76 (MQ=255)
                                                                                                              aaTACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgttgtg                                                                                        <  1:1102273/71‑1 (MQ=255)
                                                                                                              aaTACCTGATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCtgttgtg                                                                                        >  2:1102273/1‑71 (MQ=255)
                                                                                                                     gATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                        <  1:959631/112‑1 (MQ=255)
                                                                                                                     gATGCTGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGag                                        >  2:959631/1‑112 (MQ=255)
                                                                                                                          tGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                                  >  2:395091/1‑49 (MQ=255)
                                                                                                                          tGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCccgccg                                                                                                  <  1:395091/49‑1 (MQ=255)
                                                                                                                          tGGTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtca        >  1:1151373/1‑139 (MQ=255)
                                                                                                                           ggTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACa                                                                    >  2:122757/1‑78 (MQ=255)
                                                                                                                           ggTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACa                                                                    <  1:122757/78‑1 (MQ=255)
                                                                                                                            gTAGCGGTGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCACTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtcatc      >  1:822338/1‑139 (MQ=255)
                                                                                                                            gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                               >  1:745208/1‑82 (MQ=255)
                                                                                                                            gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGa                                                               <  2:745208/82‑1 (MQ=255)
                                                                                                                            gTAGCGGCGACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCtcatc      >  2:100783/1‑139 (MQ=255)
                                                                                                                                cggcgACAGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACAt  <  1:663142/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                   
CGCCGACCTGGATAATGGCGTCGTAACCTGACAGCAGGCGCACGAAGTCGGTGAATCCCTGGGCGATGGCGATATTGCGCAGCTTGCCAGTGTCGGTGACGCGTGAGAGCAATACCTGATGCTGGTAGCGGCGACGGAGGACTTTTTTAACGCGCCCGACAACGCCCGCCGCGCTGTTGTGTTGTTTCATTTGCAGGAACAGCGGATCGCCCATTACCGGGCGGTTGAGCAGCCAGGAAGAACTGACCGGATAGCGGCTCATCACAT  >  NC_000913/2118029‑2118295

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