Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,808,192 G→C 11.8% intergenic (‑301/‑254) PP_4234 ← / → dsbD‑II bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,808,1920GC11.8% 88.9 / 5.2 34intergenic (‑301/‑254)PP_4234/dsbD‑IIbifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (13/17);  new base C (2/2);  total (15/19)
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.05e-01

ATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAA  >  NC_002947/4808047‑4808324
                                                                                                                                                 |                                                                                                                                    
gttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCgcgggcg                                                                                                                                   <  5:264535‑M1/22‑1 (MQ=255)
     tcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagacacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                              <  6:283210‑M1/27‑1 (MQ=255)
     tcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                              <  6:287605‑M1/27‑1 (MQ=255)
     tcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattgaccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCt                                                                                                                              <  1:134482‑M1/27‑1 (MQ=255)
      cctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcggggttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                             <  5:246323‑M1/28‑1 (MQ=255)
       ctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTc                                                                                                                             <  8:205839‑M1/28‑1 (MQ=255)
                          gccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGGAGCGCCGCCCGCGCGGCgcttc                                                                                                                            <  2:110020‑M1/29‑3 (MQ=255)
                          gccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGGAGCGCCGCCCGCGCGGCgcttc                                                                                                                            >  1:110020‑M1/102‑128 (MQ=255)
                             ccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                                <  8:18974‑M1/25‑1 (MQ=255)
                             ccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGgcggcg                                                                                                                                >  7:18974‑M1/99‑123 (MQ=255)
                                tggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCCCAATCTca                                                                                                                      <  3:300310‑M1/35‑1 (MQ=255)
                                tggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggccagaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCCCAATCTca                                                                                                                      >  4:300310‑M1/96‑130 (MQ=255)
                                           ctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCt                                                                                        <  2:85328‑M1/65‑1 (MQ=255)
                                            tcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCt                                                                                        <  6:36703‑M1/65‑1 (MQ=255)
                                                  tcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAAt                                                                                  >  1:85328‑M1/78‑148 (MQ=255)
                                                    gcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTc                                                                                <  3:217772‑M1/73‑1 (MQ=255)
                                                    gcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTc                                                                                >  4:217772‑M1/76‑148 (MQ=255)
                                                             aaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTc                                                                      <  4:86633‑M1/83‑1 (MQ=255)
                                                             aaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTc                                                                      <  1:49859‑M1/83‑1 (MQ=255)
                                                                       gtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAAt                                                            >  3:86633‑M1/57‑149 (MQ=255)
                                                                                  ctggtgcaggagggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAAcgccgc                                                 >  5:36703‑M1/46‑149 (MQ=255)
                                                                                           gagggactgtccgaaacaaacgtaggagcgggtttaCGCTCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaaa                                        >  5:44934‑M1/37‑149 (MQ=255)
                                                                                            agggactggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaaaa                                       >  2:49859‑M1/36‑149 (MQ=255)
                                                                                                actggccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTc                                    <  2:81492‑M1/117‑1 (MQ=255)
                                                                                                    gccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAAcg                               <  6:6489‑M1/122‑1 (MQ=255)
                                                                                                     ccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                    <  1:240885‑M1/37‑1 (MQ=255)
                                                                                                     ccgaaacaaacgtaggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTcaca                                                                                                                    >  2:240885‑M1/27‑63 (MQ=255)
                                                                                                                   ggagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGa                 >  1:81492‑M1/13‑148 (MQ=255)
                                                                                                                    gagcgggtttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGaa                <  5:275164‑M1/137‑1 (MQ=255)
                                                                                                                     agcggggttaCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCGCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGaaaa              >  6:275164‑M1/11‑149 (MQ=255)
                                                                                                                                CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa         >  7:163450‑M1/1‑143 (MQ=255)
                                                                                                                                CCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa         <  8:163450‑M1/143‑1 (MQ=255)
                                                                                                                                 CGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa  >  5:6489‑M1/1‑149 (MQ=255)
                                                                                                                                  atcgAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGAAAATCCCCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaaaa                                       >  7:80239‑M1/5‑111 (MQ=255)
                                                                                                                                  atcgAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCTGAAAATCCCCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaaaa                                       <  8:80239‑M1/107‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                    
ATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAA  >  NC_002947/4808047‑4808324

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