Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,931,581 A→T 12.2% F855Y (TTC→TAC) ‡ PP_2566 ← leucine‑rich repeat‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,931,5810AT12.2% 108.4 / 7.5 41F855Y (TTC→TAC) ‡PP_2566leucine‑rich repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (19/17);  new base T (3/2);  total (22/19)
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.48e-01

GTTGATGGGCTGTTGGCCAAGCAGTGAGGCGCAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTCCCGTTCGGCGTTGCTGATGTTGGCAAGCAACTCTTCAAGCGCGCTGTCC  >  NC_002947/2931433‑2931727
                                                                                                                                                    |                                                                                                                                                  
gTTGATGGGCTGTTGGCCAAGCAGTGAGGCGCAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGaa                                                                                                                                                   <  1:558733/150‑1 (MQ=255)
            ttGGCCAAGCAGTGAGGCGCAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCgaggag                                                                                                                                       <  2:116951/150‑1 (MQ=255)
             tGGCCAAGCAGTGAGGCGCAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGaggagg                                                                                                                                      <  1:544435/150‑1 (MQ=255)
             gggCGAACCAGGGATGCGCAGGGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGaggagg                                                                                                                                      <  1:393232/149‑1 (MQ=255)
                ccAAGCAGTGAGGCGCAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGTTCAGGcccagc                                                                                                                                         >  2:287996/1‑140 (MQ=255)
                ccAAGCAGTGAGGCGCAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGTTCAGGCCcagc                                                                                                                                         <  1:287996/144‑2 (MQ=255)
                  aaGCAGTGAGGCGCAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACg                                                                                                                                 >  1:845878/1‑150 (MQ=255)
                          aGGCGCAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCTGAGGAGGGCACGATCAAGcc                                                                                                                         <  1:537978/150‑1 (MQ=255)
                            gcgcAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCcgc                                                                                                                       >  2:442464/1‑150 (MQ=255)
                                     tGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCt                                                                                                              >  2:423095/1‑150 (MQ=255)
                                       gCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTTGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTg                                                                                                            >  2:154731/1‑150 (MQ=255)
                                            cACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGaa                                                                                                       <  2:99724/150‑1 (MQ=255)
                                             acgtAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGAGCGACGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAg                                                                                                      >  1:450598/4‑150 (MQ=255)
                                               cTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGc                                                                                                    >  1:581406/1‑150 (MQ=255)
                                                   ggCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGTTCAGGccca                                                                                                                                           <  1:384308/107‑3 (MQ=255)
                                                   ggCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGTTCAGGccca                                                                                                                                           >  2:384308/1‑105 (MQ=255)
                                                   ggCGAGTTGCAGCAACCGAGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCATCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTAGGCCTGAACAGTCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGc                                                                                                >  2:563488/1‑150 (MQ=255)
                                                    gCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCt                                                                                               >  1:351520/1‑150 (MQ=255)
                                                       aGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTcc                                                                                            <  2:581406/150‑1 (MQ=255)
                                                         ttGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTc                                                                                          <  1:727090/150‑1 (MQ=255)
                                                            gctgaaCCGGCCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGACGTGCTTCCTCGGc                                                                                       >  1:441878/5‑150 (MQ=255)
                                                              gcaACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGcccat                                                                                                                                          >  1:448522/1‑94 (MQ=255)
                                                              gcaACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGcccat                                                                                                                                          <  2:448522/97‑4 (MQ=255)
                                                                     ggCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAg                                                                              <  2:434844/150‑1 (MQ=255)
                                                                      gCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGc                                                                             >  2:451641/1‑150 (MQ=255)
                                                                           ttCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCAc                                                                        <  1:423095/150‑1 (MQ=255)
                                                                              aGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCg                                                                     >  1:818912/1‑150 (MQ=255)
                                                                                 cTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTcc                                                                  <  2:845878/150‑1 (MQ=255)
                                                                                     gggTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGttt                                                              <  2:743426/150‑1 (MQ=255)
                                                                                             ggCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACg                                                      >  2:352518/1‑150 (MQ=255)
                                                                                                ttGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTc                                                   <  2:415139/150‑1 (MQ=255)
                                                                                                   cAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCgg                                                                                        <  2:170454/110‑1 (MQ=255)
                                                                                                   cAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCgg                                                                                        >  1:170454/1‑110 (MQ=255)
                                                                                                   cAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGCCGGTTTGTTAGAAGTTcccc                                                >  2:49872/1‑149 (MQ=255)
                                                                                                         cTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTCCCGTTCGGc                                          <  1:49872/150‑1 (MQ=255)
                                                                                                            cTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTTCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTTAGACGTTCCCGTTAGGCGtt                                       >  2:405757/1‑150 (MQ=255)
                                                                                                                   cGTAGACTGTCGTCGTTGGCGAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTCCCGTTCGGCGTTGCTGATg                                >  2:435310/1‑150 (MQ=255)
                                                                                                                      aGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTCCCGTTCGGCGTTGCTGATGTTg                             <  2:150717/150‑1 (MQ=255)
                                                                                                                                    ggCCATGTTGGGCCTGAACAGTCCGAGGATGCCACGATACAACCGCCCTCGACCTTGCAGAAGGCGTGCTTCATTACCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTCCCGTTCGGCGTTGCTGATGTTGGCAAGCAACTCTTc               <  1:405757/150‑1 (MQ=255)
                                                                                                                                       cctcctcggCCTGTTCAGGCCCAGGAGGGCACGATCAAGCCGTGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTCCCGTTCGGCGTTGCTGATGTTGGCAAGCAACTCTTCAAg            >  2:54831/8‑150 (MQ=255)
                                                                                                                                       cAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTCCCGTTCGGCGTTGCTGATGTTGGCAAGCAACTCTTCAAg            >  1:533518/1‑150 (MQ=255)
                                                                                                                                                 cTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTCCCGTTCGGCGTTGCTGATGTTGGCAAGCAACTCTTCAAGCGCGCTGTcc  >  1:495668/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GTTGATGGGCTGTTGGCCAAGCAGTGAGGCGCAGTGCTGGCGATCACCTAGGGCGAGTTGCAGCAACCGGGCCGGTTCAGCCTGCGGGTGCGCGGCTTGCAGCGCCTGCTTGATGCGTAGACTGTCGTCGTTGGCCAGGCTGGGCCTGAACAGGCCGAGGAGGGCACGATCAAGCCGCGCTCGAGCTTGCAGAAGGCGTGCTTCCTCGGCAATGCGAAGCGGCACTCGTCCGTTTGTCAGACGTTCCCGTTCGGCGTTGCTGATGTTGGCAAGCAACTCTTCAAGCGCGCTGTCC  >  NC_002947/2931433‑2931727

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