Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 357,591 C→T intergenic (+137/+200) codA → / ← cynR cytosine/isoguanine deaminase/transcriptional activator of cyn operon, autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913357,5910CT100.0% 28.2 / NA 12intergenic (+137/+200)codA/cynRcytosine/isoguanine deaminase/transcriptional activator of cyn operon, autorepressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/5);  total (7/5)

CGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCA  >  NC_000913/357486‑357725
                                                                                                         |                                                                                                                                      
cGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCGTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCAt                                                                                                       <  2:114073/139‑1 (MQ=12)
                                                              cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGat                                                                                                                               <  1:106053/53‑1 (MQ=255)
                                                              cGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGat                                                                                                                               >  2:106053/1‑53 (MQ=255)
                                                                   cgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATg                                    >  1:689260/1‑139 (MQ=255)
                                                                   cgcATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATg                                    >  2:26005/1‑139 (MQ=255)
                                                                          ggCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg                             >  1:158087/1‑139 (MQ=17)
                                                                          ggCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTg                             <  2:904262/139‑1 (MQ=17)
                                                                                    gCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTAt                   <  1:972523/139‑1 (MQ=255)
                                                                                         tgttCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGc              >  1:531268/4‑139 (MQ=18)
                                                                                                tgtCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTa       <  2:531268/139‑1 (MQ=18)
                                                                                                     gCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCa  >  1:601098/1‑139 (MQ=255)
                                                                                                      cGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACg         >  2:1070302/1‑131 (MQ=255)
                                                                                                         |                                                                                                                                      
CGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCA  >  NC_000913/357486‑357725

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