Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,006,478 C→T 11.3% P470S (CCG→TCG)  PP_0867 → FecA‑like outer membrane receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,006,4780CT11.3% 91.0 / 3.2 37P470S (CCG→TCG) PP_0867FecA‑like outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base C (15/16);  new base T (2/2);  total (19/18)
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.93e-01

GCGCAACTACCACACCTTCGCCATTGAACCCCGTGTGTCGCGCATCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGACCATCACCCCTGGCGTGCGCTTCGAGCACATCAATA  >  NC_002947/1006330‑1006619
                                                                                                                                                    |                                                                                                                                             
gcgcAACTACCACACCTTCGCCATTGAACCCCGTGTGTCGCGCATCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCATGGCGCTGATCGACAACGTGc                                                                                                                                               >  4:30617/1‑149 (MQ=255)
         ccACACCTTCGCCATTGAACCCCGTGTGTCGCGCATCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACg                                                                                                                                          <  5:89579/145‑1 (MQ=255)
         ccACACCTTCGCCATTGAACCCCGTGTGTCGCGCATCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGAAGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACg                                                                                                                                          >  6:89579/1‑145 (MQ=255)
                   gCCATTGAACCCCGTGTGTCGCGCATCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACtt                                                                                                                             <  2:330734/148‑1 (MQ=255)
                          aaCCCCGTGTGTCGCGCATCTTCTTTGCCGGCCCAACCAACCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTAGACTGGCGCTGATCGACAAAGTGACAACGGTGACCCCGACTTCCGAAGGc                                                                                                                     >  6:177243/1‑149 (MQ=255)
                             cccGTGTGTCGCGCATCTTCTTTGCCGGCACAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGAGCGAACACTCGACTCGACAGGCGCAGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCAc                                                                                                                  >  1:129276/1‑149 (MQ=255)
                             cccGGGTGGCGCGCAGCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCAc                                                                                                                  <  8:92931/149‑1 (MQ=255)
                                    gtCGCGCATCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTcc                                                                                                           <  3:30617/149‑1 (MQ=255)
                                       gcgcATCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAgg                                                                                                        >  2:315070/1‑149 (MQ=255)
                                           aTCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGAACCCGACTTCCGACGGCCACGTGTTCCAGGACCg                                                                                                    >  4:1670/1‑149 (MQ=255)
                                                cttTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGGTGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCg                                                                                               >  8:213884/1‑149 (MQ=255)
                                                cttTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCg                                                                                               >  8:217755/1‑149 (MQ=255)
                                                     ccGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTgg                                                                                           >  3:236744/1‑148 (MQ=255)
                                                          ccAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGTTGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGa                                                                                     <  7:213884/149‑1 (MQ=255)
                                                          ccAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGa                                                                                     <  7:217755/149‑1 (MQ=255)
                                                           cAACCACCCAGGAAGTCATTGTCGGTTACCGTTACCTGAAGGAAGCGATGCTCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAg                                                                                    <  3:171907/149‑1 (MQ=255)
                                                             cccaccCCGGAAGTCATTGTCGGTTACCGTTACCTGTAGGTAGCGTTGCGCGAACAGTCGACTCGCCTGTCGCTGTTCGACAACGTGCCGACGGTGACCCGGACTTCCGACGGGCACGTGTTCCAGGACCGCTGCGGTGGCACCGAGGc                                                                                  <  6:82880/148‑1 (MQ=255)
                                                                 cccAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAgc                                                                              <  1:315070/149‑1 (MQ=255)
                                                                    aGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGcc                                                                           >  1:216990/1‑149 (MQ=255)
                                                                          tCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCtactac                                                                     <  4:236744/149‑1 (MQ=255)
                                                                             gtgtCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATc                                                                  <  5:136530/149‑1 (MQ=255)
                                                                                  ggTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCTACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATcgacga                                                             <  4:200140/149‑1 (MQ=255)
                                                                                              cttacggaagCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGACTACTACATCGACGACAAGATCGACGt                                                 <  2:129276/144‑1 (MQ=255)
                                                                                               tGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGTCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTg                                                >  4:296265/1‑149 (MQ=255)
                                                                                               tGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGTCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTg                                                >  4:291328/1‑149 (MQ=255)
                                                                                                           tGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGTCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGAcc                                    <  3:296265/149‑1 (MQ=255)
                                                                                                           tGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGAcc                                    <  2:216990/149‑1 (MQ=255)
                                                                                                            gcgcGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGTCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGAcc                                    <  3:291328/148‑1 (MQ=255)
                                                                                                              gcgAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGCGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTg                                                >  7:222773/1‑134 (MQ=255)
                                                                                                              gcgAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGCGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTg                                                <  8:222773/134‑1 (MQ=255)
                                                                                                                      tCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGACCATCACCCCTgg                         >  2:296782/1‑149 (MQ=255)
                                                                                                                        gACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCGGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGACCATCACCCCTGGCg                       >  3:241587/1‑149 (MQ=255)
                                                                                                                        gACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGACCATCACCCCTGGCg                       >  7:57782/1‑149 (MQ=255)
                                                                                                                         aCTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGACCATCACCCCTGGCGt                      >  7:91020/1‑149 (MQ=255)
                                                                                                                              ccTGGCGCTGATCGACAACGTGACGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCCCCGAGGCCAGCGCCAACTACATCGACGACAAGGGCGCCGTGGGCAAGTGGCCCATCACCCCTGGCGAgcg                   >  7:211468/1‑147 (MQ=255)
                                                                                                                                     cgagtCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGACCATCACCCCTGGCGTGCGCTCCGAGca          <  8:8629/145‑1 (MQ=255)
                                                                                                                                         tCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGAACATCACCCCTGGCGTGCGCTTCGAGCACAt       >  7:145390/1‑148 (MQ=255)
                                                                                                                                             cAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGACCATCACCCCTGGCGTGCGCTTCGAGCACATCAATa  >  5:168815/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GCGCAACTACCACACCTTCGCCATTGAACCCCGTGTGTCGCGCATCTTCTTTGCCGGCCCAACCACCCAGGAAGTCAGTGTCGGTTACCGTTACCTGAAGGAAGCGATGCGCGAACAGTCGACTCGCCTGGCGCTGATCGACAACGTGCCGACGGTGACCCCGACTTCCGACGGCCACGTGTTCCAGGACCGCAGCGGTGGCACCGAGGCCAGCGCCTACTACATCGACGACAAGATCGACGTGGGCAACTGGACCATCACCCCTGGCGTGCGCTTCGAGCACATCAATA  >  NC_002947/1006330‑1006619

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: ‑