Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 212,034 T→A 16.7% 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,0340TA16.7% 78.1 / 11.4 36L5847Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (15/15);  new base A (2/4);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

CGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGT  >  NC_002947/211893‑212182
                                                                                                                                             |                                                                                                                                                     
cGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGa                                                                                                                                               <  2:314779/150‑1 (MQ=34)
 gTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGag                                                                                                                                              <  1:91739/150‑2 (MQ=21)
    ggTTACCGAGGAGACAGGTACGCCGGGAAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGCCAACCAGCCGCTGGGTCACGACCAGGTCGTGAgcca                                                                                                                                           <  2:169441/150‑2 (MQ=11)
      ttACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                         >  2:1121958/1‑150 (MQ=35)
          cgacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                   >  1:102802/1‑136 (MQ=37)
          cgacgaGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                   <  2:102802/136‑1 (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:539344/1‑150 (MQ=32)
                                                                        aGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAgg                                                                      <  1:228677/150‑1 (MQ=32)
                                                                            ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtga                                                                  <  2:795224/150‑1 (MQ=32)
                                                                            ttCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGtga                                                                  >  2:665575/1‑150 (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                                                                      <  1:1336664/89‑1 (MQ=255)
                                                                                                                                    gtCACGACCAGGTCGTGAGCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAgg                                                                      >  2:1336664/2‑90 (MQ=255)
                                                                                                                                            cTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGt  <  1:665575/150‑1 (MQ=35)
                                                                                                                                            cTGGTCGTGACCCTGAACAACAATGCCCAGGTCACCCTCAAGCCCGTTGAACCCAGCGCGCCGTACACCGACGCTGCGCAGGGTGACGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGt  >  2:923878/1‑150 (MQ=11)
                                                                                                                                             |                                                                                                                                                     
CGTCGGTTACCGACGAGCCAGGTACGCCGGGCAACGAAGGCGACCTGGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAACCGACCTTCACTGTTCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAATGCCCAGGTCACCATCAAGGCCGGTGAAACCAGCGCGCCGTACACCCACGCTGCGCAGGGTGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGTATCAACTCGGCGGTGGACGCCACTGGT  >  NC_002947/211893‑212182

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: