Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 202,318 Δ1 bp 15.4% coding (7824/26049 nt) PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947202,3180C.15.4% 62.8 / 8.5 26coding (7824/26049 nt)PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (13/9);  new base . (2/2);  total (15/11)
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.86e-01

CACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCC  >  NC_002947/202186‑202462
                                                                                                                                      |                                                                                                                                                
cACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTacac                                                                                                                                 >  2:323594/1‑150 (MQ=32)
        cTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCAccc                                                                                                                         >  1:73271/1‑150 (MQ=32)
           tgGATGTAGACGGTCACATCTTCGAGCCCCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGa                                                                                                                      <  2:73271/150‑1 (MQ=12)
                 tAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACcaaca                                                                                                                <  2:14099/150‑1 (MQ=32)
                          acacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCAc                                                                                                                                       >  2:378740/1‑118 (MQ=33)
                          acacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCAc                                                                                                                                       <  1:378740/118‑1 (MQ=33)
                           cacCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCt                                                                                                      >  1:115426/1‑150 (MQ=35)
                               ttCGAGAACCTGCAACTGGGCGGCCACGCTTCGGTTCAAGTGACTGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCg                                                                                                  >  1:9793/1‑150 (MQ=25)
                                              cTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGCCGTTTCCCGTTCGATCCAGGTGG‑‑ACCGAAGGCGGGGAGATCACCTACAACACCACCGTTATAGACATAGTCGGTCTGCCGATCAACAACCACTCg                                                                                   <  1:323594/150‑1 (MQ=2)
                                                    ggcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTg                                                                             >  1:169181/1‑150 (MQ=35)
                                                     gcAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGt                                                                            >  1:57417/1‑150 (MQ=35)
                                                                      gTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGAc                                                           <  2:115426/150‑1 (MQ=35)
                                                                        gACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACgg                                                         <  1:314526/150‑1 (MQ=35)
                                                                        gACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACgg                                                         <  2:57417/150‑1 (MQ=34)
                                                                           cGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACgg                                                                                                         >  1:244644/1‑99 (MQ=38)
                                                                           cGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACgg                                                                                                         <  2:244644/99‑1 (MQ=38)
                                                                              cacCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGcaccac                                                   >  2:314791/1‑150 (MQ=35)
                                                                                  cTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGt                                                                                                                                               <  1:302174/54‑1 (MQ=38)
                                                                                  cTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGt                                                                                                                                               >  2:302174/1‑54 (MQ=38)
                                                                                   tCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACcgtcg                                              >  1:147392/1‑150 (MQ=38)
                                                                                                               tcaccgccTTCGGTCACCGAAGGGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGcac                           <  1:347778/136‑1 (MQ=12)
                                                                                                               tcaccgccTTCGGTCACCGAAGGGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGcac                           >  2:347778/7‑142 (MQ=255)
                                                                                                                 accgccTTCGGTCACCGAAGGGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAg                                                                    <  2:46608/95‑1 (MQ=255)
                                                                                                                 accgccTTCGGTCACCGAAGGGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAg                                                                    >  1:46608/5‑99 (MQ=255)
                                                                                                                       ttCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCAc            >  1:101445/1‑150 (MQ=32)
                                                                                                                                 gAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGcc  >  2:272152/1‑150 (MQ=11)
                                                                                                                                      |                                                                                                                                                
CACGTCCGCTGTGGATGTAGACGGTCACACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCCTCGACGAGGTGGTGGCGAAGCTGACCG‑‑CGACCCCTTCGGTCACCGAAGGCGGTGAGATCACCTACACCATCACCCTGACCAACAAAGACGGTCTGCCGATCAACAACCACTCGGAGCTGTACTTCAAGCTGACCGACGGCACCACCGTCGTCGTGGCAGCCAACAGCACCACCGGTTCGGCCACTGCAACCGCC  >  NC_002947/202186‑202462

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