Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,287,276 T→C 14.3% intergenic (+17/+24) PP_2889 → / ← xthA transmembrane regulator PrtR/exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,287,2760TC14.3% 59.4 / 3.9 28intergenic (+17/+24)PP_2889/xthAtransmembrane regulator PrtR/exodeoxyribonuclease III/apurinic/apyrimidinic endodeoxyribonuclease VI
Reads supporting (aligned to +/- strand):  ref base T (16/8);  new base C (2/2);  total (18/10)
Fisher's exact test for biased strand distribution p-value = 6.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

AGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGGTC  >  NC_002947/3287150‑3287419
                                                                                                                               |                                                                                                                                               
aGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGACGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCt                                                                                                                            >  1:313819/1‑148 (MQ=255)
                                 ataACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGt                                                                                          >  3:163118/1‑149 (MQ=255)
                                 ataACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCGAATCAGCTCAACTCAAGCCAGATCGGCGCATGGt                                                                                          >  7:161545/1‑149 (MQ=255)
                                          ccTTGGTCAGCGCGGCGGACGATGTGCGTGGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGc                                                                                  >  1:43628/1‑149 (MQ=255)
                                            ttGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCttt                                                                               >  7:274006/1‑149 (MQ=255)
                                                tCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCa                                                                           >  7:259706/1‑149 (MQ=255)
                                                  aGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATc                                                                         >  3:98758/1‑149 (MQ=255)
                                                                tgtgCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCa                                                             >  7:19691/1‑147 (MQ=255)
                                                                          ggCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGcc                                                  >  6:156454/1‑149 (MQ=255)
                                                                                      ggAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTc                                                                             >  4:337655/1‑110 (MQ=255)
                                                                                      ggAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTc                                                                             <  3:337655/110‑1 (MQ=255)
                                                                                          ccGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATgcgc                                  <  8:19691/148‑1 (MQ=255)
                                                                                            gTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATa                                                           >  8:228210/1‑122 (MQ=255)
                                                                                            gTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATa                                                           <  7:228210/122‑1 (MQ=255)
                                                                                                        gTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAgct                        >  7:169591/1‑145 (MQ=255)
                                                                                                        gTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAgct                        <  8:169591/145‑1 (MQ=255)
                                                                                                           gTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTcc                                                                            <  6:4901/90‑1 (MQ=255)
                                                                                                           gTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTcc                                                                            >  5:4901/1‑90 (MQ=255)
                                                                                                            tAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAg                >  6:78075/1‑149 (MQ=255)
                                                                                                                tgtgGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCa            <  4:163118/149‑1 (MQ=255)
                                                                                                                tgggggCCCAAGGGGTCCCCCATGGGCCCGGCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACACCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCa            <  8:161545/146‑1 (MQ=255)
                                                                                                                   gggCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGa         >  3:138759/1‑149 (MQ=255)
                                                                                                                     gCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATc       <  8:274006/149‑1 (MQ=255)
                                                                                                                     gCCCAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGc                               <  2:306287/125‑1 (MQ=255)
                                                                                                                     gCCCAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGAt        >  2:119531/1‑148 (MQ=255)
                                                                                                                      cccAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGc              >  6:127815/1‑141 (MQ=255)
                                                                                                                      cccAAGGTGCACACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGc              <  5:127815/141‑1 (MQ=255)
                                                                                                                          aGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGGTc  >  4:111704/1‑149 (MQ=255)
                                                                                                                               |                                                                                                                                               
AGCTGGAGGCGCGGTACTGGTCGCATGGGGCTTATAACTTTGCCTTGGTCAGCGCGGCGGACGATGTGCGTGGGGCGGGGGTGGGGGAGGCGTTGAGGTTGGGGTTGTAAGTGTGGGCCCAAGGTGTGCACCTTGGGCCCGTCACCCTAATCAGCTCAACTCCAGCCAGATCGGCGCATGGTCCGACGGCTTTTCCATCCCGCGCAGCTCATAGTCCACGCCCGCTGCCTTGATGCGCGGCACCAGCTGCTGTGAGGCCATGATCAGGTC  >  NC_002947/3287150‑3287419

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