Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 782,158 C→A 7.4% P1257P (CCG→CCT)  PP_0672 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947782,1580CA 7.4% 165.9 / 3.1 54P1257P (CCG→CCT) PP_0672sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (29/21);  new base A (2/2);  total (31/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.83e-01

GCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGCCCGGGCAATCTGCGCGT  >  NC_002947/782016‑782303
                                                                                                                                              |                                                                                                                                                 
gCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGgccgcc                                                                                                                                            >  1:680685/1‑150 (MQ=255)
gCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGgccgcc                                                                                                                                            >  1:283393/1‑150 (MQ=255)
  aGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCg                                                                                                                                          >  2:901400/1‑150 (MQ=255)
    gAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCGTCGGCCGCCCGaa                                                                                                                                        >  2:1255293/1‑150 (MQ=255)
    gAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGaa                                                                                                                                        <  1:1238114/150‑1 (MQ=255)
    gAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGaa                                                                                                                                        <  2:680685/150‑1 (MQ=255)
      aTGCCGCGCGGTGGATGGCAGCGTGTTTACCGGTGTTCGGGGGTAGACCCGCGCCTGCGGGGGCCTCGTGTGCGGGTCAACGGAACATGACCCAGGTCTCGTTGGAGAATTGCGCCTCGAACTGTGTCTCTGGCATCGGGCGGCCGAACa                                                                                                                                      <  2:464864/150‑1 (MQ=255)
           gcgcGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAg                                                                                                                                 <  1:901400/150‑1 (MQ=255)
            cgcgGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGACCGAACAGATatc                                                                                                                                <  1:670790/150‑3 (MQ=255)
                         agcGGCGTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCtcgtc                                                                                                                   <  1:220355/148‑1 (MQ=255)
                         aCCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATAGGCCGCCCGAACAGATAGCCCTGCACCtcgtc                                                                                                                   >  2:496399/1‑150 (MQ=255)
                         aCCGGCAGTGCCGGTGGTCGCGGGTGAACCCGCGCCCGAGGGGGCCTTGCGTGGGGGTCAATGGAACATGAACAAAGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCtcgtc                                                                                                                   >  2:295654/1‑150 (MQ=255)
                           cGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGc                                                                                                                 >  2:283777/1‑150 (MQ=255)
                            ggCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCa                                                                                                                >  2:1225046/1‑150 (MQ=255)
                                         ttCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGGCAATGGAAAATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATAGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCgcgc                                                                                                   >  2:291541/1‑150 (MQ=255)
                                         agcgcgGGTTTACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCgcgc                                                                                                   >  2:882521/3‑150 (MQ=255)
                                          tCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAg                                                                                                                                 <  1:453950/119‑1 (MQ=255)
                                          tCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAg                                                                                                                                 >  2:453950/1‑119 (MQ=255)
                                          gcgcgGGTTTACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCa                                                                                                  >  1:879643/2‑150 (MQ=255)
                                                  aaaCCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTcc                                                                                          >  2:393308/1‑150 (MQ=255)
                                                   aaCCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATAGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCa                                                                                         <  1:496399/150‑1 (MQ=255)
                                                     cccGCGCCTACGGGGGCCTTGTGTCCGGGGCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGc                                                                                       >  1:174533/1‑150 (MQ=255)
                                                      ccGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCt                                                                                      >  1:1080821/1‑150 (MQ=255)
                                                            cTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCGTCGGCCGCCCGAACCGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGt                                                                                <  1:1255293/150‑1 (MQ=255)
                                                                      cTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCAc                                                                      <  2:1080821/150‑1 (MQ=255)
                                                                          tgtgCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCa                                                                                                  <  1:814063/118‑1 (MQ=255)
                                                                          tgtgCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCa                                                                                                  >  2:814063/1‑118 (MQ=255)
                                                                               gggTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTACCACGCCCTCGGc                                                             >  2:188343/1‑150 (MQ=255)
                                                                                 gTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGtt                                                                           >  1:751292/1‑134 (MQ=255)
                                                                                 gTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGtt                                                                           <  2:751292/134‑1 (MQ=255)
                                                                                  tCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGAt                                                          <  2:649858/150‑1 (MQ=255)
                                                                                     aTGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCATGCTGTTCATCGCAACCATGCTCGCGCAGGAACTCCAGCTGCTCGTGGGTTTCCACCCCCTCGGCGATGAc                                                       >  1:987508/1‑150 (MQ=255)
                                                                                         aaCATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGcc                                                   <  1:283777/150‑1 (MQ=255)
                                                                                                aCAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGa                                            >  1:370929/1‑150 (MQ=255)
                                                                                                    ggTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAg                                                 >  2:836565/1‑141 (MQ=255)
                                                                                                    ggTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAg                                                 >  2:836491/1‑141 (MQ=255)
                                                                                                    ggTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAg                                                 <  1:836565/141‑1 (MQ=255)
                                                                                                    ggTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAg                                                 <  1:836491/141‑1 (MQ=255)
                                                                                                       ctcGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCtgtg                                     >  1:818380/1‑150 (MQ=255)
                                                                                                        tcGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCtgtgt                                    <  1:295654/150‑2 (MQ=255)
                                                                                                        tcGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGACCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTgg                                    >  2:467619/1‑150 (MQ=255)
                                                                                                          gTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGc                                  >  2:1141133/1‑150 (MQ=255)
                                                                                                                  aTTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCga                          <  1:882521/150‑1 (MQ=255)
                                                                                                                          tCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGcc                  >  2:514212/1‑150 (MQ=255)
                                                                                                                          tCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGcc                  >  1:878065/1‑150 (MQ=255)
                                                                                                                            gAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGCCCg                >  2:711307/1‑150 (MQ=255)
                                                                                                                               cTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATCATCGccatggc             <  1:467619/150‑6 (MQ=255)
                                                                                                                                tGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGCCCGGGCa            >  1:234032/1‑150 (MQ=255)
                                                                                                                                  aTGCGCTGGCATAGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGCCCGGGCAAt          <  1:291541/150‑1 (MQ=255)
                                                                                                                                    gcgcTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGATGAGGCTGTGGGCCATGGCGATGATCGCCCGGGCAATCt        <  2:987508/150‑1 (MQ=255)
                                                                                                                                        tGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCTCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGCGGCTGTGGGCCATGGCGATGATCGCCCGGGCAATCTgcgc    <  1:1225046/150‑1 (MQ=255)
                                                                                                                                        tGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGCCCGGGCAATCTgcgc    >  2:265511/1‑150 (MQ=255)
                                                                                                                                        tGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGCCCGGGCAATCTgcga    <  2:879643/150‑2 (MQ=255)
                                                                                                                                          gCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGCCCGGGCAATCTGCGCGt  >  1:100642/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GCAGGAATGCGGCGCGGTGGATGGCACCGGCTTTGCCGGTGTTCGCGGGTAAACCCGCGCCTACGGGGGCCTTGTGTGCGGGTCAATGGAACATGAACAAGGTCTCGTTGGAGAATTGCGCCTCGAACTGATGCGCTGGCATCGGCCGCCCGAACAGATAGCCCTGCACCTCGTCGCAACCATGCTCGCGCAGGAACTCCAGCTGTTCGTGGGTTTCCACGCCCTCGGCGATGACTGCCAGGTTGAGGCTGTGGGCCATGGCGATGATCGCCCGGGCAATCTGCGCGT  >  NC_002947/782016‑782303

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