Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,196,790 A→T intergenic (+417/+77) icd → / ← ymfD isocitrate dehydrogenase/e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,7900AT100.0% 65.8 / NA 19intergenic (+417/+77)icd/ymfDisocitrate dehydrogenase/e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (11/8);  total (11/8)

CCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGT  >  NC_000913/1196717‑1196865
                                                                         |                                                                           
cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAgg                                                                      >  1:215270/1‑81 (MQ=255)
  cGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGcc                                                                    >  1:215263/1‑81 (MQ=255)
   gggACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCcg                                                                   >  2:215265/1‑81 (MQ=255)
          gCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAgcgc                                                            <  1:215242/81‑1 (MQ=255)
              ggggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGcc                                                        >  2:215262/1‑81 (MQ=255)
                ggACAGCCAACAGGCATTGGGTGCAATCACCATAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGcccc                                                      <  2:215243/81‑1 (MQ=255)
                 gACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCt                                                     >  1:215271/1‑81 (MQ=255)
                          cAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCacataca                                            >  1:215268/1‑81 (MQ=255)
                            ggCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATc                                          <  2:215244/81‑1 (MQ=255)
                              cATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCtt                                        <  2:215245/81‑1 (MQ=255)
                                      gCAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACt                                 >  2:215266/1‑80 (MQ=255)
                                       cAATCACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACt                                 >  1:215260/1‑79 (MQ=255)
                                           cACCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACt                           >  2:215267/1‑81 (MQ=255)
                                            aCCTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTg                          <  2:215246/81‑1 (MQ=255)
                                              cTTAGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGtt                        <  2:215247/81‑1 (MQ=255)
                                                 aGCGTTCAGGTACATGCGGAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGttt                       <  2:215248/79‑1 (MQ=255)
                                                                  ggAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATc    >  2:215269/1‑81 (MQ=255)
                                                                   gAATGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg   >  2:215274/1‑81 (MQ=255)
                                                                    aaTGTTAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGt  <  1:215249/81‑1 (MQ=255)
                                                                         |                                                                           
CCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGT  >  NC_000913/1196717‑1196865

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