Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,908,447 C→T 11.8% L118L (CTG→TTG)  PP_3450 → TPR repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,908,4470CT11.8% 84.4 / 3.6 34L118L (CTG→TTG) PP_3450TPR repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (14/16);  new base T (2/2);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.05e-01

GCAGCCGCTGGCCGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCG  >  NC_002947/3908305‑3908590
                                                                                                                                              |                                                                                                                                               
gCAGCCGCTGGCCGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCtggccat                                                                                                                                          <  2:513137/150‑1 (MQ=255)
      gCTGGCCGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCaa                                                                                                                                    <  2:575575/150‑1 (MQ=255)
              acacAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGGCGCGGTTGGCCATGGCCAGCcgcgtct                                                                                                                            <  2:280191/150‑4 (MQ=255)
                         ttGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGc                                                                                                                 >  1:52074/1‑150 (MQ=255)
                              tgccCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGACCAACCGCGCCATTGCCGTCGGccacca                                                                                                            >  2:408610/1‑150 (MQ=255)
                                      aCTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCGGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAgg                                                                                                    <  2:183918/150‑1 (MQ=255)
                                           cAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCaa                                                                                               <  1:193332/150‑1 (MQ=255)
                                               gcatGGCCAACACCTAGATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACcct                                                                                           >  2:500292/4‑150 (MQ=255)
                                                      cAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCTCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccac                                                                                                           >  1:420132/1‑127 (MQ=255)
                                                      cAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGccaccac                                                                                                           <  2:420132/127‑1 (MQ=255)
                                                            cTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTTGCCAGGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCt                                                                              >  1:434248/1‑150 (MQ=255)
                                                            cTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                         <  1:376036/139‑1 (MQ=255)
                                                            cTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACgctgg                                                                                         >  2:376036/1‑139 (MQ=255)
                                                                   cAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCt                                                                       >  2:371363/1‑150 (MQ=255)
                                                                     gCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTg                                                                      >  1:519151/1‑149 (MQ=255)
                                                                     gCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTg                                                                      <  2:519151/149‑1 (MQ=255)
                                                                          tCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACG‑CCCCTCCCCCGCAGCGGCCAACCAGACGCGGTtggccatggccagc                                                                                                                                   <  1:356837/82‑3 (MQ=255)
                                                                          tCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACG‑CCCCTCCCCCGCAGCGGCCAACCAGACGCGGTtggccatggccagc                                                                                                                                   >  2:356837/1‑80 (MQ=255)
                                                                                         tGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGagc                                                 >  1:479815/1‑150 (MQ=255)
                                                                                             tGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAgcg                                             >  1:36096/1‑150 (MQ=255)
                                                                                                    ctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGc                                                                         <  1:4360/115‑1 (MQ=255)
                                                                                                    ctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGc                                                                         >  2:4360/1‑115 (MQ=255)
                                                                                                    ctgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCaa                                      <  1:336514/150‑1 (MQ=255)
                                                                                                     tgGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAAc                                     <  1:178160/150‑1 (MQ=255)
                                                                                                       gCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACtt                                   <  2:36096/150‑1 (MQ=255)
                                                                                                        cccACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACttt                                  <  1:371363/150‑1 (MQ=255)
                                                                                                               ccccTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACt                           >  1:121044/1‑150 (MQ=255)
                                                                                                                cccTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTg                          <  2:15317/150‑1 (MQ=255)
                                                                                                                            gCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGc              <  1:227089/150‑1 (MQ=255)
                                                                                                                                       gACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                 >  2:311586/1‑120 (MQ=255)
                                                                                                                                       gACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTg                                 <  1:311586/120‑1 (MQ=255)
                                                                                                                                        gcgcgGTTGGCCATGGCCAGCCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAgcg  >  2:115165/2‑150 (MQ=255)
                                                                                                                                            ggCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTgc                                                             >  2:329518/1‑87 (MQ=255)
                                                                                                                                            ggCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTgc                                                             <  1:329518/87‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
GCAGCCGCTGGCCGACACAGGCCAGTTGCCTGCCCAGGACTGGCAACGCCTGGCCAATGCCTATATCCAGCAAGTCGACGAAAAGGGCCTGACTGCGCTGCTGGCCCACGCCCCCTCCCCCGCAGCGGCCAACCAGACGCGGCTGGCCATGGCCAACCGCGCCATTGCCGTCGGCCACCACCAGCAGGCCCAACGCTGGCTGCAGTCGCTACCGCCTGATCAACTGCAGCAACCCGAGCAGCGCCAGCAACTTTGGCAACTGGCCCGCGAAGGCGACGATGCAGCG  >  NC_002947/3908305‑3908590

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: ‑