Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,777,922 G→T 14.2% intergenic (‑115/‑67) PP_2425 ← / → calA AraC family transcriptional regulator/coniferyl alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,777,9220GT14.2% 85.8 / 2.0 28intergenic (‑115/‑67)PP_2425/calAAraC family transcriptional regulator/coniferyl alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (14/10);  new base T (2/2);  total (16/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GTGACGGGTCAGTTGCATGGTGGACCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGCCATCGGTTATGCCGCTCAAACCCCCACCAGCCCCCTGGCCCCCATGACGTT  >  NC_002947/2777790‑2778050
                                                                                                                                    |                                                                                                                                
gTGACGGGTCAGTTGCATGGTGGACCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGc                                                                                                                               <  2:203350/136‑1 (MQ=255)
     gggTCAGTTGCATGGTGGACCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAggc                                                                                                                          >  1:234767/1‑136 (MQ=255)
           gTTGCATGGTGGACCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGaa                                                                                                                    <  2:523203/136‑1 (MQ=255)
           gTTGCATGGTGGACCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGaa                                                                                                                    >  2:416289/1‑136 (MQ=255)
           gTTGCATGGTGGACCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCAAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGaa                                                                                                                    >  2:13176/1‑136 (MQ=255)
               cATGGTGGACCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAacgac                                                                                                                >  1:401519/1‑136 (MQ=255)
                   gtggACCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACgcgc                                                                                                            <  1:223027/136‑1 (MQ=255)
                       aCCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTTGGCGAGGCGGCGAACGACGCGCaaaa                                                                                                        >  1:133557/1‑136 (MQ=255)
                         cAAGGGCAGAGAATCAGACTGCTTGATCTTCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATGAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAAtg                                                                                                      <  2:121367/136‑1 (MQ=255)
                         cAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAAtg                                                                                                      >  1:438512/1‑136 (MQ=255)
                                     aTCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAAt                                                                                          >  1:279895/1‑136 (MQ=255)
                                                 gATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTttcacttca                                                                              >  1:313923/1‑136 (MQ=255)
                                                 gATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTttcacttca                                                                              >  1:95896/1‑136 (MQ=255)
                                                 gATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTttcacttca                                                                              >  1:86351/1‑136 (MQ=255)
                                                    catcGCAGATTGAGGGCATGGCGCCAGCCATGGATTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTg                                                                           <  2:335212/136‑1 (MQ=255)
                                                         cAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAgg                                                                      <  1:533790/136‑1 (MQ=255)
                                                                  ggCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATg                                                             <  2:281235/136‑1 (MQ=255)
                                                                        gcgcCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTTGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATAcc                                                       <  2:133557/136‑1 (MQ=255)
                                                                         cgcCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCg                                                      >  2:122675/1‑136 (MQ=255)
                                                                          gccagccaTGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGc                                                     <  2:313923/136‑1 (MQ=255)
                                                                           ccagccaTGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGcc                                                    >  2:199541/1‑136 (MQ=255)
                                                                                      ttGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTTGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGCCATCGGTTATGc                                         >  2:508028/1‑136 (MQ=255)
                                                                                      ttGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGCGGCGGCGAACGACGCGCAAAATGAGTCGCCAGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTAAACCGCCATCGGTTATga                                         >  2:321937/1‑135 (MQ=255)
                                                                                                    gAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGCCATCGGTTATGCCGCTCAAACCCcca                           <  1:13176/136‑1 (MQ=255)
                                                                                                            catgcatgCCGGAGAATCAGTCGTTGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGCCATCGGTTATGCCGCTCAAACCCCCACCAGcccc                   <  1:508028/136‑1 (MQ=255)
                                                                                                             atgcatgcCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGCCATCGGTTATGCCGCTCAAACCCCCACCCGccccc                  >  2:144790/1‑136 (MQ=255)
                                                                                                             atgcatgcCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGCCATCGGTTATGCCGCTCAAACCCCCACCAGccccc                  <  2:401519/136‑1 (MQ=255)
                                                                                                                             cAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGCCATCGGTTATGCCGCTCAAGCCCCCACCAGCCCCCGGGCCCCCCTGACGtt  >  2:401633/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                
GTGACGGGTCAGTTGCATGGTGGACCAAGGGCAGAGAATCAGACAGCTTGATCATCGCAGATTGAGGGCATGGCGCCAGCCATGGCTTGGCATCTTCGTAGAATCAGGCATGCATGCCGGAGAATCAGTCGTGGGCGAGGCGGCGAACGACGCGCAAAATGTGTCGCCTGAATCGTTTCACTTCACTGCGAGGCCCTGCATGTATACCGCCATCGGTTATGCCGCTCAAACCCCCACCAGCCCCCTGGCCCCCATGACGTT  >  NC_002947/2777790‑2778050

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: