Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 213,354 T→A 12.5% 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,3540TA12.5% 103.4 / 8.1 40L6287Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (21/14);  new base A (3/2);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.05e-01

TCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTG  >  NC_002947/213207‑213503
                                                                                                                                                   |                                                                                                                                                     
ccggctcGTCGGTCACCGACGTGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTgg                                                                                                                                                     >  2:858324/7‑150 (MQ=255)
   gctcGTCGGTCACCGACGTGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCg                                                                                                                                                  >  1:843273/4‑150 (MQ=255)
   gCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCg                                                                                                                                                  >  1:918701/1‑150 (MQ=255)
    cACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGt                                                                                                                                                 <  2:755332/150‑1 (MQ=255)
    cACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACATCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGt                                                                                                                                                 >  2:626309/1‑150 (MQ=255)
      cGTCGGTCACCGACGTGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGa                                                                                                                                               <  2:843273/150‑1 (MQ=255)
             gACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAg                                                                                                                                        <  1:511304/150‑1 (MQ=255)
                cgacgaGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGag                                                                                                                                              >  1:235359/1‑140 (MQ=255)
                cgacgaGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCAGGTCGTGag                                                                                                                                              <  2:235359/141‑2 (MQ=255)
                 gacgaGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTTGCTCACGACCAGGTCGTGagc                                                                                                                                             >  2:984090/1‑139 (MQ=255)
                 gacgaGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTTGCTCACGACCAGGTCGTGagc                                                                                                                                             <  1:984090/141‑3 (MQ=255)
                               cccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGtcat                                                                                                                      >  1:604593/1‑150 (MQ=255)
                               cccAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTCAGCAACAATGCCCAGGtcat                                                                                                                      >  2:438775/1‑150 (MQ=255)
                                            ggCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGc                                                                                                         >  1:413745/1‑150 (MQ=255)
                                             gCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCg                                                                                                        >  1:612045/1‑150 (MQ=255)
                                                 cTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCCAACCGAGGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAAc                                                                                                    >  1:413253/1‑150 (MQ=255)
                                                 cTTGGTCAAGGTCACCATCACGGCCGACCAGACATCGGAGGCACCCGACGTTGACTCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAAc                                                                                                    <  1:805569/150‑1 (MQ=17)
                                                                         cGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTgc                                                                            <  1:914694/150‑1 (MQ=255)
                                                                                         gCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCTGGCGAGACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTc                                                            <  2:1024936/150‑1 (MQ=11)
                                                                                          ccGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCt                                                           <  1:438775/150‑1 (MQ=255)
                                                                                            gagaACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGGCTAc                                                         >  1:886274/1‑150 (MQ=255)
                                                                                                   tccaGCCGACGTTCACCGTGACCCTCGCCACCGCTGTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATg                                                  <  2:528785/147‑1 (MQ=17)
                                                                                                   tCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATg                                                  <  2:413745/150‑1 (MQ=255)
                                                                                                         cGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGcc                                            >  2:604473/1‑150 (MQ=255)
                                                                                                             gTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGAt                                        >  2:79953/1‑150 (MQ=255)
                                                                                                                    gTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTg                                 <  2:429381/150‑1 (MQ=255)
                                                                                                                    gTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTg                                 >  2:803548/1‑150 (MQ=255)
                                                                                                                     tGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTgg                                >  2:860296/1‑150 (MQ=255)
                                                                                                                      gCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTggg                               <  1:858324/150‑1 (MQ=255)
                                                                                                                          aTCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCAtt                           >  1:281284/1‑150 (MQ=255)
                                                                                                                               cacCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTc                      >  2:600675/1‑150 (MQ=255)
                                                                                                                               cacCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTc                      <  1:856332/150‑1 (MQ=255)
                                                                                                                                   gCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCa                  <  1:600675/150‑1 (MQ=255)
                                                                                                                                         ggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGAGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGAc            >  1:295503/3‑150 (MQ=255)
                                                                                                                                           tCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAgg                                                                      <  1:436586/90‑1 (MQ=255)
                                                                                                                                           tCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAgg                                                                      >  2:436586/1‑90 (MQ=255)
                                                                                                                                           tCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGGCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGc          >  2:345915/1‑150 (MQ=255)
                                                                                                                                             aCGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAGGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAACCTGGGCATTAATTCGGCAGTGGACGcgc        >  2:524183/1‑148 (MQ=255)
                                                                                                                                             aCGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCa        >  2:939266/1‑150 (MQ=255)
                                                                                                                                               gACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACt      >  2:332072/1‑150 (MQ=255)
                                                                                                                                                   tGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTg  <  2:938921/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCAGCACGTCGGTGACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACTTGGTCAAGGTCACCATCACGGCCGACCAGACCTCGGTGGCCGAGAACGTCAAGCCGACGTTCACCGTGCACATCAACACCGCCTTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCATCATCAAGGCTGGCGAAACCAGCGCGCCATACACCCACGCTGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATTAATTCGGCAGTGGACGCCACTGGTG  >  NC_002947/213207‑213503

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: