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% 143.9 / NA 41T126A (ACA→GCA) yafPGNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (21/20);  total (21/20)

GTTGACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACA  >  NC_000913/252955‑253219
                                                                                                                                   |                                                                                                                                       
gttgACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAagc                                                                                                                                >  1:1297712/1‑139 (MQ=255)
   gACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagc                                                                                                                             >  1:677890/1‑139 (MQ=255)
    aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                            >  2:954540/1‑139 (MQ=255)
    aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                            >  2:884647/1‑139 (MQ=255)
    aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                            >  1:131014/1‑139 (MQ=255)
    aCCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                            >  2:30069/1‑139 (MQ=255)
           aTACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                                     <  1:741767/139‑1 (MQ=255)
             acacCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGa                                                                                                                   <  1:30069/139‑1 (MQ=255)
             acacCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGa                                                                                                                   <  2:289103/139‑1 (MQ=255)
                    ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                                            >  2:1141897/1‑139 (MQ=255)
                    ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                                            >  1:15530/1‑139 (MQ=255)
                    ccgTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCgg                                                                                                            >  2:1208304/1‑139 (MQ=255)
                          ggTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGATTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCGGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                                      <  1:1141897/139‑1 (MQ=255)
                          ggTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                                      <  2:1198599/139‑1 (MQ=255)
                          ggTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCg                                                                                                      <  1:954540/139‑1 (MQ=255)
                                 aGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTa                                                                                               <  2:677890/139‑1 (MQ=255)
                                   cgcTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                                     >  2:238659/1‑115 (MQ=255)
                                   cgcTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                                     <  1:238659/115‑1 (MQ=255)
                                                    tGATTAAGTCTGAATCCGAACTTACGGTGGTCGCAAGTATAACCGCAAAACCC‑‑TTTTATGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAgcg                                                                                                                         <  2:871680/94‑1 (MQ=255)
                                                    tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGTATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAgcg                                                                                                                         >  1:871680/1‑94 (MQ=255)
                                                    tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                            <  2:1329273/93‑1 (MQ=255)
                                                    tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCCTTTTTTTTGAACGTTATGGTTTTCAGGCAGTTAagcagca                                                                                                                            >  1:1329273/1‑93 (MQ=255)
                                                    tGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGata                                                                            >  2:1135928/1‑139 (MQ=255)
                                                           gTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCg                                                                     >  2:934436/1‑139 (MQ=255)
                                                             cTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCa                                                                   <  2:1297712/139‑1 (MQ=255)
                                                                aaTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCa                                                                   <  1:342107/136‑1 (MQ=255)
                                                                aaTCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCa                                                                   >  2:342107/1‑136 (MQ=255)
                                                                      aaCTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTtat                                                                                    <  1:185923/113‑1 (MQ=255)
                                                                      aaCTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTtat                                                                                    >  2:185923/1‑113 (MQ=255)
                                                                          tACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAg                                                      >  2:1157502/1‑139 (MQ=255)
                                                                             ggTGGACGCAAACATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCtt                                                   >  1:1024720/1‑139 (MQ=255)
                                                                                  aCGCAAACATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAAt                                              <  2:1024720/139‑1 (MQ=255)
                                                                                            aaCCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcg                                    <  2:15530/139‑1 (MQ=255)
                                                                                              ccGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                                     <  1:570239/56‑1 (MQ=255)
                                                                                              ccGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGtt                                                                                                                     >  2:570239/1‑56 (MQ=255)
                                                                                              ccGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAcgcc                                  <  1:1135928/139‑1 (MQ=255)
                                                                                                  aaaaCCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTatat                                                                                  <  1:828569/87‑1 (MQ=255)
                                                                                                  aaaaCCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTatat                                                                                  >  2:828569/1‑87 (MQ=255)
                                                                                                    aaCCC‑‑TTTTTTGAACGTTATGGTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAAAGCCCGCCTg                            >  2:874239/1‑139 (MQ=255)
                                                                                                                           gTTTTCAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGAtt       <  1:1157502/139‑1 (MQ=255)
                                                                                                                                cAGGCAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCcaca  <  1:874239/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                       
GTTGACCCTGAATACACCCGCCGTGGGGTTGCCAGCGCTTTGTTAAAACCTTTGATTAAGTCTGAATCCGAACTTACGGTGGACGCAAGCATAACCGCAAAACCC‑‑TTTTTTGAACGTTATGGTTTTCAGACAGTTAAGCAGCAGCGCGTTGAATGCCGGGGAGCGTGGTTTACTAATTTTTATATGCGATATAAACCGCAACATTAAATCCAGCTTGCAATGAAAATAACGCCCGCCTGGTATGTGCGCGTTTTCTGATTCCACA  >  NC_000913/252955‑253219

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: