Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,131,667 T→C D135G (GAC→GGC)  wcaC ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,131,6670TC94.9% 121.1 / ‑2.3 39D135G (GAC→GGC) wcaCputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (19/18);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

TGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGA  >  NC_000913/2131533‑2131801
                                                                                                                                      |                                                                                                                                      
tGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTg                                                                                                                                    <  1:935676/139‑1 (MQ=255)
                       cgATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGt                                                                                                             >  1:238503/1‑139 (MQ=255)
                                    ccGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGt                                                                                                             <  2:1308893/126‑1 (MQ=255)
                                    ccGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGt                                                                                                             <  2:516008/126‑1 (MQ=255)
                                    ccGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGt                                                                                                             >  1:1308893/1‑126 (MQ=255)
                                    ccGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGt                                                                                                             >  1:516008/1‑126 (MQ=255)
                                                         tCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAgcg                                                                           >  1:892808/1‑139 (MQ=255)
                                                            ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGc                                                                                                                                <  2:314867/83‑1 (MQ=255)
                                                            ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGc                                                                                                                                >  1:314867/1‑83 (MQ=255)
                                                            ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGt                                                                                                                              >  1:757247/1‑85 (MQ=255)
                                                            ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGt                                                                                                                              <  2:757247/85‑1 (MQ=255)
                                                            ggCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcac                                                                        >  1:292788/1‑139 (MQ=255)
                                                              cATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcacca                                                                      <  2:892808/139‑1 (MQ=255)
                                                                                      ccTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCCGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGc                                              >  2:825869/1‑139 (MQ=255)
                                                                                        ttCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcacca                                                                      <  1:449278/113‑1 (MQ=255)
                                                                                        ttCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcacca                                                                      >  2:449278/1‑113 (MQ=255)
                                                                                          cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTccagacc                                                                                                                       >  1:1104988/1‑62 (MQ=255)
                                                                                          cGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTccagacc                                                                                                                       <  2:1104988/62‑1 (MQ=255)
                                                                                                  tCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                    <  2:1154186/105‑1 (MQ=255)
                                                                                                  tCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGATTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                    >  1:1154186/1‑105 (MQ=255)
                                                                                                            cgcAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGcc                                                    >  2:769366/1‑111 (MQ=255)
                                                                                                            cgcAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGcc                                                    <  1:769366/111‑1 (MQ=255)
                                                                                                            cgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCATTTTTTCGCAAAAGCGCACCACGCTCTTAAg                                                            <  2:22842/103‑1 (MQ=255)
                                                                                                            cgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCATTTTTTCGCAAAAGCGCACCACGCTCTTAAg                                                            >  1:22842/1‑103 (MQ=255)
                                                                                                            cgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTgcagca                                     <  2:1282168/126‑1 (MQ=255)
                                                                                                            cgcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTgcagca                                     >  1:1282168/1‑126 (MQ=255)
                                                                                                             gcAGCGTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAgg                       <  2:238503/139‑1 (MQ=255)
                                                                                                                  gTCCGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGaccac                  <  2:120436/139‑1 (MQ=255)
                                                                                                                    ccGGTAACACTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCg                <  2:292788/139‑1 (MQ=255)
                                                                                                                          acacTCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTCTCTCAAAAGCGCACCACGCTCTTAGGATTTAGCCTGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCg          <  1:825869/139‑1 (MQ=255)
                                                                                                                              tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcacca                                                                      >  1:1002189/1‑75 (MQ=255)
                                                                                                                              tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGcacca                                                                      <  2:1002189/75‑1 (MQ=255)
                                                                                                                              tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                    <  1:313434/77‑1 (MQ=255)
                                                                                                                              tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACg                                                                    >  2:313434/1‑77 (MQ=255)
                                                                                                                              tCCAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCtgt      >  1:994213/1‑139 (MQ=255)
                                                                                                                                cAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCg          <  2:1152070/133‑1 (MQ=255)
                                                                                                                                cAGTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCg          >  1:1152070/1‑133 (MQ=255)
                                                                                                                                  gTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGa  >  2:522313/1‑139 (MQ=255)
                                                                                                                                  gTGGCCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGa  >  2:552008/1‑139 (MQ=255)
                                                                                                                                      |                                                                                                                                      
TGCCCGCCACCAGTTGGTGTGCGCGATCAATCTTCACCGGCGGATAGTTATTTAAGGTCGGGCATTTCTGGCAGCCTGTTTTCCAGCCTTCGCAACCGTCGGTAAAGGCGCAGCGTCCGGTAACACTCCAGTGGTCGTGCAGCGTCCAGACCAGAGTGACGTCCGGTTTATGGTTTTTCACTTTTTCGCAAAAGCGCACCACGCTCTTAAGATTTAGCCAGTAGCTGTGCAGCACATGAAAATGCAGGACCACCGGACCCGCTGTGCGA  >  NC_000913/2131533‑2131801

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