Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 212,034 T→A 9.3% L5847Q (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_002947212,0340TA 9.3% 100.5 / 4.7 43L5847Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (16/23);  new base A (2/2);  total (18/25)
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.92e-01

GGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTG  >  NC_002947/211885‑212183
                                                                                                                                                     |                                                                                                                                                     
ggTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCt                                                                                                                                                       <  1:39008/150‑1 (MQ=35)
  tCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCACGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTgg                                                                                                                                                     <  2:1491787/150‑1 (MQ=25)
  tCAGAACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTgg                                                                                                                                                     >  2:410498/1‑150 (MQ=25)
      cACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGt                                                                                                                                                 <  1:878038/150‑1 (MQ=35)
       aCGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTg                                                                                                                                                >  2:193800/1‑150 (MQ=25)
         gTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGag                                                                                                                                              <  1:728976/150‑2 (MQ=25)
           cGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGAccc                                                                                                                                            <  1:276069/150‑1 (MQ=35)
            ggTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCt                                                                                                                                           <  1:607709/150‑1 (MQ=35)
            gcttACCGACCAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGATCACCAGCGCGGCCGTCCATACTTCGGTATCCGAGAACGTCAAACCGACCCTCACGGTTCACGTCAACCAGCCTCTGGCTCACGACCTGGTCGTGACCCt                                                                                                                                           <  2:909128/148‑1 (MQ=11)
              ttACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGagccagc                                                                                                                                         <  1:193800/150‑7 (MQ=255)
                  cgacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCCAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                                     >  1:924593/1‑150 (MQ=21)
                   gacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                    >  2:1473728/1‑150 (MQ=34)
                   gacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                    >  1:564554/1‑150 (MQ=34)
                            ggTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGGCCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAg                                                                                                                           <  2:1145626/150‑1 (MQ=25)
                                 gCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCAc                                                                                                                      <  1:161556/150‑1 (MQ=33)
                                    gggCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCAt                                                                                                                   >  1:743020/1‑150 (MQ=31)
                                                                                   cTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtg                                                                    <  2:564554/150‑1 (MQ=32)
                                                                                    ttCGGTAGCCGAGCACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtga                                                                   >  1:288896/1‑150 (MQ=17)
                                                                                     tCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtgat                                                                  <  1:419666/150‑1 (MQ=32)
                                                                                     tCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtgat                                                                  <  2:482728/150‑1 (MQ=32)
                                                                                                   cGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGa                                                    <  1:410498/150‑1 (MQ=32)
                                                                                                    gTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGAt                                                   <  1:1599916/150‑1 (MQ=33)
                                                                                                     tCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTgcg                                                                           <  1:1892696/125‑1 (MQ=11)
                                                                                                     tCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTgcg                                                                           >  2:1892696/1‑125 (MQ=11)
                                                                                                        aaCNGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTg                                               >  2:1189091/1‑150 (MQ=25)
                                                                                                        aaCCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCGGTTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTg                                               <  2:1591632/150‑1 (MQ=25)
                                                                                                          ccGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGc                                             <  1:211848/150‑1 (MQ=35)
                                                                                                           cGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGcc                                            >  2:377452/1‑150 (MQ=35)
                                                                                                           cGACCTTCACTGTTCACGTCAACCAAACGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGcc                                            >  1:52250/1‑150 (MQ=25)
                                                                                                                ttCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTAGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATc                                       <  1:108565/150‑1 (MQ=25)
                                                                                                                       ttCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTgg                                <  2:923951/150‑1 (MQ=35)
                                                                                                                               aaCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAAc                        >  2:1238418/1‑150 (MQ=35)
                                                                                                                                       gCTGGCTCACGAGCTGGGCGTGACCGTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGt                <  2:924593/150‑1 (MQ=255)
                                                                                                                                       gCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGt                <  2:705094/150‑1 (MQ=35)
                                                                                                                                       gCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGt                <  2:743020/150‑1 (MQ=35)
                                                                                                                                       gCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGt                <  1:1189091/150‑1 (MQ=35)
                                                                                                                                           ggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGAc            >  2:1885481/3‑150 (MQ=255)
                                                                                                                                             tCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGc          >  2:1151303/1‑150 (MQ=35)
                                                                                                                                             tCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGc          >  1:1605544/1‑150 (MQ=35)
                                                                                                                                             tCACGACCTCGGCGTGACCCTGAGCAACAATGACCAGGTCACCATCAAGGCAGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGc          <  2:912119/150‑1 (MQ=255)
                                                                                                                                               aCGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCa        >  2:370132/1‑150 (MQ=35)
                                                                                                                                               aCGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTAAAACTGTGCTGCGCAGGGTGAAGCCGCTTATAAGCATTCGTGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGccc        >  1:1869321/1‑149 (MQ=11)
                                                                                                                                                 gACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCACAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACt      >  1:1411704/1‑150 (MQ=17)
                                                                                                                                                   ccTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTgg    <  1:1187213/150‑1 (MQ=35)
                                                                                                                                                     tGGNCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTg  >  2:1193752/1‑150 (MQ=25)
                                                                                                                                                     |                                                                                                                                                     
GGTCAGCACGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTG  >  NC_002947/211885‑212183

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