Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 468,464 A→G 9.6% C831C (TGT→TGC)  PP_0386 ← sensory box protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947468,4640AG 9.6% 137.0 / 5.5 53C831C (TGT→TGC) PP_0386sensory box protein
Reads supporting (aligned to +/- strand):  ref base A (24/23);  new base G (3/2);  total (27/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

ACGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGTGATCT  >  NC_002947/468321‑468612
                                                                                                                                                |                                                                                                                                                    
aCGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGt                                                                                                                                                 <  5:116145/149‑1 (MQ=255)
aCGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGt                                                                                                                                                 <  3:19774/149‑1 (MQ=255)
   tCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                               <  1:42366/149‑1 (MQ=255)
   tCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGTTTGGACACCTGCGCACAGGTGAg                                                                                                                                              >  8:111274/1‑149 (MQ=255)
           cTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGCCCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCt                                                                                                                                      >  7:240260/1‑149 (MQ=255)
              ctctTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTg                                                                                                                                     >  6:225553/1‑147 (MQ=255)
                ctTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGc                                                                                                                                 <  6:156786/149‑1 (MQ=255)
                      ggCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGttt                                                                                                                           >  2:316857/1‑149 (MQ=255)
                        cgcTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                         <  1:49943/149‑1 (MQ=255)
                             ggCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACa                                                                                                                    <  1:316857/149‑1 (MQ=255)
                                 aTGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGt                                                                                                                <  3:276804/149‑1 (MQ=255)
                                     ccTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGa                                                                                                            >  7:61286/1‑149 (MQ=255)
                                       ttGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGt                                                                                                                                                 >  4:104075/1‑110 (MQ=255)
                                       ttGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGt                                                                                                                                                 <  3:104075/110‑1 (MQ=255)
                                            cagcagGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGAACAGGTGAGTACGCCCTGGAGCCGCTTTTCGAACAAGGTGCGATTggt                                                                                                       >  8:200475/1‑146 (MQ=255)
                                              gcagGTCGCAGACCGGGTGGCCAAGGCTGTCGTTTATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCTCTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAgc                                                                                                   <  7:200475/149‑2 (MQ=255)
                                                        gACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCg                                                                                         >  2:77574/1‑149 (MQ=255)
                                                         aCCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGt                                                                                                                <  7:197137/125‑1 (MQ=255)
                                                         aCCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGt                                                                                                                >  8:197137/1‑125 (MQ=255)
                                                                  ccAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCgg                                                                                        >  2:38517/1‑140 (MQ=255)
                                                                  ccAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCgg                                                                                        <  1:38517/140‑1 (MQ=255)
                                                                           tCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGa                                                                        <  4:212797/147‑1 (MQ=255)
                                                                           tCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATa                                                                      <  7:238560/149‑1 (MQ=255)
                                                                              ttGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCGCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTc                                                                   <  5:225553/149‑1 (MQ=255)
                                                                              ttGATGTGCGTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTc                                                                   <  8:240260/149‑1 (MQ=255)
                                                                                gATGTGCTTGAAGCGGTCGAGGTCCAGGAACCGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTc                                                                   >  5:81012/1‑147 (MQ=255)
                                                                                gATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCg                                                                  <  7:39203/148‑1 (MQ=255)
                                                                                      cTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTgg                                                           >  6:3646/1‑149 (MQ=255)
                                                                                         gAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCt                                                        >  5:85515/1‑149 (MQ=255)
                                                                                         gAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCt                                                        >  2:139220/1‑149 (MQ=255)
                                                                                           aGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTgg                                                      >  8:68095/1‑149 (MQ=255)
                                                                                                tCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTg                                                 >  1:254480/1‑149 (MQ=255)
                                                                                                     gTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGa                                            <  2:254480/149‑1 (MQ=255)
                                                                                                       ccAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATg                                          >  6:310682/1‑149 (MQ=255)
                                                                                                        cAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTc                                                                                     >  1:200741/1‑105 (MQ=255)
                                                                                                        cAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTc                                                                                     <  2:200741/105‑1 (MQ=255)
                                                                                                                          ccTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCaaaa                        <  1:77574/149‑1 (MQ=255)
                                                                                                                              gCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCt                              <  7:37590/139‑1 (MQ=255)
                                                                                                                              gCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCt                              >  8:37590/1‑139 (MQ=255)
                                                                                                                              gCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAAcc                     <  6:85515/148‑1 (MQ=255)
                                                                                                                                ggTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGga                                     >  3:221454/1‑130 (MQ=255)
                                                                                                                                ggTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGga                                     <  4:221454/130‑1 (MQ=255)
                                                                                                                                ggTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGcc                  >  5:209104/1‑149 (MQ=255)
                                                                                                                                   tactcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGga                                     >  6:263929/5‑127 (MQ=255)
                                                                                                                                   tactcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGga                                     <  5:263929/123‑1 (MQ=255)
                                                                                                                                    actcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTgg                                                      >  2:96552/4‑109 (MQ=255)
                                                                                                                                    actcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTgg                                                      <  1:96552/106‑1 (MQ=255)
                                                                                                                                      acacCTGCGCCCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGGCGAGTTTGGCCTGGGCGTGCtt                                              <  1:265612/115‑1 (MQ=255)
                                                                                                                                      acacCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGGGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCtg                                              >  2:265612/1‑114 (MQ=255)
                                                                                                                                        acCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTg          >  8:41930/1‑149 (MQ=255)
                                                                                                                                        acCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTg          >  8:117169/1‑149 (MQ=255)
                                                                                                                                         gtatctcacaCGTGTCTACGCCCTGGATCCGTTTTTCGAACTATGTGCGATTGGGCAGTCCGGTCAGCGGGTCTTGGTGGGATTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTgg         <  4:166845/143‑1 (MQ=255)
                                                                                                                                          cTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTggg        >  6:1489/1‑149 (MQ=255)
                                                                                                                                          cTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGGCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTggg        <  8:61286/149‑1 (MQ=255)
                                                                                                                                                acaGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCACCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCACCACGAAGTgggacatcg  >  3:166845/1‑143 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ACGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTG‑CGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGTGATCT  >  NC_002947/468321‑468612

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