Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 213,354 T→A 20.9% L6287Q (CTG→CAG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947213,3540TA20.9% 86.2 / 20.7 43L6287Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (17/17);  new base A (7/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

AGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGC  >  NC_002947/213209‑213495
                                                                                                                                                 |                                                                                                                                             
aGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGCCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTc                                                                                                                                           <  2:67374/150‑1 (MQ=255)
       cGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGAccc                                                                                                                                    <  2:385714/150‑1 (MQ=255)
        ggTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGGTCACGACCAGGTCGTGAgcca                                                                                                                                   >  1:144620/1‑149 (MQ=255)
            aCCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGGTCACGACCAGGTCGTGAgccaaggc                                                                                                                               <  2:144620/150‑6 (MQ=18)
                   aGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATg                                                                                                                        <  1:676567/150‑1 (MQ=255)
                        ggTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAg                                                                                                                   <  1:356713/150‑1 (MQ=255)
                              ccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGtcatc                                                                                                             >  2:425224/1‑150 (MQ=255)
                                         aGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTgg                                                                                                  <  1:425224/150‑1 (MQ=255)
                                                ttGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAAcc                                                                                           >  2:658523/1‑150 (MQ=255)
                                                 tGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCGACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCa                                                                                          <  2:773759/150‑1 (MQ=255)
                                                   gTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAgc                                                                                        <  2:217388/150‑1 (MQ=255)
                                                    tCAAGGTCACCATCACGGCCCACTAGACCTCGGTGGCCGAGGACGTCATGCCGACGTTCACCGTGAACATCTACACCTCCTTGGCCAACGACCTGTTCGTGACCATGAGCAAGAATACCCAGACCGTCATCAAGGCTCACTAAACCAgca                                                                                       >  2:109603/1‑149 (MQ=12)
                                                      aaGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAgcgcg                                                                                     >  1:97678/1‑150 (MQ=255)
                                                            aCCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATac                                                                               >  2:725820/1‑150 (MQ=255)
                                                                    ggCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGc                                                                       >  2:229576/1‑150 (MQ=255)
                                                                                 tCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGAt                                                          <  1:725820/150‑1 (MQ=255)
                                                                                      ggCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGt                                                     <  1:229576/150‑1 (MQ=255)
                                                                                       gCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTgg                                                                                                  >  1:446676/1‑104 (MQ=255)
                                                                                       gCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTgg                                                                                                  <  2:446676/104‑1 (MQ=255)
                                                                                         cGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGa                                                           >  1:484684/1‑141 (MQ=255)
                                                                                         cGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGa                                                           <  2:484684/141‑1 (MQ=255)
                                                                                         cGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTc                                                    <  1:238241/148‑1 (MQ=255)
                                                                                         cGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTc                                                    >  2:238241/1‑148 (MQ=255)
                                                                                         cGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                  >  2:443904/1‑150 (MQ=255)
                                                                                         cGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                  >  2:818483/1‑150 (MQ=255)
                                                                                         cGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCACACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTa                                                  >  2:106204/1‑150 (MQ=255)
                                                                                          gagaACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTAc                                                 <  1:139551/150‑1 (MQ=255)
                                                                                                gTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACg                                                      <  2:732127/139‑1 (MQ=255)
                                                                                                gTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACg                                                      >  1:732127/1‑139 (MQ=255)
                                                                                                 tCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGCCCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCTCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATg                                          <  2:97678/150‑1 (MQ=255)
                                                                                                     gCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTgg                                      >  2:722674/1‑150 (MQ=255)
                                                                                                                ccGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGcc                           >  1:441827/1‑150 (MQ=255)
                                                                                                                 cGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCAAGATCAGCCt                          <  1:818483/150‑1 (MQ=255)
                                                                                                                  gTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTg                         <  1:658523/150‑1 (MQ=255)
                                                                                                                            acacCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAAtt               >  2:717331/1‑150 (MQ=255)
                                                                                                                               ccGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCgg            >  2:449468/1‑150 (MQ=255)
                                                                                                                                   tcagggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGt        >  1:29305/7‑150 (MQ=255)
                                                                                                                                     agggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTgg      >  1:609507/5‑150 (MQ=255)
                                                                                                                                      gggtCACGACCAGGTCGTGACCCTGAGCAACAATACCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGa     >  1:379304/4‑150 (MQ=255)
                                                                                                                                        gtCACGACCAGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCaa                                                                <  2:74745/88‑1 (MQ=255)
                                                                                                                                        gtCACGACCAGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCaa                                                                >  1:74745/2‑89 (MQ=255)
                                                                                                                                        gtCACGACCAGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACg   >  2:393622/2‑150 (MQ=255)
                                                                                                                                         tCACGACCAGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGc  >  1:342742/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
AGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGC  >  NC_002947/213209‑213495

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