Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,366,953 G→A G10R (GGA→AGA)  yhcF → putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,366,9530GA100.0% 135.7 / NA 39G10R (GGA→AGA) yhcFputative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/17);  total (22/17)

TTAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCGGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATT  >  NC_000913/3366826‑3367082
                                                                                                                               |                                                                                                                                 
ttAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCttt                                                                                                                        <  2:362941/139‑1 (MQ=255)
     cTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGc                                                                                                                   <  2:70764/139‑1 (MQ=255)
        acacTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCAt                                                                                                                >  1:588475/1‑139 (MQ=255)
        acacTGACAATGGCTCCAACTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCAt                                                                                                                >  2:355356/1‑139 (MQ=255)
              aCAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAgc                                                                                                          >  1:402580/1‑139 (MQ=255)
                        cAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCt                                                                                                >  2:312731/1‑139 (MQ=255)
                                     tGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCaa                                                                                   >  2:479758/1‑139 (MQ=255)
                                          gggAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGAc                                                                              <  1:355356/139‑1 (MQ=255)
                                             aaCTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCaa                                                                                                      >  1:230654/1‑112 (MQ=255)
                                             aaCTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCaa                                                                                                      <  2:230654/112‑1 (MQ=255)
                                                 aaCCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAgg                                                                       <  2:588475/139‑1 (MQ=255)
                                                       cATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAAt                                                                                  <  2:117310/122‑1 (MQ=255)
                                                       cATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAAt                                                                                  >  1:117310/1‑122 (MQ=255)
                                                                    aCAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGc                                                    >  2:299101/1‑139 (MQ=255)
                                                                      aGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                    <  1:205232/89‑1 (MQ=255)
                                                                      aGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGc                                                                                                    >  2:205232/1‑89 (MQ=255)
                                                                       ggCAATCTAACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                 >  1:63155/1‑139 (MQ=255)
                                                                          aaTCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGc                                                                                                                   <  1:527830/70‑1 (MQ=255)
                                                                          aaTCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGc                                                                                                                   >  2:527830/1‑70 (MQ=255)
                                                                           aTCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACt                                             >  1:542144/1‑139 (MQ=255)
                                                                            tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTg                                                        <  1:399260/127‑1 (MQ=255)
                                                                            tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTg                                                        >  2:399260/1‑127 (MQ=255)
                                                                            tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                 >  2:261558/1‑134 (MQ=255)
                                                                            tCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTg                                                 <  1:261558/134‑1 (MQ=255)
                                                                                      cTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGc                                                                                                                   <  2:219068/58‑1 (MQ=255)
                                                                                      cTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGc                                                                                                                   >  1:219068/1‑58 (MQ=255)
                                                                                      cTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCa                                                                                                       <  2:510480/70‑1 (MQ=255)
                                                                                      cTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCa                                                                                                       >  1:510480/1‑70 (MQ=255)
                                                                                      cTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTAt                                                                                      >  2:216664/1‑87 (MQ=255)
                                                                                      cTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTAt                                                                                      <  1:216664/87‑1 (MQ=255)
                                                                                              ggCAAAATGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAAcggcg                          >  2:149754/1‑139 (MQ=255)
                                                                                                   aaTGAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAgtt                     >  1:530662/1‑139 (MQ=255)
                                                                                                      gAACAATGTAAAATTACTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAgttgtt                  <  1:299101/139‑1 (MQ=255)
                                                                                                                      cTGATTGCTAGAAGTGCCTTTTTTTCCATGTCAGCTCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTg                                                                     <  2:392223/72‑1 (MQ=255)
                                                                                                                      cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTg                                                                     >  1:392223/1‑72 (MQ=255)
                                                                                                                      cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTa                                               <  2:279614/94‑1 (MQ=255)
                                                                                                                      cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTa                                               >  1:279614/1‑94 (MQ=255)
                                                                                                                      cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCAtt  >  2:587350/1‑139 (MQ=255)
                                                                                                                      cTGATTGCCAGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCAtt  <  1:149754/139‑1 (MQ=255)
                                                                                                                               |                                                                                                                                 
TTAGCCTCACACTGACAATGGCTCCAGCTTTTAAGGATGAAGGGGAACTAACCGACATGACAGATATTACAGGCAATCTGACGGTCCTGGTGGAGGCAAAATGAACAATGTAAAATTACTGATTGCCGGAAGTGCCTTTTTTGCCATGTCAGCGCAAGCCGCTGATAGAGTATCAATTGACGTTAAGGTGACTCTGGAAGCTGCAGCTTGTACTCCAATACTAAGCAACGGCGGAGTTGTTAACTTCGGTAGTCATT  >  NC_000913/3366826‑3367082

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