Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,678,604	G→T	A60S (GCA→TCA)	ydgH →	DUF1471 domain‑containing protein YdgH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,678,604	0	G	T	100.0%	166.0 / NA	48	A60S (GCA→TCA)	ydgH	DUF1471 domain‑containing protein YdgH
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (22/26);  total (22/26)

CTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAA  >  NC_000913/1678470‑1678740                                                                                                                                       |                                                                                                                                         cTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGAt                                                                                                                                     <  1:902376/140‑1 (MQ=255) cTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGAt                                                                                                                                     <  1:50992/140‑1 (MQ=255)         ttCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAgg                                                                                                                             >  2:42232/1‑140 (MQ=255)         ttCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAgg                                                                                                                             >  1:121164/1‑140 (MQ=255)                ttAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgcc                                                                                                                      <  2:668990/140‑1 (MQ=255)                ttAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgcc                                                                                                                      >  2:217673/1‑140 (MQ=255)                ttAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTgccgcc                                                                                                                      <  2:723776/140‑1 (MQ=255)                    cGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCtt                                                                                                                  >  1:241547/1‑140 (MQ=255)                       agcaACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCttttt                                                                                                               <  1:848931/140‑1 (MQ=255)                       agcaACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCttttt                                                                                                               <  2:770123/140‑1 (MQ=255)                         caACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTAt                                                                                                             <  2:241604/140‑1 (MQ=255)                         caACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTAt                                                                                                             <  1:842146/140‑1 (MQ=255)                          aaCAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTAtg                                                                                                            <  1:881671/140‑1 (MQ=255)                                  tGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGac                                                                                                    <  2:224111/140‑1 (MQ=255)                                    acacCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGATGTCGacac                                                                                                  >  2:89194/1‑140 (MQ=255)                                         ggAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTg                                                                                             <  2:554327/140‑1 (MQ=255)                                         ggAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTg                                                                                             >  2:657038/1‑140 (MQ=255)                                         ggAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTg                                                                                             >  2:763824/1‑140 (MQ=255)                                         ggAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTg                                                                                             <  2:825729/140‑1 (MQ=255)                                          gAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGa                                                                                            >  1:241336/1‑140 (MQ=255)                                          gAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGa                                                                                            >  2:718726/1‑140 (MQ=255)                                                 cGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGt                                                                                     <  2:623833/140‑1 (MQ=255)                                                 cGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGt                                                                                     <  1:619937/140‑1 (MQ=255)                                                 cGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTCAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGt                                                                                     >  1:242490/1‑140 (MQ=255)                                                                  tttGACCGTGTAGAGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCg                                                                    >  1:836939/1‑140 (MQ=255)                                                                     gACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCgtgt                                                                 <  1:763824/140‑1 (MQ=255)                                                                          tgtAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCg                                                            <  2:121164/140‑1 (MQ=255)                                                                            tAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCt                                                          <  1:42232/140‑1 (MQ=255)                                                                            tAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCt                                                          <  1:657038/140‑1 (MQ=255)                                                                               tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       <  2:241547/140‑1 (MQ=255)                                                                               tGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGATGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       <  1:89194/140‑1 (MQ=255)                                                                                 gTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACct                                                     >  1:721613/1‑140 (MQ=255)                                                                                    aCCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       >  1:381189/1‑135 (MQ=255)                                                                                    aCCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGAc                                                       <  2:381189/135‑1 (MQ=255)                                                                                         tCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCt                                                          >  2:776743/1‑127 (MQ=255)                                                                                         tCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCt                                                          <  1:776743/127‑1 (MQ=255)                                                                                           gTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGcc                                           >  1:91413/1‑140 (MQ=255)                                                                                           gTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGcc                                           >  1:596454/1‑140 (MQ=255)                                                                                                  tGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTg                                    >  1:413548/1‑140 (MQ=255)                                                                                                  tGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTg                                    >  2:195929/1‑140 (MQ=255)                                                                                                       ttGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGaaaaa                               >  1:892161/1‑140 (MQ=255)                                                                                                        tGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAg                              <  1:215665/140‑1 (MQ=255)                                                                                                         ggCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGa                                   >  2:150985/1‑134 (MQ=255)                                                                                                         ggCGAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGa                                   <  1:150985/134‑1 (MQ=255)                                                                                                          gCGAAGCGGTGAAAGCCGTTTCTGGTCGCTCAGATAAAGAAGGTGCCGCCTCGTTTTATGGTGGCGACACGTCTGATTGTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCa                            <  2:422565/140‑1 (MQ=255)                                                                                                            gAAGCGGTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCaga                          <  1:868414/140‑1 (MQ=255)                                                                                                                 ggTGAAAGCCGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGaaa                     >  2:363113/1‑140 (MQ=255)                                                                                                                          cGTTTCTCGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATc            >  1:862579/1‑140 (MQ=255)                                                                                                                              tctcGTCGCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGt        <  2:241336/140‑1 (MQ=255)                                                                                                                                   tcgCTCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTa   >  1:167042/1‑140 (MQ=255)                                                                                                                                    cgctCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTaa  <  1:751964/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                         CTATGGCATTCTCCGTTAACGCAGCAACAGAACTGACACCGGAGCAAGCGGCAGCGGTTAAACCTTTTGACCGTGTAGTGGTTACCGGTCGTTTTAATGCTATTGGCGAAGCGGTGAAAGCCGTTTCTCGTCGCGCAGATAAAGAAGGTGCCGCCTCTTTTTATGTTGTCGACACTTCTGATTTTGGTAACAGCGGTAACTGGCGTGTGGTCGCTGACCTCTATAAAGCCGATGCTGAAAAAGCAGAAGAAACAAGTAATCGCGTAATTAA  >  NC_000913/1678470‑1678740

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