Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,273,070 C→G 20.0% intergenic (+35/+256) PP_2003 → / ← PP_2004 hypothetical protein/AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,273,0700CG20.0% 49.2 / 8.7 25intergenic (+35/+256)PP_2003/PP_2004hypothetical protein/AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (9/11);  new base G (2/3);  total (11/14)
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.26e-01

GTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAG  >  NC_002947/2272928‑2273192
                                                                                                                                                |                                                                                                                          
gtggATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGaaa                                                                                                                      <  4:135502/149‑1 (MQ=255)
    aTGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGCCCCATTCGCCGGCAGCCCGCGCAAGGGc                                                                                                                  >  2:105006/1‑149 (MQ=255)
        cAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGCCTCGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAGGGTCg                                                                                                                <  1:212771/149‑1 (MQ=38)
         aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAg                                                                                                                     <  8:186278/141‑1 (MQ=255)
         aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGGCCCTTTCGCGGGCTGCCGGCGAAAg                                                                                                                     >  7:186278/1‑141 (MQ=255)
         aGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGtt                                                                                                             <  7:83841/149‑1 (MQ=255)
           ccAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGCGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGCCCCTTTCGCCGGCAGCCCGCGCAAGGGGCGGTTcc                                                                                                           >  3:200447/1‑148 (MQ=255)
             aGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGACCCCTCCCTTTGATGGGTCGCCGGAACCCCCCCTTTTGCCGGCCGCCCGCGGAAGGGGCGGTTCAgg                                                                                                         >  4:22278/1‑149 (MQ=255)
                    gCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCCGCCCTTTAGCCCGCAGCCCGCGAAAGGGCCGGGTATGGCTGTaa                                                                                                   >  1:190532/1‑148 (MQ=255)
                        ccTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATc                                                                                              <  1:105006/149‑1 (MQ=255)
                            tgGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATcccc                                                                                           <  5:100439/148‑1 (MQ=255)
                                              cTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGAAGGGTCGCCTGAACCCACCCTTTCGCCCCCACCCCCCCAAAGCGCCCGTTAAGCCTGTAAACAACCCCTGTTACAAAACCCCcgcg                                                                         >  8:239093/1‑148 (MQ=255)
                                               tGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGa                                                                        >  8:235780/1‑148 (MQ=255)
                                                  aCCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGaaa                                                                                                                      >  5:177875/1‑99 (MQ=255)
                                                  aCCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGaaa                                                                                                                      <  6:177875/99‑1 (MQ=255)
                                                   cccATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCGGGCTGCCGGCGaaa                                                                                                                      <  6:149962/98‑1 (MQ=255)
                                                   cccATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCGGGCTGCCGGCGaaa                                                                                                                      >  5:149962/1‑98 (MQ=255)
                                                    ccATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCATAACCCACCCAATACCCCGCACCCCGCCAAAAGGCCGGTTAAGGCTGTAAACATACCCAGTTACAAAACCCCCCCGATAAcat                                                                  >  6:174857/1‑147 (MQ=255)
                                                            gTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCGCCCTTTGATGGGTCGCCTGAACCGCCCCATTGGCCGGCCGCCCGCGCAAGGGCCGGGTCAGGCTGGAAACAACCCCGGGTACAGACGCCCCGCGCGAACtttgagacaa                                                          >  5:103666/1‑142 (MQ=255)
                                                               cagcagGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCt                                                       <  7:186097/149‑1 (MQ=255)
                                                                    ggATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTAc                                                  <  7:235780/149‑1 (MQ=255)
                                                                              ccggccgAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACt                                        <  5:265059/149‑1 (MQ=255)
                                                                                     aGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGc                                 <  2:213298/149‑1 (MQ=255)
                                                                                                                  tCGCCTGAAACGACCCTTTCGCCTGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCTATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGc      <  7:239093/149‑1 (MQ=255)
                                                                                                                      cTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAg  <  4:200447/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                          
GTGGATGACAGCCAGGGTTCGCAGCCTGTGGCCCCGCCGGCCGAAGCTGCACCCATTCAGGTGCAGCAGGATGCCCAGCCGGCCGAGGCCGC‑‑CGCCCCCTCCCTTTGATGGGTCGCCTGAACCGACCCTTTCGCCGGCAGCCCGCGAAAGGGCCGGTTCAGGCTGTAAACATCCCCTGTTACAGACGCCCCGCGATAACTTTCTGAACACCTGCTACGTTTACCACTCATAAGCTTGTCATCATCCCTTCATCCCACCTGCTTAG  >  NC_002947/2272928‑2273192

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: