Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,125,096 A→T 13.8% intergenic (+70/‑112) PP_4511 → / → PP_4512 AraC family transcriptional regulator/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,125,0960AT13.8% 76.2 / 4.5 36intergenic (+70/‑112)PP_4511/PP_4512AraC family transcriptional regulator/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base A (15/16);  new base T (3/2);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCA  >  NC_002947/5124949‑5125239
                                                                                                                                                   |                                                                                                                                               
aGTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTa                                                                                                                                                 >  5:25733/1‑148 (MQ=255)
    ccTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAcc                                                                                                                                             <  6:25733/148‑1 (MQ=255)
       ggtctggtggccagtttggcgaggacgccggtaaggccgaagaacagcgcgccaatgtgcagtgcgccgagggcggtgcggttgttcatggtgggtccttgatggagatctctacagtcagtgccggccctTTCGCGGGTAAACCCGc                                                                                                                                          >  8:234203‑M1/132‑148 (MQ=255)
        ttCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGc                                                                                                                                          <  8:115782/147‑1 (MQ=255)
         tCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTc                                                                                                                                        <  5:210331/148‑1 (MQ=255)
         tCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTc                                                                                                                                        >  3:229017/1‑148 (MQ=255)
          cAGCGGCGATTCTGCGGGCGTTCGGGAGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTTACCCGcgaa                                                                                                                                       <  6:227402/148‑4 (MQ=255)
           aGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGTCGGCCCCTTCGCGGGTTTACCCGcgaag                                                                                                                                      <  7:167009/148‑5 (MQ=255)
           aGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTa                                                                                                                                     >  4:51606/1‑149 (MQ=255)
               gcgATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAgg                                                                                                                                 <  3:51606/149‑1 (MQ=255)
                      tgcGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTAAAACCCGCCAAAAGGGTGcccc                                                                                                                          >  3:44843/1‑149 (MQ=255)
                       gcGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGccccc                                                                                                                         >  1:125179/1‑149 (MQ=255)
                           gCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGGCCCTTCGCGGGTAAAACCCCTCCAACACGATGCCCCCGc                                                                                                                       >  5:233342/1‑147 (MQ=255)
                           gCGTTCGGGGGGACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGGCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCt                                                                                                                     <  4:195824/149‑1 (MQ=255)
                             gTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTAAACCCCACCAAAACGATGGCccccccgg                                                                                                                    >  3:127036/1‑146 (MQ=255)
                                cGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACa                                                                                                                                   >  3:76876/1‑130 (MQ=255)
                                cGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACa                                                                                                                                   <  4:76876/130‑1 (MQ=255)
                                  ggTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAAtgtt                                                                                                              <  2:67511/149‑1 (MQ=255)
                                          ccTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCACCTACAGACTGCCCCCGCCAGAAAGTTGTGCaaaa                                                                                                      >  2:81442/1‑149 (MQ=255)
                                            cgaagaacagcgcgccaatgtgcagtgcgccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCGGAATGTTGTGCAAAAcc                                                                                                    <  1:41662‑M1/55‑1 (MQ=255)
                                              tGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGGCTGAATGTTGTGGaaaaacgg                                                                                                  >  6:216990/1‑145 (MQ=255)
                                                  tttGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAAACCCCTGCAACAAGCTGCCCCCCCCGGAAAGTGGTGCAAAAccggg                                                                                                 >  2:224930/1‑143 (MQ=255)
                                                      cggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCg                                                                                          <  2:125179/149‑1 (MQ=255)
                                                       ggcgAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGTCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCgg                                                                                         <  6:164597/149‑1 (MQ=255)
                                                                  ggCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTTAACACGCACCTACAGGATGCCGCCGCCGGTAAGTTGTGCAAAAACGGTAGGTGCCGGTTATCCcgg                                                                               >  7:184726/1‑147 (MQ=255)
                                                                         ccgagggcggtgcggtggttcatggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGc                                                                       >  2:41662‑M1/66‑149 (MQ=255)
                                                                             aCAGCTGTGTTGCTACATATTGTGTTGGCCTGTGGTTTTACACTTTGCCTTTGCCGGCCCGTTCGCGGGAAAACCCGCTCCTACTGGATGCCCCCGCCTTAATGTTTTGAAAACCCTTTTGGCGGGGTTTTACCCGCGAAAGGGCCgg                                                                    <  6:233342/148‑1 (MQ=255)
                                                                                 cagTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGc                                                                                           >  8:120911/1‑121 (MQ=255)
                                                                                 cagTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGc                                                                                           <  7:120911/121‑1 (MQ=255)
                                                                                               tggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGt                                                                                                 <  7:131664‑M1/58‑1 (MQ=255)
                                                                                               tggtgggtccttgatggagatctctacaggcagtgccggccctTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGt                                                                                                 >  8:131664‑M1/44‑101 (MQ=255)
                                                                                                        gCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGCATGTTGTGCAAAACCTGTTGGAGCGGGTTTACCCGCGGAAGGGCCGGCCCTGCCTGTAGAGATGTCCATCAAGGa                                        >  1:39007/1‑149 (MQ=38)
                                                                                                                      acTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCCCTCCAACAAGAAGCCCCCCCCGGAATTTTGTGCAAAAACTGTAGGATCGGGTTTAACCGCGAAATGGCCCGCACAGCCAGTATAGAAATCCATAAACGACCCACCAAAAaccc                          >  6:250742/1‑148 (MQ=255)
                                                                                                                       atttgcctttgccGGGCCCTTCGCTGGTAAACGCGCTCCTACTGGATTCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAaccac                         <  8:184726/148‑1 (MQ=38)
                                                                                                                            cctttgccGGCCCCTTCGCGGTTAACCCCGCTCCTACAGGATCCCCCCGCCTTATTGTTGTGCAAAACCTTTAGGAGGGGTTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACcaccgca                     <  8:24222/148‑1 (MQ=38)
                                                                                                                                              CGGGTAACCCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAACCCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGca  >  2:39007‑M1/1‑149 (MQ=255)
                                                                                                                                              CGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGca  <  4:186093‑M1/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
AGTGCCTTTTCAGCGGCGATGCTGCGGGCGTTCGGGTGTACACCTTTGGCTTTGCGGCGAGAGTCTGGCGATAACTGACAGCAGTCTTGCGACAGAATGTGTTTGCCTGAGGCTTTACACTTTGCCTTTGCCGGCCCCTTCGCGGGTAAACCCGCTCCTACAGGATGCCCCCGCCTGAATGTTGTGCAAAACCTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGCACTGCCTGTAGAGATCTCCATCAAGGACCCACCATGAACCACCGCACCGCCCTCGGCGCACTGCA  >  NC_002947/5124949‑5125239

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: