Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,845,703:1 +G 100% intergenic (‑70/+28) PP_3394 ← / ← PP_3395 3‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,845,7031.G97.5% 130.6 / ‑3.1 40intergenic (‑70/+28)PP_3394/PP_33953‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (24/15);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAA  >  NC_002947/3845557‑3845840
                                                                                                                                                    |                                                                                                                                         
gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑Tcgc‑g                                                                                                                                          <  4:40980/148‑1 (MQ=255)
gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑Tcgc‑g                                                                                                                                          <  4:94662/148‑1 (MQ=255)
gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGC‑gg                                                                                                                                         >  5:32280/1‑149 (MQ=255)
 ctgctgTCTGCATTATCACTTGCAGGCTCTGCATGCCATCTCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGTGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGC‑ggg                                                                                                                                        <  1:275904/149‑2 (MQ=255)
     tgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGC‑gggctga                                                                                                                                    >  3:228319/1‑144 (MQ=255)
      gTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGAt                                                                                                                                    >  5:189437/1‑149 (MQ=255)
       tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATc                                                                                                                                   >  6:165321/1‑149 (MQ=255)
          gCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCAt                                                                                                                                <  6:189437/149‑1 (MQ=255)
          gCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCAt                                                                                                                                <  2:39405/149‑1 (MQ=255)
                  aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGc                                                                                                                         >  2:96943/1‑148 (MQ=255)
                           ctctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGAt                                                                                                               <  4:228319/149‑1 (MQ=255)
                            tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGAt                                                                                                               >  1:281883/1‑148 (MQ=255)
                              tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGC‑gggctga                                                                                                                                    >  3:188739/1‑119 (MQ=255)
                              tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGC‑gggctga                                                                                                                                    <  4:188739/124‑6 (MQ=255)
                                cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGc                                                                                                           <  4:322633/148‑1 (MQ=255)
                                      aTCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGcc                                                                                                    >  6:341137/1‑149 (MQ=255)
                                               gCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGGCGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCa                                                                                                             <  4:41291/131‑1 (MQ=255)
                                               gCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGGCGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCa                                                                                                             >  3:41291/1‑131 (MQ=255)
                                               gCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCa                                                                                                             <  2:162908/131‑1 (MQ=255)
                                               gCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCa                                                                                                             >  1:162908/1‑131 (MQ=255)
                                                        attgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACg                                                                                  >  3:178287/1‑149 (MQ=255)
                                                         ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGC‑gggctg                                                                                                                                     <  4:321498/96‑5 (MQ=255)
                                                         ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGC‑gggctg                                                                                                                                     >  3:321498/1‑92 (MQ=255)
                                                          tgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACg                                                                                  >  1:351777/1‑147 (MQ=255)
                                                                         tCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTc                                                                   >  3:60127/1‑147 (MQ=255)
                                                                           aTCGGGGTGTCTCTGTCTGAAGCTTGGCGCCATATTATTCGGTGGATGCGGGGAAGAAAAGCGGGGTTGGGGG‑GGCTGATCCATGGTCAGCCACGTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAg                                                                <  4:60127/148‑1 (MQ=255)
                                                                                ggTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAgg                                                                        >  8:299925/1‑135 (MQ=255)
                                                                                ggTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAgg                                                                        <  7:299925/135‑1 (MQ=255)
                                                                                      tctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                    >  5:174231/1‑149 (MQ=255)
                                                                                         atcggAAGCGTTGCGCCAGAATAAGCGCTGGATGACGGGCAGAAAAGACTTGTT‑TCGCGGGCTGATCCATCGTCTGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagccga                                                 <  7:79008/145‑1 (MQ=255)
                                                                                          tctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagcc                                                   <  6:174231/146‑1 (MQ=255)
                                                                                          tctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGCGCAGAAAAGACGTGTT‑TCGCGGGCAGATCCCTCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCCGCAACAGGTGGTCAACCAGGCcccccgccgc                                                 >  5:25718/1‑147 (MQ=255)
                                                                                            tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                                <  4:263909/83‑1 (MQ=255)
                                                                                            tgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                                >  3:263909/1‑83 (MQ=255)
                                                                                              aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGGCCCGATGGGTCAGCCCGCGCCCTGGCAGGCCGGTGGACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGGCAGCg                                            >  5:231374/1‑149 (MQ=255)
                                                                                              aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGCTGCACGCAGCAACAGGTGGTCAAGCAGGCgcccagccca                                                 >  8:79008/1‑142 (MQ=255)
                                                                                                       gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgcag                                                                               <  3:231519/105‑1 (MQ=255)
                                                                                                       gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgcag                                                                               >  4:231519/1‑105 (MQ=255)
                                                                                                       gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt                                   >  1:93052/1‑149 (MQ=255)
                                                                                                        ccAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGtgtg                                  >  3:117311/1‑149 (MQ=255)
                                                                                                                 gagCTGGATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGGCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgcg                         <  1:96943/149‑1 (MQ=255)
                                                                                                                      ggATGACGGGCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                     >  5:360712/1‑148 (MQ=255)
                                                                                                                             gggCAGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc             >  5:111729/1‑149 (MQ=255)
                                                                                                                                 aGAAAAGAGTTGTT‑TCGCGGGCTGATCCATCGTCAGCCACTATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                            <  8:247323/130‑1 (MQ=255)
                                                                                                                                 aGAAAAGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                            >  7:247323/1‑130 (MQ=255)
                                                                                                                                   aaaaGACGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGTCGGTCGCGCATGCAGATGACGAATTc       >  7:25091/1‑149 (MQ=255)
                                                                                                                                        aCGTGTT‑TCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCaa  <  5:341137/149‑1 (MQ=255)
                                                                                                                                             tt‑tCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAg                  >  3:231735/1‑128 (MQ=255)
                                                                                                                                             tt‑tCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAg                  <  4:231735/128‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                         
GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTT‑TCGC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAA  >  NC_002947/3845557‑3845840

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: ‑