Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,361,724 C→G 15.0% R718P (CGG→CCG)  PP_2076 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,361,7240CG15.0% 102.2 / 10.1 40R718P (CGG→CCG) PP_2076hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/14);  new base G (3/3);  total (23/17)
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.91e-01

CTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCC  >  NC_002947/2361583‑2361862
                                                                                                                                             |                                                                                                                                           
ctTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCGGTACGGCGCGaa                                                                                                                                      <  1:260780/149‑2 (MQ=255)
              aaTCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCcaca                                                                                                                         >  4:182590/1‑148 (MQ=255)
                                  accacGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAg                                                                                                    >  3:3190/1‑149 (MQ=255)
                                               cAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCg                                                                                       >  8:201815/1‑149 (MQ=255)
                                                    ccccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGac                                                                                                                                      >  1:246709/1‑97 (MQ=255)
                                                    ccccGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGac                                                                                                                                      <  2:246709/97‑1 (MQ=255)
                                                                  aaaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAg                                                                                          <  4:178807/127‑1 (MQ=38)
                                                                  aaaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAg                                                                                          >  3:178807/1‑127 (MQ=38)
                                                                   aaCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCa                                                                     <  3:182590/147‑1 (MQ=37)
                                                                    aCTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGa                                                                  <  7:201815/149‑1 (MQ=38)
                                                                              cATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCtt                                                        >  5:19722/1‑149 (MQ=39)
                                                                                        aaCGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGa                                                                  <  7:185773/129‑1 (MQ=40)
                                                                                        aaCGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGa                                                                  >  8:185773/1‑129 (MQ=40)
                                                                                         aCGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCAt                                             <  6:19722/149‑1 (MQ=39)
                                                                                          cGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATc                                            >  5:9833/1‑149 (MQ=39)
                                                                                            cccTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCAc                                          >  3:262694/1‑149 (MQ=39)
                                                                                                 gagCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCa                                     <  4:3190/149‑1 (MQ=38)
                                                                                                          ctgctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCAGCGCATCACTGTCAGGCAGTTcg                           >  6:188104/1‑149 (MQ=255)
                                                                                                           tgctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgc                          >  8:36496/1‑149 (MQ=39)
                                                                                                            gctCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGTTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGtt                                                               <  6:17852/112‑1 (MQ=25)
                                                                                                             ctcctcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTcgc                          <  6:9833/147‑1 (MQ=39)
                                                                                                                 tcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGcc                                                                                        <  2:158059/82‑1 (MQ=38)
                                                                                                                 tcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGcc                                                                                        >  1:158059/1‑82 (MQ=38)
                                                                                                                 tcGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCtt                    >  6:219631/1‑149 (MQ=39)
                                                                                                                      cGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGaa                <  7:36496/148‑1 (MQ=39)
                                                                                                                       ggTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAg                                    <  1:247103/128‑1 (MQ=38)
                                                                                                                       ggTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAg                                    >  2:247103/1‑128 (MQ=38)
                                                                                                                          cacaCTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCa           >  3:66618/1‑149 (MQ=39)
                                                                                                                             ggtggtCGCGCGGTACGGCGCGAACACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCAt          >  8:66807/6‑147 (MQ=255)
                                                                                                                               tggtCGCGCGGTACGGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAg       >  3:257129/4‑148 (MQ=255)
                                                                                                                               tGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAg       >  6:77197/1‑148 (MQ=39)
                                                                                                                                 gtCGCGCGGTACGGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCAt    >  2:160787/2‑149 (MQ=255)
                                                                                                                                  tCGCGCGGTACGGCGCGAACACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATc   <  7:66807/149‑1 (MQ=255)
                                                                                                                                   cgcgCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCt                     >  6:88885/1‑130 (MQ=40)
                                                                                                                                   cgcgCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCt                     <  5:88885/130‑1 (MQ=40)
                                                                                                                                   cgcgCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATcc  >  5:255030/1‑149 (MQ=39)
                                                                                                                                   cgcgCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATcc  >  2:24907/1‑149 (MQ=39)
                                                                                                                                    gcgcGGTACGGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCCAGTTCGCGCGCTTTGAAGGCCATCAGCATc   <  1:160787/148‑1 (MQ=255)
                                                                                                                                    gcgcCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCAccc                                                   >  7:78088/1‑100 (MQ=40)
                                                                                                                                    gcgcCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCAccc                                                   <  8:78088/100‑1 (MQ=40)
                                                                                                                                             |                                                                                                                                           
CTTTTGTAATGCTCAATCGAGTCAGTGAAATCCCACCACGCGAGTACCAACCCCCCGGTGCCAGCTAAACTGGCGCCCCATAGTTTCAAACGCCCTAGAGCGGCCGCTGCTCCTCGTCCGGTCACACTGTTCGCGCCGTACCGCGCGACCACCTGGTCCACATAACTGGCTCCAACTTCGAAGCCCGCCGCAGCCGTAGTCATCGCGGCCGCCAGGAGTTCGGTCTTCACCCGCGCATCACTGTCAGGCAGTTCGCGCGCTTTGAAGGCCATCAGCATCC  >  NC_002947/2361583‑2361862

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