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.G95.6% 143.3 / ‑0.0 45intergenic (‑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 (24/19);  total (24/21)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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.

GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAA  >  NC_002947/3845557‑3845840
                                                                                                                                                      |                                                                                                                                          
gctgctGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCgg                                                                                                                                          >  7:4564/1‑149 (MQ=255)
  tgctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCgggc                                                                                                                                        <  3:361885/149‑3 (MQ=255)
    ctgTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCAGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCgggctg                                                                                                                                      <  5:224777/149‑5 (MQ=255)
      gTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGAGTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGAt                                                                                                                                    >  1:274772/1‑149 (MQ=255)
       tCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGAt                                                                                                                                    <  5:106647/148‑1 (MQ=255)
           cATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATc                                                                                                                               <  8:4564/149‑1 (MQ=255)
              gATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                                             >  1:111801/1‑148 (MQ=255)
                tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCGATGGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGt                                                                                                                             >  1:303011/1‑148 (MQ=255)
                tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAg                                                                                                                          >  6:177183/1‑149 (MQ=255)
                tCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAg                                                                                                                          >  2:45712/1‑149 (MQ=255)
                  aCTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGc                                                                                                                         <  5:254163/148‑1 (MQ=255)
                      gCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCAc                                                                                                                      <  4:9301/147‑1 (MQ=255)
                      gCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                    >  3:71781/1‑149 (MQ=255)
                       cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGAt                                                                                                                   <  4:71781/149‑1 (MQ=255)
                       cAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTCTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGa                                                                                                                    <  5:318440/148‑1 (MQ=255)
                        aGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATg                                                                                                                  >  2:100655/1‑149 (MQ=255)
                          gCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGa                                                                                                                <  6:246624/149‑1 (MQ=255)
                              tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATc                                                                                                              <  1:45712/147‑1 (MQ=255)
                              tGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATc                                                                                                              <  5:177183/147‑1 (MQ=255)
                                cAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGcc                                                                                                          >  6:106150/1‑149 (MQ=255)
                                 aaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATg                                                                                                                  >  7:277944/1‑140 (MQ=255)
                                 aaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATg                                                                                                                  >  7:283135/1‑140 (MQ=255)
                                 aaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATg                                                                                                                  <  8:277944/140‑1 (MQ=255)
                                 aaGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATg                                                                                                                  <  8:283135/140‑1 (MQ=255)
                                   gCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCcgc                                                                                                       <  6:338211/149‑1 (MQ=255)
                                         gcgcAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTg                                                                                                 >  5:170052/1‑149 (MQ=255)
                                          cgcAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGc                                                                                                >  6:66065/1‑149 (MQ=255)
                                               ggcGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGCCTAAGCGCTGGATGACGGGCAGAAAAGCCGTCTTTCGCGGGCTGATCCATCGAGAGACACGATCGATCAGCCAGCTCCCTGACAGGc                                                                                           >  8:119397/3‑149 (MQ=255)
                                                   aCTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGt                                                                                       <  5:106150/149‑1 (MQ=255)
                                                       cATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCAc                                                                                   >  1:24338/1‑149 (MQ=255)
                                                             ttgttattgttGTCATCG‑‑CTGTCTCTTTGTGTGAAGGGTTGGGCAATACTATGAGGTGTATTAGGGGCTCAAAAGCCGTGTTTCGGGGGCTTATCCTTCGTCAGCCACGATGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACgcagc                                                                              <  2:24338/148‑1 (MQ=255)
                                                                  tttgttgTAATCG‑‑CGGTCTGTCTTGCTTAAGCGTTCCGCAAGATTATGAGTTGTATTAGGGGCTGAAACCCGGTTTTGCGGGGCCGGATCCAACGATCCCCACTTTGTATCTGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAg                                                                         <  2:303011/147‑1 (MQ=255)
                                                                     gttgtCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAggtgg                                                                     >  3:191423/1‑149 (MQ=255)
                                                                                   gTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCacactg                                                         >  1:325012/1‑144 (MQ=255)
                                                                                    tctctcTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgcc                                                        <  6:170052/149‑1 (MQ=255)
                                                                                      tctctGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCgcc                                                        >  6:265666/1‑147 (MQ=255)
                                                                                                aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                             <  5:379212/148‑1 (MQ=255)
                                                                                                aaGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGc                                             >  8:103399/1‑148 (MQ=255)
                                                                                                         gcCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCgtgt                                   <  1:100655/149‑1 (MQ=255)
                                                                                                                        ggATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCTTCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATg                     >  5:78077/1‑148 (MQ=255)
                                                                                                                         gATGACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCCCGCCCTGCCAGGCCGGTGCAAGCCGCAACAGGTGGTCAAGGAGGCGCCGAGCCGCCAGCGAGAGCGTGTGCCGGTCGCGCCTg                     >  4:146629/1‑147 (MQ=255)
                                                                                                                            gACGGGCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATAAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCCCCGAGCCGACAGCGAGACCGTGTGCCGGTCGCGCATGCAGAt                >  5:253862/1‑149 (MQ=255)
                                                                                                                             aCGGGAATAAATGCCGTGTTGCGCGGGCTGATCCATGGTCAGCCACTTTGTATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGAt                <  2:325012/148‑1 (MQ=255)
                                                                                                                                ggCAGAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGCCAGCGAGAGCGGGTGCCGGTCGCGCATGCAGATGAc             >  3:33529/1‑148 (MQ=255)
                                                                                                                                    gAAAAGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAt         >  6:230390/1‑148 (MQ=255)
                                                                                                                                     aaaaGACGTGTTTCGCGGGCTGATCCCTCGTCAGCCCCGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCCCCGAGCCGCCAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAAtt        >  3:301485/1‑148 (MQ=255)
                                                                                                                                     aaaaGACGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTc       <  7:103399/149‑1 (MQ=255)
                                                                                                                                          aCGTGTTTCGCGGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCaa  <  2:274772/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                          
GCTGCTGTCGGCATGATCACTTGCAGGCTCTGCAAGCCATCGCGCAGGCCGACTTCATTGATTGTTATTGTTGTCATCG‑‑CGGTCTCTCTGTCTGAAGCGTTGCGCCAGACTAAGAGCTGGATGACGGGCAGAAAAGACGTGTTTCGCGGCTGATCCATCGTCAGCCACGATGGATCAGCCCGCGCCCTGCCAGGCCGGTGCACGCAGCAACAGGTGGTCAAGCAGGCGCCGAGCCGACAGCGAGAGCGTGTGCCGGTCGCGCATGCAGATGACGAATTCACCAA  >  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: