Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 783,248 G→A 7.2% A894V (GCC→GTC)  PP_0672 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947783,2480GA 7.2% 157.8 / 2.1 55A894V (GCC→GTC) PP_0672sensory box protein
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base A (2/2);  total (26/29)
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.97e-01

TGGATGGCGCGGTTCAGGGCCATTTCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGT  >  NC_002947/783101‑783390
                                                                                                                                                    |                                                                                                                                              
tGGATGGCGCGGTTCCGGCCCATCTCTCGCGTGGGTGGTCGTTGCAGCAACAGGGTAAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGCCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCgccggcc                                                                                                                                              <  2:9822/150‑1 (MQ=255)
       cgcgGTTCAGGGCCATTTCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGc                                                                                                                                       <  1:364339/150‑1 (MQ=255)
        gcgGTTCAGGGCCATTTCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCCGGAGAGACTACATTCTCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGcc                                                                                                                                      >  2:870302/1‑150 (MQ=255)
            ttCAGGGCCATTTCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAgg                                                                                                                                  <  1:956396/150‑1 (MQ=255)
            ttCAGGGCCATTTCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAgg                                                                                                                                  <  2:425901/150‑1 (MQ=255)
               aGGGCCATTTCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAg                                                                                                                               >  1:378266/1‑150 (MQ=255)
                  gCCATTTCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCg                                                                                                                            <  2:632176/150‑1 (MQ=255)
                   ccATTTCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCCGGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGt                                                                                                                           <  2:31009/150‑1 (MQ=255)
                      tttCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGCCGGCGACCGCATgctcaa                                                                                                                                    <  2:495433/138‑6 (MQ=255)
                      tttCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGCCGGCGACCGCATgctcaa                                                                                                                                    >  1:495433/1‑133 (MQ=255)
                         cACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCg                                                                                                                     <  1:429806/150‑1 (MQ=255)
                          aCGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCgg                                                                                                                    >  2:1107630/1‑150 (MQ=255)
                          aCGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGCCGGCGACCGCATgctcaag                                                                                                                                   >  2:333018/1‑129 (MQ=255)
                          aCGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGCCGGCGACCGCATgctcaag                                                                                                                                   <  1:333018/135‑7 (MQ=255)
                             ggTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGttt                                                                                                                 <  1:302044/150‑1 (MQ=255)
                             ggTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGttt                                                                                                                 <  1:639130/150‑1 (MQ=255)
                                ggggCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCGGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGaa                                                                                                              <  1:785179/150‑1 (MQ=255)
                                  ggCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCCTGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGc                                                                                                            <  1:197808/150‑1 (MQ=255)
                                   gcgcGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCg                                                                                                           <  1:36301/150‑1 (MQ=255)
                                    cgcgGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCgg                                                                                                          <  1:1090229/150‑1 (MQ=255)
                                            cagcaACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGt                                                                                                                   <  1:139569/133‑1 (MQ=255)
                                            cagcaACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGt                                                                                                                   >  2:139569/1‑133 (MQ=255)
                                             agcaACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCATTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCggtcgaggtcg                                                                                                 >  1:407936/1‑150 (MQ=255)
                                               caACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCGGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGgtcgaggtcgag                                                                                               <  2:214034/150‑1 (MQ=255)
                                                aaCAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGtcgaggtcgagg                                                                                              >  1:82405/1‑150 (MQ=255)
                                                  cAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGaa                                                                                            <  2:378266/150‑1 (MQ=255)
                                                     ggTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACAt                                                                                         >  2:32771/1‑150 (MQ=255)
                                                            tcgtcgCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGcacc                                                                                  >  1:494934/1‑150 (MQ=255)
                                                               tcgCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCaccacc                                                                               >  2:921077/1‑150 (MQ=255)
                                                               tcgCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCaccacc                                                                               >  2:188061/1‑150 (MQ=255)
                                                                      ccATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAg                                                                                      >  1:500603/1‑136 (MQ=255)
                                                                      ccATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAg                                                                                      <  2:500603/136‑1 (MQ=255)
                                                                      acATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCt                                                                        >  2:1127375/2‑150 (MQ=255)
                                                                       cATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCtt                                                                       <  1:1127375/150‑1 (MQ=255)
                                                                       cATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCtt                                                                       <  1:197665/150‑1 (MQ=255)
                                                                        acgCGTGCCACCGTGGCGTCGTAATCCAAGCAGGCCAGCAGGCGAAGGGCCATGTCCTTGAGCATGCGGTCGCAGGCCGCATGGCCGAGGGAGTCGTTGATAGGTTTGAAGCAGACGAGGTCGAGGAACATCAGCACCACCCACGCCTTc                                                                      >  2:228320/3‑150 (MQ=255)
                                                                                 aCCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCg                                                             >  1:1094447/1‑150 (MQ=255)
                                                                                  ccGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAAACACATCTCCACCCAGGCCTTCTGCCGCTCgg                                                            <  1:870302/150‑1 (MQ=255)
                                                                                     tgtCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCct                                                         >  2:775499/1‑150 (MQ=255)
                                                                                                  cACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGgtcgaggtcgag                                                                                               >  2:1057384/1‑99 (MQ=255)
                                                                                                  cACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGgtcgaggtcgag                                                                                               <  1:1057384/99‑1 (MQ=255)
                                                                                                  cACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGtt                                            >  2:1084405/1‑150 (MQ=255)
                                                                                                               cagGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTGTGAAGCGGTCGAGGTCGAGGAACATCATCACCACCCAGGCCTTCTGCCGCTCAGCCTGCTGCAGGGCGTTGTTCAAActttt                                <  2:515551/150‑5 (MQ=255)
                                                                                                                  gcgcAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTg                             >  1:247048/1‑150 (MQ=255)
                                                                                                                       gggCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACa                        <  1:1107630/150‑1 (MQ=255)
                                                                                                                          ccATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCggtcgaggtcga                                                                                                >  1:133425/1‑75 (MQ=255)
                                                                                                                          ccATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCggtcgaggtcga                                                                                                <  2:133425/75‑1 (MQ=255)
                                                                                                                            aTGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGtgtg                   >  1:792857/1‑150 (MQ=255)
                                                                                                                            aTGTCCTTGAGCATGCGGGCGCCGGCCGCAGGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCGCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGtgtg                   <  2:80820/150‑1 (MQ=255)
                                                                                                                                 cTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGtt              >  2:1128752/1‑150 (MQ=255)
                                                                                                                                  ttGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTg             >  1:93344/1‑150 (MQ=255)
                                                                                                                                        aTGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAAc                                                                                           <  1:916911/66‑1 (MQ=255)
                                                                                                                                        aTGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAAc                                                                                           >  2:916911/1‑66 (MQ=255)
                                                                                                                                        aTGCGGTCGCCAGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAgg       >  1:425477/1‑150 (MQ=255)
                                                                                                                                         tGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGt      <  2:247048/150‑1 (MQ=255)
                                                                                                                                             gTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGt  <  1:189488/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TGGATGGCGCGGTTCAGGGCCATTTCACGGGTGGGGCGCGGTTGCAGCAACAGGGTGAACTCGTCGCCGCCCATGCGTGCCACCGTGTCGTCGTCATCCACGCAGGCCAGCAGGCGCAGGGCCATGTCCTTGAGCATGCGGTCGCCGGCCGCATGGCCGAGGGAGTCGTTGATCGGTTTGAAGCGGTCGAGGTCGAGGAACATCAGCACCACCCAGGCCTTCTGCCGCTCGGCCTGCTGCAGGGCGTTGTACAAACGGTCCTGGAACAGTGTGCGGTTGGGCAGGTGGGT  >  NC_002947/783101‑783390

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