Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,845,720 G→C 18.9% intergenic (‑87/+11) 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,7200GC18.9% 68.8 / 12.4 37intergenic (‑87/+11)PP_3394/PP_33953‑hydroxy‑3‑methylglutaryl‑CoA lyase/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (21/9);  new base C (2/5);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 8.01e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

TTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGAC  >  NC_002947/3845577‑3845868
                                                                                                                                                 |                                                                                                                                                    
ttGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCAc                                                                                                                                                  >  1:168324/1‑149 (MQ=255)
  gCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGa                                                                                                                                                >  1:178745/1‑149 (MQ=255)
   cAGTCTCTGCATGCAATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCTGTATCTCTGTATGAAGCGTTGCGCAAGACTAAGAGCTGGATGACGGGCAGAAAAGACGGGTTTC‑GCTGGCTGATCCATCGTCAGCCACGAt                                                                                                                                               <  8:59334/149‑1 (MQ=255)
   cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGAt                                                                                                                                               <  3:205292/149‑1 (MQ=255)
         ctGAAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTGGCTGACGATGGATCa                                                                                                                                         <  7:126078/149‑1 (MQ=255)
          tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAg                                                                                                                                        >  5:316105/1‑149 (MQ=255)
            cAAGCCATCGCGCAGGCCGACTTCATTGAGTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTGGCTGACGATGGATCa                                                                                                                                         <  3:239404/146‑1 (MQ=255)
                  aTCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTGGCTGACGATGGATCAg                                                                                                                                        <  3:224682/141‑1 (MQ=38)
                  aTCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTGGCTGACGATGGATCAg                                                                                                                                        >  4:224682/1‑141 (MQ=38)
                  aTCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGcc                                                                                                                                >  6:182233/1‑149 (MQ=255)
                             cGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCg                                                                                                                     >  4:181255/1‑149 (MQ=255)
                                                 gttgtCATCGCGGTTTGTCTGTCTTATGCGTTGCGCCAGACTATTAGCTGTTTGACGGGCAGAAATGAGGTTTTTCGGG‑GGCTTATCGATGGTCTGCAACGTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtg                                                                                                  <  2:267151/148‑1 (MQ=255)
                                                     tCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTc                                                                                               >  3:68565/1‑147 (MQ=255)
                                                     tCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCAAGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTc                                                                                               >  8:154609/1‑147 (MQ=255)
                                                       aTCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGc                                                                                           <  7:154609/149‑1 (MQ=255)
                                                            ggTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGgcg                                                                                      >  5:384074/1‑149 (MQ=255)
                                                                                   gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATACATCGTAAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTAAAGCAGGCGACGAGCCGACAGCGAGAGCgtgt                                                               >  4:185390/1‑149 (MQ=255)
                                                                                       gACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCg                                                           <  6:384074/149‑1 (MQ=255)
                                                                                            agagCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgc                                                      >  5:277294/1‑149 (MQ=255)
                                                                                            agagCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTcgc                                                      >  3:354639/1‑149 (MQ=255)
                                                                                            agagCTGGATGACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTc                                                        <  5:182233/147‑1 (MQ=255)
                                                                                                      gACGGGCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt                                            >  4:44013/1‑149 (MQ=255)
                                                                                                       aggggCACAAATCACGTGTTTCGGG‑GGCTGATCCATCGTCTCCCACTATGTTTAAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt                                            <  2:330060/146‑1 (MQ=255)
                                                                                                         gggCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGAc                                         >  5:357771/1‑149 (MQ=255)
                                                                                                           gCAGAAAAGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACg                                        >  8:109674/1‑148 (MQ=255)
                                                                                                               aaaaGACGTGTTTC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGGCTGCAGATGACGAAtag                                   >  3:100512/1‑147 (MQ=255)
                                                                                                                         tttC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgcgc                          <  7:113107/148‑1 (MQ=255)
                                                                                                                         tttC‑GCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAgcgc                          >  8:113107/1‑148 (MQ=255)
                                                                                                                                   cTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGACAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtca                <  6:357771/149‑1 (MQ=255)
                                                                                                                                    tGATCCATCGTGGCTGACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGt                                                         <  5:48512/107‑1 (MQ=38)
                                                                                                                                    tGATCCATCGTGGCTGACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGt                                                         >  6:48512/1‑107 (MQ=38)
                                                                                                                                      aTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagt              >  4:173977/1‑148 (MQ=255)
                                                                                                                                      aTCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCgtcagt              >  3:151681/1‑148 (MQ=255)
                                                                                                                                        ccATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGtcagtcag           >  6:234611/1‑149 (MQ=255)
                                                                                                                                        ccATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGtcagtcag           >  1:189565/1‑149 (MQ=255)
                                                                                                                                           tCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAgcg         >  7:243191/1‑148 (MQ=255)
                                                                                                                                              tCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAAt     <  8:243191/149‑1 (MQ=255)
                                                                                                                                                 gCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGAc  <  3:44013/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGCGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTC‑GC‑GGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAAGCGCCCAGTCGTCAGTCAGCGCAATGAC  >  NC_002947/3845577‑3845868

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