Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,388,541 A→G 12.9% intergenic (‑46/‑236) PP_1209 ← / → dps cold‑shock protein/DNA‑binding stress protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,388,5410AG12.9% 76.1 / 3.6 31intergenic (‑46/‑236)PP_1209/dpscold‑shock protein/DNA‑binding stress protein
Reads supporting (aligned to +/- strand):  ref base A (16/11);  new base G (2/2);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

TCGCATCGGTTTTCAGGCTGGGTATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGCAGCCTACTCAAACCGCAGGCCGCCCAGTGGGCAATGTATTG  >  NC_002947/1388417‑1388690
                                                                                                                            |                                                                                                                                                      
tCGCATCGGTTTTCAGGCTGGGTATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAgc                                                                                                                               <  2:374092/150‑1 (MQ=255)
               ggCTGGGTATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCCCAGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAg                                                                                                                >  1:37211/1‑150 (MQ=255)
               ggCTGGGTATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATgaga                                                                                                                                       <  1:245293/127‑1 (MQ=255)
               ggCTGGGTATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATgaga                                                                                                                                       >  2:245293/1‑127 (MQ=255)
                gCTGGGTATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATga                                                                                                                                         <  2:470585/124‑1 (MQ=255)
                gCTGGGTATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATga                                                                                                                                         >  1:470585/1‑124 (MQ=255)
                  tGGGTATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCGGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGcc                                                                                                             >  1:476640/1‑150 (MQ=255)
                      tATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGcc                                                                                                         >  1:467241/1‑150 (MQ=255)
                        tGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTa                                                                                                       <  1:300549/150‑1 (MQ=255)
                        tGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTa                                                                                                       >  2:300549/1‑150 (MQ=255)
                                     caAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAAt                                                                                          <  2:467241/150‑1 (MQ=255)
                                      aaaGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATg                                                                                         >  1:282335/1‑150 (MQ=255)
                                               aCGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGc                                                                               >  1:138626/1‑150 (MQ=255)
                                                        aCCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCCGCAATGACAGTCCGCACGGGTGGc                                                                      <  2:193736/150‑1 (MQ=255)
                                                        aCCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGGTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGc                                                                      >  2:326974/1‑150 (MQ=255)
                                                          cAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGcac                                                                    >  1:555213/1‑150 (MQ=255)
                                                                     tCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTTGGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTAttt                                                         >  1:344812/1‑150 (MQ=255)
                                                                     tCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTAttt                                                         >  2:10029/1‑150 (MQ=255)
                                                                         gAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACa                                                     <  2:282335/150‑1 (MQ=255)
                                                                                  gTTGGCTAACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCg                                            <  1:260604/150‑1 (MQ=255)
                                                                                   ttGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGc                                           <  2:555213/150‑1 (MQ=255)
                                                                                         cACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCGGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGCAGCCTa                                     <  2:476640/150‑1 (MQ=255)
                                                                                           cTGTTGATAAGGCCGATCACAGGAAGAGCTGGCGGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGCa                                          <  2:162613/143‑1 (MQ=255)
                                                                                           cTGTTGATAAGGCCGATCACAGGAAGAGCTGGCGGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGCa                                          >  1:162613/1‑143 (MQ=255)
                                                                                                  tAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGTACACTGCCTATTTAACAGTAAAGCCGCAGCCTACTCAAACCg                            >  1:122365/1‑150 (MQ=255)
                                                                                                    aGGCCGATCACAGGAAGCGCTGGCAGCACGGTGCGCATGAGATATAGCGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATCACAGTCCGCACGGGTGGCACACTGCCCATTTAACACTAAAGCCGCAGCCTACTCAAACCGCa                          <  2:159566/150‑1 (MQ=255)
                                                                                                     ggCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGCAGCCTACTCAAACCGCAg                         >  1:51770/1‑150 (MQ=255)
                                                                                                        cGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGc                                           >  1:469422/1‑129 (MQ=255)
                                                                                                        cGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGc                                           <  2:469422/129‑1 (MQ=255)
                                                                                                         gATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGCAGGAAGCAGTCCGCACGGTTGGCACACTGCATATTTCCCCGTATCCCCATAGCCTACGCAAACCGCAGgccg                     >  1:434282/1‑150 (MQ=255)
                                                                                                          aTCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGCAGCCTACTCAAACCGCAGgccgc                    >  1:325802/1‑150 (MQ=255)
                                                                                                                            agcaCGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGCAGCCTACTCAAACCGCAGGCCGCCCAGTGGGCAATGTATTg  >  2:217381/1‑150 (MQ=255)
                                                                                                                            |                                                                                                                                                      
TCGCATCGGTTTTCAGGCTGGGTATGAGCGATAGCAGCAAAGCTGCCACGCCCGTCACCAGATGGACGATCTTGAACATGGGGTTGGCTCACTGTTGATAAGGCCGATCACAGGAAGAGCTGGCAGCACGGTGCGCATGAGATATAGAGCGCGAACACGTGGCAGGCCAGCCTATGCACCCGGCAATGACAGTCCGCACGGGTGGCACACTGCCTATTTAACAGTAAAGCCGCAGCCTACTCAAACCGCAGGCCGCCCAGTGGGCAATGTATTG  >  NC_002947/1388417‑1388690

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