Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 253,084 A→G T126A (ACA→GCA)  yafP → GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913253,0840AG100.0% 120.3 / NA 35T126A (ACA→GCA) yafPGNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (17/18);  total (17/18)

ACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCC  >  NC_000913/252959‑253187
                                                                                                                             |                                                                                                       
aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAGCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                            <  2:698756/139‑1 (MQ=255)
aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                            >  1:1050402/1‑139 (MQ=255)
aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                            >  2:447691/1‑139 (MQ=255)
           acCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAAt                                                                                 >  1:220726/1‑139 (MQ=255)
             ccgccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGc                                                                               >  1:1038039/1‑139 (MQ=255)
              cgccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGcc                                                                              <  2:814236/139‑1 (MQ=255)
                ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAg                                                                                                      >  2:422834/1‑113 (MQ=255)
                ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAg                                                                                                      <  1:422834/113‑1 (MQ=255)
                ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                            >  1:989767/1‑139 (MQ=255)
                    ggggTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAg                                                                        <  2:1050402/139‑1 (MQ=255)
                       gTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGt                                                                     >  1:537279/1‑139 (MQ=255)
                                  tttGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGa                                                                                   <  1:796484/114‑1 (MQ=255)
                                  tttGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGa                                                                                   >  2:796484/1‑114 (MQ=255)
                                          aaCCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTatat                                                  <  2:990943/139‑1 (MQ=255)
                                                       gTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                     >  1:137601/1‑91 (MQ=255)
                                                       gTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                     <  2:137601/91‑1 (MQ=255)
                                                            aaTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGata                                            >  1:624408/1‑127 (MQ=255)
                                                            aaTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGata                                            <  2:624408/127‑1 (MQ=255)
                                                                       aCGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATg                                                                                <  1:135119/80‑1 (MQ=255)
                                                                       aCGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATg                                                                                >  2:135119/1‑80 (MQ=255)
                                                                             gACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAg                                                                                                      <  2:994656/52‑1 (MQ=255)
                                                                             gACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAg                                                                                                      >  1:994656/1‑52 (MQ=255)
                                                                             gACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTGAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCaa               >  1:773537/1‑139 (MQ=255)
                                                                                       tAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                            <  2:801243/52‑1 (MQ=255)
                                                                                       tAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                            >  1:801243/1‑52 (MQ=255)
                                                                                        aaCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcg    <  2:1075049/139‑1 (MQ=255)
                                                                                        aaCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcg    <  2:1038039/139‑1 (MQ=255)
                                                                                         aCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTatat                                                  >  2:343964/1‑92 (MQ=255)
                                                                                         aCCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTatat                                                  <  1:343964/92‑1 (MQ=255)
                                                                                          ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc  <  2:220726/139‑1 (MQ=255)
                                                                                          ccGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc  <  2:537279/139‑1 (MQ=255)
                                                                                                            aaCGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGa                                                                                   <  1:866979/40‑1 (MQ=255)
                                                                                                            aaCGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGa                                                                                   >  2:866979/1‑40 (MQ=255)
                                                                                                                 tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCg                                               <  2:417363/71‑1 (MQ=255)
                                                                                                                 tATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTAAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCg                                               >  1:417363/1‑71 (MQ=255)
                                                                                                                             |                                                                                                       
ACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCC  >  NC_000913/252959‑253187

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