Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,319,181 G→A 13.2% intergenic (‑102/+90) treSA ← / ← PP_2919 trehalose synthase A/LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,319,1810GA13.2% 64.0 / 2.3 30intergenic (‑102/+90)treSA/PP_2919trehalose synthase A/LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (13/13);  new base A (2/2);  total (15/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

CGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTGCGC  >  NC_002947/3319043‑3319328
                                                                                                                                          |                                                                                                                                                    
cGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGc                                                                                                                                           <  1:106070/150‑1 (MQ=255)
                             tGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTAGCGCTGCAc                                                                                                             >  1:154125/1‑150 (MQ=255)
                               ggTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCt                                                                                                           >  1:202791/1‑150 (MQ=255)
                                           tCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGCCGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGgc                                                                                               <  2:281318/150‑1 (MQ=255)
                                                           gCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGc                                                                               <  2:202791/150‑1 (MQ=255)
                                                              ccGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCgg                                                                            >  1:455603/1‑150 (MQ=255)
                                                                     aGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCTCGAAGAGGCCGGGCCTGCa                                                                     >  2:402458/1‑150 (MQ=255)
                                                                       gTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCaca                                                                   <  2:455603/150‑3 (MQ=255)
                                                                          aTGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCCAATGCGGACATACCGACGCCAGTTACGATGTGTGCTGTGGGATCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGc                                                                >  1:367670/1‑150 (MQ=255)
                                                                                          gatctcGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCGGCACCGGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGc                                                <  1:330196/145‑1 (MQ=255)
                                                                                                gTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGGTCCCACAGGCGGCCGTACCGCTGAAGATTTGCGCTGCACTTGTGGGAAAGGGCGCGGCCGCGAAGAGGCCGGGCCTGCGAAGGCAGGATCCGCGGCGGGCCggatc                                          >  2:330196/1‑146 (MQ=255)
                                                                                                       tCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCACACAGGCGGACATACCGCCGAAGATTTGGGCTGCCCCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACtgctgc                                   >  2:389713/1‑150 (MQ=255)
                                                                                                              gggTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGa                                                                                                                                        <  2:231080/43‑1 (MQ=255)
                                                                                                              gggTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGa                                                                                                                                        >  1:231080/1‑43 (MQ=255)
                                                                                                                  gtTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTgg                                                                                                       >  2:125051/1‑71 (MQ=255)
                                                                                                                  gtTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTgg                                                                                                       <  1:125051/71‑1 (MQ=255)
                                                                                                                  gtTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGaaa                        <  1:102495/150‑1 (MQ=255)
                                                                                                                   ttGGTTTCTTCGCGGGCTTGCCCACTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCgg                                                  >  1:440556/1‑123 (MQ=255)
                                                                                                                   ttGGTTTCTTCGCGGGCTTGCCCACTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCgg                                                  <  2:440556/123‑1 (MQ=255)
                                                                                                                   ttGGTTTCTTCGCGGGCTTGCCCACTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCgg                                                  >  1:440555/1‑123 (MQ=255)
                                                                                                                   ttGGTTTCTTCGCGGGCTTGCCCACTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCgg                                                  <  2:440555/123‑1 (MQ=255)
                                                                                                                         tcttcGCGGGCGCGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTgg                                                                                                       <  1:503768/64‑1 (MQ=37)
                                                                                                                         tcttcGCGGGCGCGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTgg                                                                                                       >  2:503768/1‑64 (MQ=37)
                                                                                                                           ttcGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGGGGAGGGCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATg               <  1:177932/150‑1 (MQ=255)
                                                                                                                           ttcGCGGGCTTGCCCGCGCCCACAGGCGGACAGACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGCCAACTGCTGCGGTGCCAGAAACGCATCATg               >  2:6802/1‑150 (MQ=255)
                                                                                                                              gcgGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCTCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGaa            <  1:402458/150‑1 (MQ=255)
                                                                                                                              gcgGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGaa            <  2:327280/150‑1 (MQ=255)
                                                                                                                              gcgGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTCCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGaa            <  2:367670/150‑1 (MQ=255)
                                                                                                                                       cccGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTgcg   >  1:172867/1‑150 (MQ=255)
                                                                                                                                        ccGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTTCTGCGGTGCCAGAAACGCATCATGGAAGTAATTgcgc  >  2:119558/1‑150 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
CGAGCCATTTGACGTATGACGGGTCGGGCTGGGTCATTTCGCGTCCTCTTTGCGTGTCGGCACCGGCTCAGGTTATGACGGCGCTGCGGCCTGAACGTTCGGCTCAAGATGGGTGTTGGTTTCTTCGCGGGCTTGCCCGCTCCCACAGGCGGACATACCGCTGAAGATTTGCGCTGCACCTGTGGGAGCGGGCGCGCCCGCGAAGAGGCCGGGCCTGCAAAGGCAGGATCAGCGCCGGGCCGACAACTGCTGCGGTGCCAGAAACGCATCATGGAAGTAATTGCGC  >  NC_002947/3319043‑3319328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: