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% 26.7 / 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)

GCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGC  >  NC_000913/357473‑357713
                                                                                                                      |                                                                                                                          
gCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGcc                                                                                                        >  1:772833/1‑139 (MQ=255)
          tGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcgcg                                                                                              <  2:660427/139‑1 (MQ=16)
               gATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCAttt                                                                                         >  1:256300/1‑138 (MQ=2)
                                       gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCca                                                                                                       >  1:855590/1‑101 (MQ=255)
                                       gCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCca                                                                                                       <  2:855590/101‑1 (MQ=255)
                                                                ggCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACcg                                                                                                >  2:1081466/1‑83 (MQ=255)
                                                                      aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCcaaca                                                                                                    <  1:823640/73‑1 (MQ=11)
                                                                      aTAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCcaaca                                                                                                    >  2:823640/1‑73 (MQ=255)
                                                                                     cAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGattca                                                                                                                <  1:869353/46‑1 (MQ=255)
                                                                                     cAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGattca                                                                                                                >  2:869353/1‑46 (MQ=255)
                                                                                     cAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATgcggcg                   >  1:275058/1‑139 (MQ=255)
                                                                                                      tGATCTGTGTCAGCGGTTACCGTGATTCATTTCCGCCAACAACCGCGCATTTATCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGc  <  2:275058/139‑1 (MQ=18)
                                                                                                                      |                                                                                                                          
GCTTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACGTAGAGCACTGAACTCGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTCCAGCCGCTGATCTGTGTCAGCGGCTACCGTGATTCATTCCCGCCAACAACCGCGCATTCCTCCAACGCCATGTGCAAAAATGCCTTCGCAGCGGCTGTCTGCCAGCTGTAGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGC  >  NC_000913/357473‑357713

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: