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.G94.3% 108.7 / 0.6 35intergenic (‑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 (20/13);  total (20/15)
Fisher's exact test for biased strand distribution p-value = 1.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.50e-01

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  NC_002947/3845557‑3845842
                                                                                                                                                     |                                                                                                                                           
gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTcgcg                                                                                                                                            <  5:147913/148‑1 (MQ=255)
gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTcgcg                                                                                                                                            <  3:354308/148‑1 (MQ=255)
 ctgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCggg                                                                                                                                          <  6:147469/149‑2 (MQ=255)
 ctgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATCGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCggg                                                                                                                                          <  6:36269/149‑2 (MQ=255)
  tgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCgggc                                                                                                                                         <  1:331328/149‑3 (MQ=255)
         ggCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATcc                                                                                                                                    <  1:173923/148‑1 (MQ=255)
         ggCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCa                                                                                                                                   <  6:92274/149‑1 (MQ=255)
               aTCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                                               >  8:250149/1‑147 (MQ=255)
                tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAg                                                                                                                            >  6:218939/1‑149 (MQ=255)
                           ctctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGAt                                                                                                                 <  6:257597/149‑1 (MQ=255)
                            tctGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGAt                                                                                                                 >  7:94571/1‑148 (MQ=255)
                              tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATAGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAg                                                                                                              >  2:325575/1‑149 (MQ=255)
                                cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGcc                                                                                                            >  3:277646/1‑149 (MQ=255)
                                     cATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCCGCAGCCACGATGGTGCAGCCCgcgc                                                                                                       >  2:339257/1‑149 (MQ=255)
                                       tCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGccc                                                                                                     >  5:363081/1‑149 (MQ=255)
                                               gCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGc                                                                                             >  8:39667/1‑149 (MQ=255)
                                                    attCATTGATTGTTTTTTTTGTCATCGCGGTCTCTGTGTCTTAAGCGTTGCGCCAGACTAATTGCTGGATGAGGGGCTGAAATGACGTGTTTCGCGGGCTGATCCATCTTCAGCCAGGATGGTTCAGCCCGCGCCCTGCCAGGCCGGt                                                                                         <  6:363081/147‑1 (MQ=255)
                                                                 tattgttGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                                           >  3:99749/1‑149 (MQ=255)
                                                                      ttgtCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtgga                                                                      >  7:91784/1‑148 (MQ=255)
                                                                           aTCGGGGTCTGTGTGTCTCAAGGGTTGCGCACTAATAATTGCTGTTTGAGGGGACTAAAAGACGTGTTGGGGGGGGTCACCAAGCGTCAGCCACGTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGc                                                                 <  4:277646/149‑1 (MQ=255)
                                                                           aTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGc                                                                 <  5:382162/149‑1 (MQ=255)
                                                                                    tctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccga                                                        >  7:164673/1‑149 (MQ=255)
                                                                                        tgtctgAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgccgagc                                                      <  2:116021/147‑1 (MQ=255)
                                                                                             gAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCAATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                               >  2:116796/1‑149 (MQ=255)
                                                                                              aaGCGTTGCGCCAGACTATGAGCTGGATGACGGGCATAAAAGAGGTGTTTCGGGGGCTGATCCATCGTCAGCCACTATGTATCAGCCCGCGCCCt                                                                                                    <  2:228639/95‑1 (MQ=255)
                                                                                              aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCt                                                                                                    >  1:228639/1‑95 (MQ=255)
                                                                                              aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAg                                                               >  7:292710/1‑132 (MQ=255)
                                                                                              aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAg                                                               <  8:292710/132‑1 (MQ=255)
                                                                                              aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATAAGCCCCCGCCCt                                                                                                    >  7:123193/1‑95 (MQ=255)
                                                                                                 cGTGGCGCCAGACTAAGAGCTGGATGTCGGGCATAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCgaga                                           <  6:293012/149‑1 (MQ=255)
                                                                                                        ccAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGtgtg                                    >  3:158037/1‑149 (MQ=255)
                                                                                                         cAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGc                                   >  3:251884/1‑149 (MQ=255)
                                                                                                                      ggATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                       >  3:238846/1‑148 (MQ=255)
                                                                                                                            cGGGCATAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGa                <  7:250149/149‑1 (MQ=255)
                                                                                                                             gggCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc               >  8:237324/1‑149 (MQ=255)
                                                                                                                                 aGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAt           >  7:25146/1‑149 (MQ=255)
                                                                                                                                        aCGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGATTTCACCaa    <  8:94571/149‑1 (MQ=255)
                                                                                                                                          gtgtTTCGCGGGCTTATCCATCGTCTCCAACTATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTTGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc  <  1:330841/149‑1 (MQ=255)
                                                                                                                                          gtgtTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc  <  4:251884/149‑1 (MQ=255)
                                                                                                                                          gtgtTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgc  <  1:325575/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGC  >  NC_002947/3845557‑3845842

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: