Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 939,928:1 +C 10.0% coding (5697/18933 nt) PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947939,9281.C10.0% 125.2 / 11.3 40coding (5697/18933 nt)PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base . (20/16);  new base C (2/2);  total (22/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.51e-01

GAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAG  >  NC_002947/939782‑940077
                                                                                                                                                    |                                                                                                                                                     
gAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGtg                                                                                                                                                    <  2:500249/150‑1 (MQ=255)
 aGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGtgc                                                                                                                                                   >  1:750429/1‑150 (MQ=255)
           tgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTc                                                                                                                                         <  1:677631/150‑1 (MQ=255)
                  tCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCgg                                                                                                                                  >  1:315478/1‑150 (MQ=255)
                          cGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGAc                                                                                                                          >  2:218429/1‑150 (MQ=255)
                            tGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCa                                                                                                                        <  1:446424/150‑1 (MQ=255)
                                ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGt                                                                                                                             <  2:694804/141‑1 (MQ=255)
                                ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGt                                                                                                                             >  1:694804/1‑141 (MQ=255)
                                ttGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGACCCAACTATCGGGCCCCCCGCACAGATATTGCGGTGATCGGCCGGCCACCGTCCGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCt                                                                                                                    >  1:278181/1‑150 (MQ=9)
                                     gCACCTCGCCGTCTATCTGTGGCGGGGCAAGGGTGATGCTGAAGCAGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGGGACCACGCTGCCGt                                                                                                               >  1:583116/1‑150 (MQ=255)
                                              cGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGa                                                                                                      <  1:605393/150‑1 (MQ=35)
                                                          gCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTtcgg                                                                                          >  2:794873/1‑150 (MQ=35)
                                                                cGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGGCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTg                                                                                    <  2:848147/150‑1 (MQ=25)
                                                                 gAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGc                                                                                   >  2:25432/1‑150 (MQ=35)
                                                                   gggTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCgg                                                                                 >  1:507513/1‑150 (MQ=35)
                                                                           cTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCg                                                                         <  1:218429/150‑1 (MQ=35)
                                                                            tGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCgg                                                                        <  1:25432/150‑1 (MQ=35)
                                                                             gAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                       >  2:180113/1‑150 (MQ=35)
                                                                                           cgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcgg                                                         >  2:660501/1‑150 (MQ=35)
                                                                                               aaCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGAt                                                     >  2:639389/1‑150 (MQ=35)
                                                                                                 cTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGAttt                                                   <  1:580534/150‑1 (MQ=35)
                                                                                                  tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAg                                                                              <  1:764673/122‑1 (MQ=35)
                                                                                                  tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAg                                                                              >  2:764673/1‑122 (MQ=35)
                                                                                                     cACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggt                                               <  2:507513/150‑1 (MQ=32)
                                                                                                            cccAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGGGGGCAGc                                        >  1:346093/1‑150 (MQ=17)
                                                                                                            cccAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGGGATTTGGGTGGGCAGc                                        >  1:188309/1‑150 (MQ=17)
                                                                                                             ccAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCACGGGCGGTGATATCCGGTGCTGTGATTTGGGTGGGCAGCg                                       >  2:400989/1‑150 (MQ=17)
                                                                                                              cAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTAAGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGTCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGa                                      >  1:699900/1‑150 (MQ=17)
                                                                                                                     gtggtgCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGATGATATCcggtgcggt                                                        >  1:690866/1‑125 (MQ=18)
                                                                                                                     gtggtgCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGATGATATCcggtgcggt                                                        <  2:690866/125‑1 (MQ=18)
                                                                                                                       tggtgCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAacac                              <  1:639389/150‑1 (MQ=32)
                                                                                                                       tggtgCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAacac                              >  2:725208/1‑150 (MQ=32)
                                                                                                                            ccGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                       >  1:210801/1‑104 (MQ=33)
                                                                                                                            ccGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                       <  2:210801/104‑1 (MQ=33)
                                                                                                                                  ccgggcAGCACCGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcgg                                                         >  2:776812/6‑112 (MQ=11)
                                                                                                                                  ccgggcAGCACCGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcgg                                                         <  1:776812/107‑1 (MQ=11)
                                                                                                                                   cgggCAGCACCGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCaa                  >  1:400435/1‑150 (MQ=11)
                                                                                                                                    gggcAGCACCGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAAt                 <  2:400435/147‑1 (MQ=11)
                                                                                                                                    ggcCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAAt                 <  2:750429/150‑1 (MQ=32)
                                                                                                                                                tACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTc     <  2:315478/150‑1 (MQ=32)
                                                                                                                                                   ggTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCGCCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCcag  <  2:583116/150‑1 (MQ=17)
                                                                                                                                                    |                                                                                                                                                     
GAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACGGTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATCACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCAG  >  NC_002947/939782‑940077

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