Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,089,570 T→C E96E (GAA→GAG pgaB ← poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine (PGA) N‑deacetylase outer membrane export lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,089,5700TC96.9% 101.6 / ‑6.4 32E96E (GAA→GAGpgaBpoly‑beta‑1,6‑N‑acetyl‑D‑glucosamine (PGA) N‑deacetylase outer membrane export lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (13/18);  total (14/18)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACGCAGTGCTGATGTCCGCACTGACATAAAACGCTGGTCGGCAGCTT  >  NC_000913/1089455‑1089679
                                                                                                                   |                                                                                                             
cagGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGAc                                                                                        <  1:148890/139‑1 (MQ=255)
cagGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGAc                                                                                        >  2:549582/1‑139 (MQ=255)
   gCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCgg                                                                                     >  1:741669/1‑139 (MQ=255)
     cACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGtt                                                                                   <  1:957116/139‑1 (MQ=255)
     cACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGtt                                                                                   >  2:649810/1‑139 (MQ=255)
           aaGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAAc                                                                             <  2:840931/139‑1 (MQ=255)
                  gaagaaTTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTtctc                                                                      >  2:592953/1‑139 (MQ=255)
                     gaaTTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTcgcg                                                                   <  1:592953/139‑1 (MQ=255)
                             aaGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGTCAgg                                                           >  1:306815/1‑139 (MQ=255)
                                  gcgGGTATAAAAACTCTGGTAGACGTCATCAAAAGTCAGCACTACAGCTTCTTACGATAGCGATGTTTATCCTCGATGTGCTTCTCGAATCTGAGAAATACTGACCGGTTGATTACAGTTCTCGCGCAGCCAGCGAAAt                                                      >  2:282006/1‑139 (MQ=255)
                                        ataAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTTCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCa                                                <  1:282006/139‑1 (MQ=255)
                                          aaaaaCTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGTCAGGCAAATTGTTCACg                                              <  2:306815/139‑1 (MQ=255)
                                          aaaaaCTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACg                                              <  1:649810/139‑1 (MQ=255)
                                          aaaaaCTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACg                                              <  1:549582/139‑1 (MQ=255)
                                                tctGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACGCAGTGc                                        >  1:45001/1‑139 (MQ=255)
                                                   ggTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGc                                                                                                 <  1:510487/79‑1 (MQ=255)
                                                   ggTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGc                                                                                                 >  2:510487/1‑79 (MQ=255)
                                                   ggTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGtt                                                                         >  1:731735/1‑103 (MQ=255)
                                                   ggTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGtt                                                                         <  2:731735/103‑1 (MQ=255)
                                                   ggTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACGCAGTGCTGa                                     >  2:1021594/1‑139 (MQ=255)
                                                           tcatcaAAAGTCAGCTGTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCGCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCAGGCAGTGCTGATGTCCGcg                             <  2:489024/139‑2 (MQ=255)
                                                           tcatcaAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACGCAGTGCTGATGTCCGCa                             <  2:2680/139‑1 (MQ=255)
                                                            catcaAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCa                                                                 <  1:1011475/102‑1 (MQ=255)
                                                            catcaAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCa                                                                 >  2:1011475/1‑102 (MQ=255)
                                                            catcaAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACGCAGTGCTGATGTCCGCAc                            <  2:680790/139‑1 (MQ=255)
                                                                           cTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACGCAGTGCTGATGTCCGCACTGACATAAAACGCTg             <  2:741669/139‑1 (MQ=255)
                                                                              cAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACGCAGTGCTGATGTCCGCACTGACATAAAACGCTGGTc          >  1:519168/1‑139 (MQ=255)
                                                                                      tCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAg                                                                                                  >  2:1044285/1‑43 (MQ=255)
                                                                                      tCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAg                                                                                                  <  1:1044285/43‑1 (MQ=255)
                                                                                      tCCGGTAGCGGTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACGCAGTGCTGATGTCCGCACTGACATAAAACGCTGGTCGGCAGCtt  <  2:519168/139‑1 (MQ=255)
                                                                                                gTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCg                                                                                      >  1:833795/1‑45 (MQ=255)
                                                                                                gTTTTCCTCCTCGATGTGCCTCACGAATTTGAGCAATACTGACCg                                                                                      <  2:833795/45‑1 (MQ=255)
                                                                                                                   |                                                                                                             
CAGGCCACTGGAAGGCCTGAAGAATTGGGAAGACGCGGGTATAAAAACTCTGGTAGCCGTCATCAAAAGTCAGCACTACAGCTTTTTCCGGTAGCGGTTTTCCTCCTCGATGTGCTTCACGAATTTGAGCAATACTGACCGGTTGATAACCGTTCTCGCGCAGCCAGGCAAATTGTTCACGCAGTGCTGATGTCCGCACTGACATAAAACGCTGGTCGGCAGCTT  >  NC_000913/1089455‑1089679

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