Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 152,157 A→T 14.8% intergenic (‑51/+43) PP_0143 ← / ← PP_0144 membrane protein/insulinase family metalloprotease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947152,1570AT14.8% 57.7 / 4.0 27intergenic (‑51/+43)PP_0143/PP_0144membrane protein/insulinase family metalloprotease
Reads supporting (aligned to +/- strand):  ref base A (12/11);  new base T (2/2);  total (14/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 = 2.76e-01

CCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCC  >  NC_002947/152012‑152293
                                                                                                                                                 |                                                                                                                                        
ccGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATCTGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCaa                                                                                                                                         <  1:157391/147‑1 (MQ=255)
ccGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGTGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCaa                                                                                                                                         <  1:236802/147‑1 (MQ=255)
ccGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGc                                                                                                                                       >  1:31479/1‑149 (MQ=255)
        ggAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGAGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                 >  7:139274/1‑147 (MQ=255)
            agGAAACCGCCGAGCGCAATGATGCCGGGGTAGCGCTGGAGTAGTGCGATGAGGCGCTGGCCGATATGGGGCATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGt                                                                                                                           <  1:79249/149‑1 (MQ=255)
                     ccgAGTGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGTCGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTgg                                                                                                                  <  6:105914/149‑1 (MQ=255)
                     ccgAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTgg                                                                                                                  <  2:31479/149‑1 (MQ=255)
                          cgcACTGCTGCCCGGGGAGCGCTTCAGTAGTTCGATGATGCCCTGGCCGATATGGGGGATTTGGGGCCCTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCaaaa                                                                                                             <  2:76209/149‑1 (MQ=38)
                             aaTGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGa                                                                                                          >  5:17509/1‑149 (MQ=255)
                              atgatgCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGAt                                                                                                         <  6:17509/149‑1 (MQ=255)
                                     ccGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGGGAAAGCGTCGGTACTGGCAAAAGGATAccgcc                                                                                                   >  5:54552/1‑148 (MQ=255)
                                      cGGGTAGCGCTTCAGTAGTTCGATGATGCGCTTTACGATATGGCGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGa                                                                                                 <  2:29267/149‑1 (MQ=255)
                                                               ttgCGGTGGACGTTTTGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTTCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcaccac                                                                         <  6:54552/147‑1 (MQ=255)
                                                                     tGGACGAGATGGGGGATTTGGGGCATTGGGGGCTACGCCAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGGACTGGCAAAAGGAGACCGCCGAGTTACCGCTGCCGTGTAGCACCACTGGGCCg                                                                  >  8:218792/1‑149 (MQ=255)
                                                                         cGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGtt                                                                >  5:317792/1‑147 (MQ=255)
                                                                              tGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcac                                                                            >  3:177386/1‑130 (MQ=255)
                                                                              tGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGcac                                                                            <  4:177386/130‑1 (MQ=255)
                                                                                        ggggCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATgccgccgc                                               >  2:133522/1‑149 (MQ=255)
                                                                                                      tACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGccc                                 >  6:98978/1‑149 (MQ=255)
                                                                                                                  gTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGcc                     >  3:163038/1‑149 (MQ=255)
                                                                                                                      ggTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGc                 <  5:98978/149‑1 (MQ=255)
                                                                                                                                ccAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGc        >  7:98484/1‑148 (MQ=255)
                                                                                                                                 cAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGc        >  2:175503/1‑147 (MQ=255)
                                                                                                                                  aGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAg     <  4:163038/149‑1 (MQ=255)
                                                                                                                                     ctctTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGcc  <  6:317792/149‑1 (MQ=255)
                                                                                                                                      ctttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGc                                   <  2:27827/113‑1 (MQ=255)
                                                                                                                                      ctttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGc                                   >  1:27827/3‑115 (MQ=255)
                                                                                                                                       tttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGa          <  2:310726/138‑1 (MQ=255)
                                                                                                                                       tttCGCGGGCTTGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGa          >  1:310726/2‑139 (MQ=255)
                                                                                                                                                 |                                                                                                                                        
CCGATACCGGAGAGGAAACCGCCGAGCGCAATGATGCCCGGGTAGCGCTTGAGTAGTTCGATGATGCGCTGGACGATATGGGGGATTTGGGGCATTGGGGGCTACGACAGAAAAGTGAGGTGGCTGTGCCAGCCTCTTCGCGGGCAAGCCCGCGAAAGGGTCGGTACTGGCAAAAGGATACCGCCGATTTACCGCTGCCGTGTAGCACCACTGGGCCGTTCAGGTCCATGCCGCCGCCAGCGCAACAGCCCCGCCGCCAGGCCCAGCCCCACGAGCAAGGCC  >  NC_002947/152012‑152293

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