Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,496,668 A→G 14.9% intergenic (+18/+34) PP_3987 → / ← PP_5651 site‑specific recombinase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,496,6680AG14.9% 59.2 / 4.2 28intergenic (+18/+34)PP_3987/PP_5651site‑specific recombinase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (11/12);  new base G (2/2);  total (13/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 = 4.38e-01

CCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTG  >  NC_002947/4496527‑4496808
                                                                                                                                             |                                                                                                                                            
ccccAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTccc                                                                                                                                       <  8:104654/149‑1 (MQ=255)
   cAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGc                                                                                                                                    <  5:140489/149‑1 (MQ=255)
    aaGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCt                                                                                                                                   >  7:82331/1‑149 (MQ=255)
      ggcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGAGCGTACCCCGCttc                                                                                                                                 <  7:69115/149‑2 (MQ=255)
      ggcggcAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGAGCGTACCCCGCttc                                                                                                                                 <  6:228751/149‑2 (MQ=255)
              gCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAACCGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACaa                                                                                                                         >  5:224133/1‑149 (MQ=255)
                   aTCGGCGAGTGGGGTCACGCACGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTcac                                                                                                                    <  6:224133/149‑1 (MQ=255)
                                        aGCCGCCTATTGGGTCACACTGACAAGCTCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATAc                                                                                               <  3:98484/149‑1 (MQ=255)
                                               tATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTAACAACTGATCTATTGAGCTATACGCATTg                                                                                         >  3:176778/1‑148 (MQ=255)
                                                         cacTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTa                                                                              <  7:63390/149‑1 (MQ=255)
                                                              aCAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTtg                                                                         <  8:82331/149‑1 (MQ=255)
                                                                     cATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGctct                                                                  >  4:294948/1‑149 (MQ=255)
                                                                      aTCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTg                                                                 >  3:200971/1‑149 (MQ=255)
                                                                         acacaAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCa                                                              >  8:42830/1‑149 (MQ=255)
                                                                          cacaAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCAc                                                             >  4:162382/1‑149 (MQ=255)
                                                                             aaaCGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATACACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCAccc                                                           >  4:102693/1‑148 (MQ=255)
                                                                               aCGTTATATAGAGGAGTAGGTGATATCTTACAGCCAACACGATAGGTAATTCGAAGCGTGGTTCGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACGCAat                                                        <  8:290378/149‑1 (MQ=255)
                                                                                cGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTAGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAatc                                                       <  3:102693/149‑1 (MQ=255)
                                                                                     tataGACTAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATcatcat                                                  <  3:237091/149‑1 (MQ=255)
                                                                                        aGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGctct                                                                  <  7:205051/130‑1 (MQ=255)
                                                                                        aGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGctct                                                                  >  8:205051/1‑130 (MQ=255)
                                                                                              gTAGGTGAAATCGTACATCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCga                                         <  3:162382/149‑1 (MQ=255)
                                                                                                   tGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGaa                                    >  1:189284/1‑149 (MQ=255)
                                                                                                                 cAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCt                      <  3:294948/149‑1 (MQ=255)
                                                                                                                               aTTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTg         <  4:200971/148‑1 (MQ=255)
                                                                                                                                ttCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGa       >  2:157733/1‑149 (MQ=255)
                                                                                                                                 tcaaaGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACATCGCACGCCTGAGTAGCTCTTTTGAc      >  4:43094/4‑149 (MQ=255)
                                                                                                                                     aGCGGGGAGCGTACCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTg  >  7:156154/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                            
CCCCAAGGCGGCAAGCGAAATCGGCGATTTGGGTCACGCAAGCCGCCTATTGGGTCACACTGACAAGCGCATCACACAAACGGTATATAGACGAGTAGGTGAAATCGTACAGCCAACACGATAGGTAATTCGAAGCGGGGTACGCTCCCCGCTTTGAGCACAATTCACAACTGATCTATTGAGCTATACGCATTGTATCAATCTTAGCTTGCTGCTCTGCCACCCAATCATCATAGGTTTCGAGAGAATACCTCGCACGCCTGAGTAGCTCTTTGGACCTTG  >  NC_002947/4496527‑4496808

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: