Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,575,089 G→A 9.3% G75S (GGC→AGC)  PP_2258 → sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,575,0890GA 9.3% 112.8 / 2.9 43G75S (GGC→AGC) PP_2258sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base A (2/2);  total (21/22)
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.27e-01

TACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAGCGCCAGCTGCTGCATACCTTG  >  NC_002947/2574942‑2575238
                                                                                                                                                   |                                                                                                                                                     
tACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGc                                                                                                                                                     >  1:80185/1‑150 (MQ=255)
   ggcAATCAGCAGTTCTGCAGCATCTCGGGCTACAGACGTGAGGCACTGCTGGGGGCCAACCCCCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGaa                                                                                                                                                  >  2:462178/3‑150 (MQ=255)
   gTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGaa                                                                                                                                                  <  2:130746/150‑1 (MQ=255)
     cAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAAt                                                                                                                                                <  2:351054/150‑1 (MQ=255)
      aaTCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATc                                                                                                                                               <  2:885800/150‑1 (MQ=255)
       aTCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCt                                                                                                                                              >  1:769451/1‑150 (MQ=255)
             agTTATGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAAcc                                                                                                                                        <  2:519500/150‑1 (MQ=255)
                   gcagcaTCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGAGCGAAATCt                                                                                                                                              >  1:964451/1‑138 (MQ=255)
                   gcagcaTCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGAGCGAAATCt                                                                                                                                              >  1:964444/1‑138 (MQ=255)
                   gcagcaTCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGAGCGAAATCt                                                                                                                                              <  2:964444/138‑1 (MQ=255)
                   gcagcaTCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGAGCGAAATCt                                                                                                                                              <  2:964451/138‑1 (MQ=255)
                    cagcaTCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCCCCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCaa                                                                                                                                 <  2:477111/150‑1 (MQ=255)
                         tctcGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACg                                                                                                                            >  1:8665/1‑150 (MQ=255)
                                             gAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGAc                                                                                                        <  1:794894/150‑1 (MQ=255)
                                                  gctgGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCAc                                                                                                   <  1:241052/150‑1 (MQ=255)
                                                    tgGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCa                                                                                                 >  2:864813/1‑150 (MQ=255)
                                                        ggCCAACCACCGCATCCTCAATTCCGGGCTGCGTGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGt                                                                                             <  2:279228/150‑1 (MQ=255)
                                                                accGCAGCCTCAATTCCGGGCGGCATGAACCGACGTTTGTTGTCGGCAGGTGGCGGGCTTTGGGCGCCGGGAGGGGGTAGAATGGCGGAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGTGTCGACAGCACCATGGTGCAACtcg                                                                                     >  2:315611/1‑148 (MQ=255)
                                                                      tCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATc                                                                               <  2:769451/150‑1 (MQ=255)
                                                                         tCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATg                                                                                               >  1:270915/1‑131 (MQ=255)
                                                                         tCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATg                                                                                               <  2:270915/131‑1 (MQ=255)
                                                                             ttCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGAc                                                                        >  1:241081/1‑150 (MQ=255)
                                                                                 gggCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCAACTGGTCGATCCGCTGACAGGc                                                                    <  2:801103/150‑1 (MQ=255)
                                                                                  ggCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAACTCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGcc                                                                   >  1:250327/1‑150 (MQ=255)
                                                                                             ccGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAg                                                        >  2:523178/1‑150 (MQ=255)
                                                                                                aCGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAg                                                                      <  2:492393/133‑1 (MQ=255)
                                                                                                aCGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAg                                                                      >  1:492393/1‑133 (MQ=255)
                                                                                                aCGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTAt                                                     <  2:8665/150‑1 (MQ=255)
                                                                                                  gTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTAtgt                                                   <  1:523178/150‑1 (MQ=255)
                                                                                                          tCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATcc                                           >  1:66495/1‑150 (MQ=255)
                                                                                                                  tgGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGAc                                   <  2:66495/150‑1 (MQ=255)
                                                                                                                    gCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTAtgtg                                                  >  2:321754/1‑133 (MQ=255)
                                                                                                                    gCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTAtgtg                                                  <  1:321754/133‑1 (MQ=255)
                                                                                                                               tCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATgg                                                                                              >  1:908501/1‑78 (MQ=255)
                                                                                                                               tCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATgg                                                                                              <  2:908501/78‑1 (MQ=255)
                                                                                                                               tCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAgc                      >  2:1857/1‑150 (MQ=255)
                                                                                                                               tCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAgc                      >  2:1589/1‑150 (MQ=255)
                                                                                                                                     ggcgggTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCt                                        >  1:95424/1‑126 (MQ=255)
                                                                                                                                     ggcgggTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCt                                        <  2:95424/126‑1 (MQ=255)
                                                                                                                                      gcgggTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAGCGCCAgct               <  2:80185/150‑1 (MQ=255)
                                                                                                                                          gtgtGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAGCGCCAGctgctg           >  1:500458/1‑150 (MQ=255)
                                                                                                                                               gAAGGGCGAAATCTGCAACCGGTTCAAGGACGGTTCGTTGTACTGGGTCGAGAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAGCGCCAGCTGCTGCATAc      <  2:70211/150‑1 (MQ=255)
                                                                                                                                                aaGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAGCGCCAGCTGCTGCATAcc     >  2:672208/1‑150 (MQ=255)
                                                                                                                                                   ggCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAGCGCCAGCTGCTGCATACCTTg  <  2:500458/150‑1 (MQ=255)
                                                                                                                                                   ggCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAGCGCCAGCTGCTGCATACCTTg  >  2:264292/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TACGTCAATCAGCAGTTCTGCAGCATCTCGGGCTACAGCCGTGAGGAACTGCTGGGGGCCAACCACCGCATCCTCAATTCCGGGCTGCATGAACCGACGTTTTTTGTCGGCATGTGGCGGGCTTTGGTCGCCGGGCGGGTGTGGAAGGGCGAAATCTGCAACCGGTCCAAGGACGGTTCGTTGTACTGGGTCGACAGCACCATGGTGCCACTGGTCGATCCGCTGACAGGCCAGGTGCGCAAGTATGTGTCGATCCGCTTTGACGTCACGGAGAAGCGCCAGCTGCTGCATACCTTG  >  NC_002947/2574942‑2575238

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