Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,799,179 Δ1 bp intergenic (‑237/+47) pheS ← / ← pheM phenylalanine‑‑tRNA ligase subunit alpha/pheST‑ihfA operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,799,1770T.100.0% 71.8 / NA 18intergenic (‑235/+49)pheS/pheMphenylalanine‑‑tRNA ligase subunit alpha/pheST‑ihfA operon leader peptide
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (9/9);  total (9/9)

GCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCG  >  NC_000913/1799102‑1799252
                                                                           |                                                                           
gCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCttcagg                                                                        >  1:324784/1‑77 (MQ=255)
         acacCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGttt                                                               <  2:324766/80‑1 (MQ=255)
          cacCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTc                                                             <  2:324767/81‑1 (MQ=255)
           acCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTcc                                                            <  2:324768/81‑1 (MQ=255)
                    gaCAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTttctctt                                                   >  1:324790/1‑81 (MQ=255)
                     aCAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTttctcttc                                                  >  2:324799/1‑81 (MQ=255)
                      cAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTttctcttct                                                 <  1:324769/81‑1 (MQ=255)
                       aGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTttctcttctc                                                <  2:324770/81‑1 (MQ=255)
                           tttttATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCAcgc                                            >  1:324800/1‑81 (MQ=255)
                            ttttATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCAcgcg                                           >  2:371357/1‑81 (MQ=255)
                                 tttttCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGc                                      <  2:324771/81‑1 (MQ=255)
                                     tCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGcctc                                   <  2:324772/80‑1 (MQ=255)
                                      cAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGcctcct                                 >  1:324798/1‑81 (MQ=255)
                                         aaCGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGa                              >  2:324794/1‑81 (MQ=255)
                                             cgcATACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGAtt                            >  2:324792/1‑79 (MQ=255)
                                                aTACAAAAAAAGCCTCCACTGGGAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGt                       >  2:324805/1‑81 (MQ=255)
                                                                     ggAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGc   <  1:324774/80‑1 (MQ=255)
                                                                     ggAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGc   <  1:324773/80‑1 (MQ=255)
                                                                      gAGGCTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCg  >  2:324801/1‑80 (MQ=255)
                                                                           |                                                                           
GCGACAATTACACCTGATGAGACAGGCTTTTTATTTTTCAAAACGCGCATACAAAAAAAGCCTCCACTGGGAGGCTTTCAGGCGCTGTTTTCCGTTTCTCTTCTCACGCGCTAGCCTCCTGGATTCAGGTGCTAAAGTAAAAAAAGAAGCG  >  NC_000913/1799102‑1799252

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