Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 212,043 C→G 12.1% T5850S (ACC→AGC)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947212,0430CG12.1% 76.5 / 4.9 33T5850S (ACC→AGC) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (14/15);  new base G (2/2);  total (16/17)
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.28e-01

GGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCG  >  NC_002947/211897‑212192
                                                                                                                                                  |                                                                                                                                                      
ggTTACCGAGGAGACAGGTACGCCGGGAAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGCCAACCAGCCGCTGGGTCACGACCAGGTCGTGAgcca                                                                                                                                                     <  2:169441/150‑2 (MQ=11)
  ttACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                                   >  2:1121958/1‑150 (MQ=35)
        acgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGAcc                                                                                                                                                       >  2:29541/1‑140 (MQ=37)
        acgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGAcc                                                                                                                                                       <  1:29541/140‑1 (MQ=37)
          gaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAt                                                                                                                                           <  1:678047/150‑1 (MQ=35)
               aGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAATGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCa                                                                                                                                      >  1:201012/1‑150 (MQ=25)
               aGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCa                                                                                                                                      <  2:188171/150‑1 (MQ=35)
                ggTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAg                                                                                                                                     <  1:613035/150‑1 (MQ=35)
                     gCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACAAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGCGCAACAATGCCCAggaccc                                                                                                                                >  1:360041/1‑146 (MQ=255)
                                                                    aGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAgg                                                                                <  1:228677/150‑1 (MQ=32)
                                                                    aGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAgg                                                                                >  1:539344/1‑150 (MQ=32)
                                                                        ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtga                                                                            >  2:665575/1‑150 (MQ=32)
                                                                        ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtga                                                                            <  2:795224/150‑1 (MQ=32)
                                                                         tCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtgat                                                                           <  2:34219/150‑1 (MQ=32)
                                                                            gTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTAAACCAGCCGCTGGCTCACGAACTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGACGGTGAAACCCGCGCGACGTACACCCACGCTGCGCAGGGTGATGAc                                                                        >  2:1003137/1‑150 (MQ=11)
                                                                                 cGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAATTCTGCCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTa                                                                   >  1:374559/1‑150 (MQ=11)
                                                                                  gagaGCGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTAc                                                                  <  1:718882/150‑1 (MQ=17)
                                                                                  gagaACGTCAAACCGACCTTCACTACTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTAc                                                                  >  1:1261588/1‑150 (MQ=17)
                                                                                        gTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtgatga                                                                         <  1:851187/137‑1 (MQ=33)
                                                                                        gTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtgatga                                                                         >  2:851187/1‑137 (MQ=33)
                                                                                                   cTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACTCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGAt                                                 >  2:376222/1‑150 (MQ=25)
                                                                                                     tCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCa                                               <  1:572658/150‑1 (MQ=35)
                                                                                                         tGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCt                                           >  1:839703/1‑150 (MQ=35)
                                                                                                              aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTa                                      <  2:859721/150‑1 (MQ=35)
                                                                                                                    aCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACt                                >  1:1005300/1‑150 (MQ=35)
                                                                                                                    aCCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAAGGCGCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCTCGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACt                                <  2:360041/150‑1 (MQ=255)
                                                                                                                               ggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCAACGCTGCGCAGGGTGATGACGTCTACAACGATGCTgg                                                       >  1:638703/3‑116 (MQ=11)
                                                                                                                               ggtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCAACGCTGCGCAGGGTGATGACGTCTACAACGATGCTgg                                                       <  2:638703/114‑1 (MQ=11)
                                                                                                                                gtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAgg                                                                                >  2:1336664/2‑90 (MQ=255)
                                                                                                                                gtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAgg                                                                                <  1:1336664/89‑1 (MQ=255)
                                                                                                                                        cTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGt            <  1:665575/150‑1 (MQ=35)
                                                                                                                                        cTGGTCGTGACCCTGAACAACAATGCCCAGGTCACCCTCAAGCCCGTTGAACCCAGCGCGCCGTACACCGACGCTGCGCAGGGTGACGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGt            >  2:923878/1‑150 (MQ=11)
                                                                                                                                           gTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGAt                                                 >  2:1290494/1‑110 (MQ=38)
                                                                                                                                           gTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGAt                                                 <  1:1290494/110‑1 (MQ=35)
                                                                                                                                            tCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGa        <  2:839703/150‑1 (MQ=35)
                                                                                                                                                 aCCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTc   <  2:1261588/150‑1 (MQ=35)
                                                                                                                                                  cccTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCg  >  1:1048096/1‑150 (MQ=35)
                                                                                                                                                  |                                                                                                                                                      
GGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGTGCGACCTTCG  >  NC_002947/211897‑212192

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: