Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 3,845,705 +G intergenic (‑72/+26) 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.G95.3% 142.5 / 0.7 43intergenic (‑70/+28)PP_3394/PP_33953‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base G (25/16);  total (25/18)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCA  >  NC_002947/3845557‑3845839
                                                                                                                                                     |                                                                                                                                        
gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGATAAGACGTGTTTcgcg                                                                                                                                         >  1:141997/1‑148 (MQ=255)
  tgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCgg                                                                                                                                        >  1:283468/1‑147 (MQ=255)
       tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGAt                                                                                                                                   >  2:211837/1‑148 (MQ=255)
       tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCgggctg                                                                                                                                    <  8:264388/146‑5 (MQ=255)
          gCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATcc                                                                                                                                 >  5:144435/1‑147 (MQ=255)
             tGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCg                                                                                                                             >  1:90880/1‑148 (MQ=255)
               aTCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                                            >  1:24408/1‑147 (MQ=255)
                tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                                            >  1:43033/1‑146 (MQ=255)
                 cACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAg                                                                                                                         >  2:215759/1‑148 (MQ=255)
                       cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                   <  8:175290/148‑1 (MQ=255)
                       cAGGCTATTCAAGCCATCGCGCAGGCAGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACg                                                                                                                    <  7:30416/147‑1 (MQ=255)
                         ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATg                                                                                                                 <  2:310925/148‑1 (MQ=255)
                         ggCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATg                                                                                                                 >  4:111799/1‑148 (MQ=255)
                                   gCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTGGCTGACGATGGATCAGCCcg                                                                                                       <  1:127252/148‑1 (MQ=38)
                                           gcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTg                                                                                                <  5:79366/147‑1 (MQ=255)
                                           gcAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACATAGTCAGCCACGATGGATCAGCCCGCGCCCGGc                                                                                               >  2:166752/1‑148 (MQ=255)
                                                 cGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGc                                                                                          >  6:270472/1‑147 (MQ=255)
                                                         ttgattgTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACg                                                                                 >  8:162516/1‑148 (MQ=255)
                                                            ttttttattgttGTCATCGCGGTCTCTCTGTCGCAAGCGTGGCGCCAGACTAAGCGCTGGATGAGGGGCACAAAACCGGTGTTTCGCGGGCTGACCCATCGTCCGCCACGTTGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACgca                                                                               <  2:90880/143‑1 (MQ=255)
                                                                   ttgttgTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                           >  3:197802/1‑144 (MQ=255)
                                                                   ttgttgTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCaa                                                                           <  4:197802/144‑1 (MQ=255)
                                                                   ttgttgTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                                        <  7:162516/147‑1 (MQ=255)
                                                                    tgttgtCATGGGGGTCTCTGTGTCTGAGGCTTGGCGCCAGACTAAGTGGTGTTTGAGGGGCATAAAAGACGTGTTTGGGGGGCTGATCCATCGTCCGCCACGATGGATCTGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggt                                                                      <  8:63324/148‑1 (MQ=255)
                                                                       tgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtg                                                                     >  7:63324/1‑146 (MQ=255)
                                                                           aTCGCGGTCTCTCTGTCTTAAGCGTTGCGCCTGACTAAGAGCTGGTTGACGGGCTGAAACGCGGTGTTTGGCGGGCTTATCCTTCGTCAGCCAGGAGGCATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCaa                                                                <  1:14683/147‑1 (MQ=255)
                                                                            tCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGc                                                              <  2:43033/148‑1 (MQ=255)
                                                                             cgcgGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACCTCGTGAGCCCCGATGGATCACCCCCCGCCCGGCCAGGCCGGGGCACCGAGCAACAGGTGGTGaa                                                                >  5:224013/1‑145 (MQ=255)
                                                                                       ctgtctgAAGCGTTGCGCCTTACTAATAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgcc                                                       <  6:144435/144‑1 (MQ=255)
                                                                                       ctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTTAAACATAGTAAGCCACGATGGATAAGCCCCCCCCCTGCCAGGCCGGTGGACGCAGCAACAGGTGGCCCAGCCGGCgccgagc                                                   >  2:126504/1‑148 (MQ=255)
                                                                                       ctgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGCCAAGCAGGCgccgagc                                                   >  8:165643/1‑148 (MQ=255)
                                                                                             gAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCATGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACa                                              >  5:176766/1‑147 (MQ=255)
                                                                                               aGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                            >  8:186328/1‑147 (MQ=255)
                                                                                                        ccAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt                                  >  3:86260/1‑148 (MQ=255)
                                                                                                         cAGACTAAGAGCTGGATGACGGGCAGAAAAGCCGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGtgtg                                 <  1:211837/148‑1 (MQ=255)
                                                                                                              tAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATACATCGTCAGCCACGATGGATCAGCCCTAGCCATGCCAGGCCGGTGCAAGCAGCAACAGGTGGTAAAGCAGTCGCCGAGCCGACAGCGAGAGCGTGTGCCgg                             >  3:300794/1‑147 (MQ=255)
                                                                                                                   gCTGTATGACGGGCATAAAATATGTGTTTCGGGGGCTGATCCTTCGTCAGCCACTTTGTATCTGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgcg                        <  4:86260/147‑1 (MQ=255)
                                                                                                                     tGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCa                      >  5:331326/1‑147 (MQ=255)
                                                                                                                     tGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCAt                     >  2:205693/1‑148 (MQ=255)
                                                                                                                        ttGACGGGCATAAAAGACGTGTTTCGCGGGCTTATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                    <  6:99148/145‑1 (MQ=255)
                                                                                                                        aTGACGGGCAGAAAAGACGTGTTTCGCGGGCTTAAACATCCTGAGCCACCATGGAAAACACCCCGCCCAACCAGGCCCGTGCACGCAGCAACAAGTGGGCAAGCACGCCCCGAGCCGACACCGAGAGCGGGTTCCCGTCGAGCATGc                   >  4:179041/1‑147 (MQ=255)
                                                                                                                        aTGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                    >  5:99148/1‑146 (MQ=255)
                                                                                                                               gCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc            >  8:16756/1‑147 (MQ=255)
                                                                                                                                 aGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGa          <  5:270472/147‑1 (MQ=255)
                                                                                                                                  gAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGa             <  6:313126/143‑1 (MQ=255)
                                                                                                                                  gAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGa             >  5:313126/1‑143 (MQ=255)
                                                                                                                                         cGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCa  <  7:186328/147‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                        
GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCA  >  NC_002947/3845557‑3845839

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: