Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,040,838 C→T 13.8% intergenic (+2/‑43) gbd → / → PP_2651 4‑hydroxybutyrate dehydrogenase/MFS transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,040,8380CT13.8% 73.4 / 4.3 29intergenic (+2/‑43)gbd/PP_26514‑hydroxybutyrate dehydrogenase/MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (11/14);  new base T (2/2);  total (13/16)
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.57e-01

CCCGGACCTGTTCCCGCGCATCGTCCAGGGTGCGTTGGCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCGAACCCGTTCAGAAAGTACCAGGT  >  NC_002947/3040722‑3040963
                                                                                                                    |                                                                                                                             
cccGGACCTGTTCCCGCGCATCGTCCAGGGTGCGTTGGCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAAcc                                                                                                            <  1:15153/136‑1 (MQ=255)
        tGTTCCCGCGCATCGTCCAGGGTGCGTTGGCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATaacaa                                                                                                    <  1:325242/136‑1 (MQ=255)
        ggTTCCCGCGCATCGTCCAGGGTGCGTTGGCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATaacaa                                                                                                    <  1:355897/135‑1 (MQ=255)
               gcgcATCGTCCAGGGTGCGTTGGCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGt                                                                                             <  2:471927/136‑1 (MQ=255)
                                cGTTGGCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATAcc                                                                            <  1:253583/136‑1 (MQ=255)
                                    ggCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCaaaa                                                                        >  1:529359/1‑135 (MQ=255)
                                    ggCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCaaaa                                                                        >  1:3555/1‑135 (MQ=255)
                                    ggCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAat                                                                        >  2:459227/1‑136 (MQ=255)
                                    ggCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAat                                                                        >  1:430444/1‑136 (MQ=255)
                                      ccGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAatat                                                                      <  1:247774/136‑1 (MQ=255)
                                           cacaGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCa                                                                 >  1:127230/1‑136 (MQ=255)
                                             caGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAAc                                                               >  1:194008/1‑136 (MQ=255)
                                                      cACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAAAAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTg                                                      >  2:396810/1‑136 (MQ=255)
                                                       aCCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGc                                                     <  2:105400/136‑1 (MQ=255)
                                                          aaCCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCcg                                                  <  2:218919/136‑1 (MQ=255)
                                                                    aGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGaa                                        <  2:139460/136‑1 (MQ=255)
                                                                                   aCTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCg                         <  1:112448/136‑1 (MQ=255)
                                                                                   aCTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCACATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCg                         >  2:415183/1‑136 (MQ=255)
                                                                                            ggCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCGAACCCGTTc                <  1:390580/136‑1 (MQ=255)
                                                                                            ggCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCGAACCCGTTc                <  1:560723/136‑1 (MQ=255)
                                                                                                gctTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCGAACCCGTTCAGaa            >  1:546962/1‑136 (MQ=255)
                                                                                                  ttGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCGAACCCGTTCAGAAAg          <  1:55302/136‑1 (MQ=255)
                                                                                                    gAAGCGTCGATGTAAGTCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGt                                   <  2:382330/109‑1 (MQ=255)
                                                                                                    gAAGCGTCGATGTAAGTCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGt                                   <  2:385159/109‑1 (MQ=255)
                                                                                                    gAAGCGTCGATGTAAGTCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGt                                   >  1:382330/1‑109 (MQ=255)
                                                                                                    gAAGCGTCGATGTAAGTCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGt                                   >  1:385159/1‑109 (MQ=255)
                                                                                                    gAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCGAACCCGTTCAGAAAGTa        >  1:30271/1‑136 (MQ=255)
                                                                                                    gAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCGAACCCGTTCAGAAAGTa        >  2:480478/1‑136 (MQ=255)
                                                                                                          tCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCGAACCCGTTCAGAAAGTACCAGGt  <  2:430444/136‑1 (MQ=255)
                                                                                                                    |                                                                                                                             
CCCGGACCTGTTCCCGCGCATCGTCCAGGGTGCGTTGGCCGATCACAGCCACCGCACCAACCCGCGTGAGGCATCTGCCGAGGACTACCAGTGGCTGCTTGAAGCGTCGATGTAAGCCGTTACCTGCTGTTCAACCAATAACAACAAAAGTGAACGCCCATGAATACCAAATATGCGCAACCCGGGCTTGCCCGCGAACCTGAAGCGGTCGCCAGTTCGAACCCGTTCAGAAAGTACCAGGT  >  NC_002947/3040722‑3040963

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: