Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,295,820 C→G 8.5% intergenic (+66/+67) 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,8200CG 8.5% 145.2 / 4.8 47intergenic (+66/+67)PP_5496/PP_2023hypothetical protein/glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base G (2/2);  total (23/24)
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.39e-01

ACCCACAGCTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCG  >  NC_002947/2295677‑2295968
                                                                                                                                                |                                                                                                                                                    
acccacAGCTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACcg                                                                                                                                                 >  3:282833/1‑149 (MQ=255)
   cacaGCTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCgg                                                                                                                                              >  5:243280/1‑149 (MQ=255)
    acaGCTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGt                                                                                                                                             <  6:243280/149‑1 (MQ=255)
        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                                                                                                                                 >  8:166011‑M1/80‑132 (MQ=255)
                                  gccgagcggttccgccaggaaggcagccagctgtaccacaaagtctaaaccgcagcaccttaccctgtgggagcgggctTGCCCGCGAACACCGGCAGAGCCGGTGCCATGCACCGCAGTGCTTGCTTcgcg                                                                                                                                 <  7:166011‑M1/53‑1 (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                                                                                     <  5:195438/148‑1 (MQ=255)
                                                             gCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTa                                                                                     >  4:349038/1‑148 (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                                                >  6:394847/1‑149 (MQ=255)
                                                                                                 ttGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAagc                                                >  1:110264/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                                                           <  8:5641‑M1/121‑1 (MQ=255)
                                                                                                                   CCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAAc                                                           >  7:5641‑M1/1‑121 (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)
                                                                                                                                                |                                                                                                                                                    
ACCCACAGCTACACCGAGCAGTACCTGCGCCAGAGCGCCAATTTCCATGCGGCCCTGAGTGGCAAGCACCACCACTGATCTGGCTTGTGTGAGGTCCTTGTGGGAGCGGGCATGCCCGCGAACACCGGCAGAGCCGGTGCCATCCACCGCGGTGCCCGCTTCGCGGGCATGCCCGCTCCCACACGGACCCATGTGCAGCCCGATGTTTAATCAGCGGACCATGATTGCCTTGAACAGCTCCGAAGCCAGCCACCCCTCCAGCCAGCTCCGCACACGCGGATAAGCCGCCTCG  >  NC_002947/2295677‑2295968

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: