Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,295,833 C→T 8.3% intergenic (+79/+54) PP_5496 → / ← PP_2023 hypothetical protein/glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,295,8330CT 8.3% 133.0 / 2.5 48intergenic (+79/+54)PP_5496/PP_2023hypothetical protein/glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base C (21/23);  new base T (2/2);  total (23/25)
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.85e-01

CTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCGGCGAACCACTGCG  >  NC_002947/2295685‑2295981
                                                                                                                                                     |                                                                                                                                                    
cTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCGTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGccc                                                                                                                                                      <  4:269180/149‑1 (MQ=255)
      cGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTcg                                                                                                                                                <  4:282833/149‑1 (MQ=255)
            gTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCa                                                                                                                                          >  3:107985/1‑149 (MQ=255)
                       agagCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGGGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTccc                                                                                                                               <  7:90466/149‑1 (MQ=255)
                         agCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCcac                                                                                                                             <  2:353707/149‑1 (MQ=255)
                          gccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttaccctgtgggagcgggctTGCCCGCGAACACCGGCAGAGCCGGTGCCATGCACCGCAGTGCTTGCTTcgcg                                                                                                                                              <  7:166011‑M1/53‑1 (MQ=255)
                          gccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttaccctgtgggagcgggctTGCCCGCGAACACCGGCAGAGCCGGTGCCATGCACCGCAGTGCTTGCTTcgcg                                                                                                                                              >  8:166011‑M1/80‑132 (MQ=255)
                           cgcCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCacac                                                                                                                           <  2:307554/149‑1 (MQ=255)
                           cgcCAAATTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCacac                                                                                                                           <  8:172319/149‑1 (MQ=255)
                            gcCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACg                                                                                                                          >  5:292263/1‑149 (MQ=255)
                            cgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcacctcaccctgtgggagcgggctTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                                                             <  6:124007‑M1/54‑1 (MQ=255)
                            cgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcacctcaccctgtgggagcgggctTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCgg                                                                                                                                             >  5:124007‑M1/78‑131 (MQ=255)
                                      aTGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGATGGGGCCGGGCATGCCCCGCACCCCCGGCAAGGCCGGGCCCTACCCCCCCGTTCCCCGCTTGCCGGGCATCCCCCGGCCCCCCCGGCCCCAtgt                                                                                                                 >  1:107633/1‑149 (MQ=255)
                                      aTGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCAtgtg                                                                                                                >  1:306082/1‑149 (MQ=255)
                                       tGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGGATGCCCACTCCCACAAGGACCCCTGTGc                                                                                                               >  8:372759/1‑149 (MQ=255)
                                          ggCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCAGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGAc                                                                                                                       <  8:301658/138‑1 (MQ=255)
                                          ggCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCAGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGAc                                                                                                                       >  7:301658/1‑138 (MQ=255)
                                            cccTCAGGGGCAAGCACCACCACTGTGCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCTAACACCGGCAGCGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGccc                                                                                                          <  2:370195/149‑1 (MQ=37)
                                                gAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCGTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGAt                                                                                                       <  6:292263/148‑1 (MQ=255)
                                                 aGTGGCAAGCACCACCAGTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGTTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGAt                                                                                                       <  7:372759/147‑1 (MQ=255)
                                                 aGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGt                                                                                                     >  1:265008/1‑149 (MQ=255)
                                                  gTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGtt                                                                                                    <  2:265008/149‑1 (MQ=255)
                                                     gCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTa                                                                                                  >  4:349038/1‑148 (MQ=255)
                                                     gCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTa                                                                                                  <  5:195438/148‑1 (MQ=255)
                                                     gCAAGAACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTaa                                                                                                 <  8:4799/149‑1 (MQ=255)
                                                         gCACCCCCATTTTTCTGGCTTGTGTGTGGGCCTTGTGGGTGCGGGCTTGCCCGCTAAAACCGGCATAGCGGTTGCAATCCACCGCGGTGCCCGCTTCGCGGGAATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAg                                                                                             <  2:157203/149‑1 (MQ=255)
                                                         gCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGtt                                                                                                    >  7:2285/1‑142 (MQ=255)
                                                             caccacTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCCTGTGCAGCCCGCTGTGTAATCAGCggc                                                                                         >  4:258318/1‑148 (MQ=255)
                                                              accacTGATCTGGGTTGGGGGAGTGCCTTGGGGGCGCGGGCATCCCCGCTAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGAc                                                                                        <  3:196464/149‑1 (MQ=37)
                                                                     aTCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCCCATAGCCGGCCTGCCCGCGCCCCCACGGACCCCAGTGCAGCCCCATGTTTAATCAGCGGACCATGAtt                                                                                 >  6:38934/1‑149 (MQ=37)
                                                                         gtcTTGTTGGAGGTCCTTGTGGGAGCGGGCATGCCCGCTAACACCGGCAGAGCCGGTGCAATCCACCGCTGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGcc                                                                              <  7:149213/146‑1 (MQ=37)
                                                                            ttGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTg                                                                           >  1:6226/1‑148 (MQ=255)
                                                                                     gTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCCCCCGGACCCATGTGCAGCCCGATGTTTAATAAGCCGACCATGATTTCCTTGAACAGCTCCg                                                                 >  4:254909/1‑149 (MQ=37)
                                                                                         ttGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAagc                                                             >  1:110264/1‑149 (MQ=255)
                                                                                         ttGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAagc                                                             >  6:394847/1‑149 (MQ=255)
                                                                                           tgggggCGCGGGCATGCCCGCTAACACCGGCAGAGCGGGTGCCATCCACGGGGGTGCCCGCTTGGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAagc                                                             <  5:38934/145‑1 (MQ=38)
                                                                                             tgGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAagcca                                                           <  4:107985/148‑1 (MQ=255)
                                                                                              gggAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCAACCACCCCGACCCCTGTGCAGCCCCATGTTTAATAAGCGGACCATTATTTCCTTGAACACCAACGAagccagc                                                         >  4:233906/1‑149 (MQ=25)
                                                                                                           CCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAAc                                                                        >  7:5641‑M1/1‑121 (MQ=255)
                                                                                                           CCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAAc                                                                        <  8:5641‑M1/121‑1 (MQ=255)
                                                                                                                    acCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGca                                   <  3:233906/149‑1 (MQ=255)
                                                                                                                      cGGCAGAGCCAGTGCCATCCACCGCGGTGTCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGcaca                                 >  6:345442/1‑149 (MQ=255)
                                                                                                                         cAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACAcgc                              <  2:6226/149‑1 (MQ=255)
                                                                                                                                     cATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAgccgc                   >  1:229474/1‑148 (MQ=255)
                                                                                                                                        ccACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCg               >  5:147490/1‑149 (MQ=255)
                                                                                                                                          aCCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCg                                                                                           >  6:123085/1‑71 (MQ=255)
                                                                                                                                          aCCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCg                                                                                           <  5:123085/71‑1 (MQ=255)
                                                                                                                                             gcgGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGAGCATGATTTCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCAGACGCTGATAAGCCGCCTCGGCGaa          <  3:349038/149‑1 (MQ=255)
                                                                                                                                              cgGGGCCCGCTTCGCGGGCATGCCCGCGCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCGGCGAAc         <  3:258318/149‑1 (MQ=255)
                                                                                                                                                     cGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCGGCGAACCACTGCg  >  7:324603/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCGGCGAACCACTGCG  >  NC_002947/2295685‑2295981

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