Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,321,111 A→G 22.7% intergenic (‑303/‑115) PP_3791 ← / → PP_5629 site‑specific recombinase/protease/amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,321,1110AG22.7% 34.2 / 7.4 22intergenic (‑303/‑115)PP_3791/PP_5629site‑specific recombinase/protease/amidase
Reads supporting (aligned to +/- strand):  ref base A (8/9);  new base G (2/3);  total (10/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 = 6.52e-01

GCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACTGACACCGGGTTT  >  NC_002947/4320976‑4321260
                                                                                                                                       |                                                                                                                                                     
gCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGGATTCAATATcgaaa                                                                                                                                         <  2:521895/150‑4 (MQ=255)
                gCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCAtcct                                                                                                                         >  1:254199/1‑150 (MQ=255)
                       cTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCt                                                                                                                  >  2:82372/1‑150 (MQ=255)
                           gCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGa                                                                                                              <  1:510957/150‑1 (MQ=255)
                                     tCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAg                                                                                                    >  1:316259/1‑150 (MQ=255)
                                           tCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAAc                                                                                              <  2:520077/150‑1 (MQ=255)
                                                gTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGc                                                                                                            <  1:243637/131‑1 (MQ=255)
                                                gTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGc                                                                                                            >  2:243637/1‑131 (MQ=255)
                                                 tcaatcaaGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGCCCCTGGTGCCAGCAGCACAGTATTCTTCATTTCGATATTGAATCTAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTgct                                                                                        >  2:115246/1‑148 (MQ=255)
                                                       aaGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGt                                                                                  <  2:254199/150‑1 (MQ=255)
                                                                                             ggCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGaa                                                                                                                                  <  2:112004/64‑1 (MQ=255)
                                                                                             ggCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGaa                                                                                                                                  >  1:112004/1‑64 (MQ=255)
                                                                                                ccTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTg                                         >  2:123343/1‑150 (MQ=255)
                                                                                                          gcagcaCAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCa                               <  1:123343/150‑1 (MQ=255)
                                                                                                              cacaGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCAcgcg                           >  1:373827/1‑150 (MQ=255)
                                                                                                               acaGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCCCTAGCTGCACCGTGCTTTCTGCTTGAAGCCCCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGt                          <  2:33803/150‑1 (MQ=255)
                                                                                                                       cTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACACTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGa                  <  2:316259/150‑1 (MQ=255)
                                                                                                                         tAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATa                <  1:82372/150‑1 (MQ=255)
                                                                                                                              gcGATATTGGATTCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGa                  <  1:376243/142‑1 (MQ=255)
                                                                                                                              gcGATATTGGATTCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGa                  >  2:376243/2‑143 (MQ=255)
                                                                                                                                  tatTGGATTCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTgg                                        >  1:416214/1‑117 (MQ=255)
                                                                                                                                  tatTGGATTCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTgg                                        <  2:416214/117‑1 (MQ=255)
                                                                                                                                       aaTCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACTGACACCGGGttt  <  1:277951/150‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                                     
GCGTGTTATCTAAGATGCGTCAGCTTCGCTTTCGCCCTCATAATCGTAGTCAATCAAGCTTGCCGACACTTGACAAATTTCAAATTACCCTAAGGCCCTGGTGCCAGCAGCACAGTATTCTTAATTTCGATATTGAATCCAATATCGCGACCGGGAACATCATCCTCCCTTCTATGAGCACCGCTAGCTGAACCTTGATGTCTGTTTGAAGCCGCCGTACGCTATCGGCATTCATGGAGGAATCTGGCTGATGGCACGCGTGGGTTTGATACTGACACCGGGTTT  >  NC_002947/4320976‑4321260

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