Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,306 T→C 28.3% N604N (AAT→AAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,3060TC28.3% 34.5 / 17.3 29N604N (AAT→AACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (9/11);  new base C (5/3);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGT  >  NC_002947/196173‑196440
                                                                                                                                     |                                                                                                                                      
aGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATcc                                                                                                                                      >  2:398810/1‑136 (MQ=39)
   aaCTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGc                                                                                                                                   >  2:271802/1‑136 (MQ=37)
                 ccGACCCGGGCTCACACCACGAGCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGa                                                                                                                     <  1:133884/136‑1 (MQ=255)
                      cccGCCTCACACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGACGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTg                                                                                                                <  1:204052/135‑1 (MQ=255)
                              aGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAAACCTGCCCAGACGCCGGTAACAGTGAccactgg                                                                                                        >  1:145344/1‑131 (MQ=255)
                                  cccGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAAc                                                                                                    <  1:156586/134‑1 (MQ=255)
                                        cAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTCCCCACACCCCCGTTACACTGCCCCTGTCCACGGCCaaa                                                                                             >  2:200413/1‑136 (MQ=11)
                                             aCTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCACGGCCAAACCAtt                                                                                        >  2:15032/1‑136 (MQ=255)
                                                  gTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCg                                                                                    >  2:236208/1‑136 (MQ=40)
                                                     gACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTg                                                                                 <  1:236208/136‑1 (MQ=38)
                                                     gACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCACGGCCAAACCATTACCGTTGa                                                                                <  1:424918/136‑1 (MQ=255)
                                                          cacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCACGGCCAAACCATTACCGTTGAAGCCg                                                                           <  1:15032/136‑1 (MQ=21)
                                                          cacTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCACGGCCAAACCATTACCGTTGAAGCCg                                                                           <  1:200413/136‑1 (MQ=21)
                                                           acaaCCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                           <  1:22645/133‑1 (MQ=17)
                                                           acTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCg                                                                           <  1:398810/136‑1 (MQ=34)
                                                                cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCaag                                                                      >  1:148935/1‑135 (MQ=14)
                                                                cgcgACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCaaa                                                                      >  2:5347/1‑136 (MQ=12)
                                                                       cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAggg                                                              >  1:221837/1‑136 (MQ=21)
                                                                                   ccGTCGGTCACCGAAGGCGGCGTGATCACATACCACGTGACCCTGAGCAACCCAGCCCAGACGCCGGTGACCGTGACCCTGTCAAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCCGGGCAGCGTCGacg                                                   >  1:379043/1‑134 (MQ=11)
                                                                                       cGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACg                                                                                                                            <  2:259797/59‑1 (MQ=35)
                                                                                                      gcttGATCACTTCCACCGTTAGCGTTATAATTCCTGCCCAGACGCCGGTGACATTGACCCTGTCAACGGCCAAACCATTACCGTTTAAGCCGGCAAAACCCAGTGCAGCGTCGATTTCCAGACCCCGGCGAATGat                               <  2:306038/133‑2 (MQ=11)
                                                                                                         tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGt                             <  1:193195/136‑1 (MQ=255)
                                                                                                         tGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGt                             <  1:154978/136‑1 (MQ=255)
                                                                                                          gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGCCCCCGGCGAATGACGtg                            >  1:87146/1‑135 (MQ=11)
                                                                                                                    acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACg                  >  2:195743/1‑136 (MQ=255)
                                                                                                                           cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGc                                                          >  1:27573/1‑89 (MQ=11)
                                                                                                                           cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGc                                                          <  2:27573/89‑1 (MQ=11)
                                                                                                                                   aaTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgtgg  >  2:232876/1‑135 (MQ=255)
                                                                                                                                   aaTCCTGCCCAGACGACGGTGACAGAGACCCTGACCACGGCCAAACCATTAGCGTTGAAGCCGGCATAACCCCGGGCAGCGTCGATTTCCAGACGCCGGCGAATGACGTCTCCAACAACGGTGGGTCGGTCAgtgg  >  2:483139/1‑135 (MQ=11)
                                                                                                                                     |                                                                                                                                      
AGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGACACCACTACCGCGACCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGT  >  NC_002947/196173‑196440

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