Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,254,369	G→A	W17* (TGG→TGA)	yhaK →	bicupin‑related protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,254,369	0	G	A	100.0%	130.3 / NA	38	W17* (TGG→TGA)	yhaK	bicupin‑related protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (22/16);  total (22/16)

ACCGGGCAGAATATCTGGCTAACAATGTAGCGTCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGGTTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCTTTCAGCCGCGAACC  >  NC_000913/3254243‑3254492                                                                                                                               |                                                                                                                            aCCGGGCAGAATATCTGGCTAACAATGTAGCGTCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAg                                                                                                                       >  1:354140/1‑133 (MQ=255)  ccGGGCAGCATATCTTGCTAACTATGTAGCTTCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAgg                                                                                                                      <  1:197261/133‑1 (MQ=255)   cGGGCAGAATATCTGGCTAACAATGTAGCGTCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTA‑TACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGcc                                                                                                                    >  1:105219/1‑133 (MQ=255)      gCAGAATATCTGGCTAACAATGTAGCGTCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCg                                                                                                                  >  2:82915/1‑133 (MQ=255)                    tAACAATGTAGCGTCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCCAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCt                                                                                                    <  1:149410/133‑1 (MQ=255)                      aCAATGTAGCGTCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTA‑TACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTg                                                                                                 <  2:105219/133‑1 (MQ=255)                                gTCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACactact                                                                                        >  1:378038/1‑133 (MQ=255)                                 tCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGac                                                                                              <  2:403337/126‑1 (MQ=255)                                 tCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGac                                                                                              >  1:403337/1‑126 (MQ=255)                                    aGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGa                                                                                    >  1:82714/1‑133 (MQ=255)                                      ccccTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGAcc                                                                                  >  1:48095/1‑133 (MQ=255)                                         cTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGa                                                                               <  2:140114/133‑1 (MQ=255)                                         cTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGa                                                                               <  2:354140/133‑1 (MQ=255)                                           tACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGaaa                                                                             >  1:135059/1‑133 (MQ=255)                                            aCCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAAt                                                                            >  2:51835/1‑133 (MQ=255)                                                  aaaGGAAGTATAGAGTGGTCAAGAATTACGATAACTACCCGTACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAg                                                                      <  2:174132/133‑1 (MQ=255)                                                         gTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGc                                                               <  1:105750/133‑1 (MQ=255)                                                                 agGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGc                                                       <  1:593684/133‑1 (MQ=255)                                                                    tCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCgtg                                                    >  1:340390/1‑133 (MQ=255)                                                                         aaTTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCtt                                               >  1:464710/1‑133 (MQ=255)                                                                                  tactacCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGaa                                      >  2:323765/1‑133 (MQ=255)                                                                                     tacCCGAAATGCCAGACAATGTTGACAAGCCGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTg                                   <  1:201916/133‑1 (MQ=255)                                                                                         cGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTgg                               >  2:585244/1‑133 (MQ=255)                                                                                            aCTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCTTGCGTGTACTTGACCAGGAAGGGCTGGccc                            >  1:214742/1‑133 (MQ=255)                                                                                            aCTGCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGccc                            >  1:15531/1‑133 (MQ=255)                                                                                               gCCAGGCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAg                         >  2:156148/1‑133 (MQ=255)                                                                                                   ggCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTgc                     <  1:82915/133‑1 (MQ=255)                                                                                                   ggCAGTGTGGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTgc                     <  2:15531/133‑1 (MQ=255)                                                                                                        tgtgGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCt                >  1:43676/1‑133 (MQ=255)                                                                                                        tgtgGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCt                >  2:604258/1‑133 (MQ=255)                                                                                                        tgtgGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCt                >  2:621584/1‑133 (MQ=255)                                                                                                          tgGACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCttt              >  1:162300/1‑133 (MQ=255)                                                                                                            gACAAGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCTTTCa            <  2:48095/133‑1 (MQ=255)                                                                                                               aaGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCTTTCAGcc         <  1:51835/133‑1 (MQ=255)                                                                                                                aGCAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCTTTCAGCcg        >  2:24554/1‑133 (MQ=255)                                                                                                                  cAGACTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCTTTCAGCcgcg      >  2:102943/1‑133 (MQ=255)                                                                                                                   aGACTACGGACGATTTCAGGCCCGTTATTCTTGTTCCTTTGAACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCTTTCAGCCGCGa     <  2:217188/133‑1 (MQ=255)                                                                                                                     aCTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCTTTCAGCCGCGAAc   <  1:102943/133‑1 (MQ=255)                                                                                                                      cTACGGATGATTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCTTTCAGCCGCGAAcc  <  1:585244/133‑1 (MQ=255)                                                                                                                               |                                                                                                                            ACCGGGCAGAATATCTGGCTAACAATGTAGCGTCCAGCCCCTTACCATAAAAGGAAGTAAAGAGAGGTCAAGAATTATGATTACTACCCGAACTGCCAGGCAGTGTGGACAAGCAGACTACGGATGGTTGCAGGCCCGGTATACTTTTTCCTTTGGACACTACTTCGACCCGAAATTGTTAGGCTATGCCTCCCTGCGTGTGCTTAACCAGGAAGTGCTGGCCCCAGGTGCCGCCTTTCAGCCGCGAACC  >  NC_000913/3254243‑3254492

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