Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,580,094	C→T	D233N (GAT→AAT)	ydeM ←	putative YdeN‑specific sulfatase‑maturating enzyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,580,094	0	C	T	100.0%	150.2 / NA	43	D233N (GAT→AAT)	ydeM	putative YdeN‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/24);  total (19/24)

TTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTA  >  NC_000913/1580000‑1580188                                                                                               |                                                                                               ttGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                           <  1:512150/100‑1 (MQ=255) ttGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                           <  2:827393/100‑1 (MQ=255) ttGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGtttt                                                                                           <  1:420905/100‑1 (MQ=255)  tGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGttttt                                                                                          <  2:310244/100‑1 (MQ=255)   gCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTa                                                                                         <  2:1162830/100‑1 (MQ=255)    cAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTaa                                                                                        >  2:536700/1‑100 (MQ=255)      ggACTCCTGGAAAATACAACTGGTATGCAAATAAAAAAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAcc                                                                                      <  2:1376063/100‑1 (MQ=255)       gACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAAccc                                                                                     <  2:1434608/100‑1 (MQ=255)             tGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCa                                                                               <  1:1413390/100‑1 (MQ=255)                 aaaTACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATaaa                                                                           <  1:591301/100‑1 (MQ=255)                    tACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAAt                                                                        >  1:1571243/1‑100 (MQ=255)                      cAACTGGTATGCCCATCCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATgg                                                                      >  2:189861/1‑100 (MQ=255)                            gTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACa                                                                <  1:1351000/100‑1 (MQ=255)                                 cccATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAAc                                                           >  1:1564023/1‑100 (MQ=255)                                   cATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACtt                                                         >  2:625850/1‑100 (MQ=255)                                        cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAATTTGCCAt                                                    >  2:598300/1‑100 (MQ=255)                                        cccAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATCAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCAt                                                    <  2:854953/100‑1 (MQ=255)                                           aaaaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAg                                                 >  1:1014618/1‑100 (MQ=255)                                             aaaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGcc                                               <  1:193215/100‑1 (MQ=255)                                              aaaCGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCg                                              <  2:1283754/100‑1 (MQ=255)                                                 cGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATACACTTGACATaacaagtgg                                           >  1:680686/1‑93 (MQ=255)                                                 cGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTgg                                           >  1:1541323/1‑100 (MQ=255)                                                     tGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAgg                                       <  2:477852/100‑1 (MQ=255)                                                      gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                      <  1:17708/100‑1 (MQ=255)                                                      gACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGc                                      >  1:1353999/1‑100 (MQ=255)                                                           cgatctCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGa                                 <  2:431113/95‑1 (MQ=255)                                                              tttCAAAATGAGGGATGAAAATTTCACCCACATTGCTTTTAACCCATTGCATAAAAATGGATGACGTAAACTTGACATAAGCCGTGGGAGGCACAGacaa                              >  2:828125/1‑100 (MQ=255)                                                               ttCAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaaaaa                             >  2:719193/1‑100 (MQ=255)                                                                 cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGaagatcg                           >  1:431115/1‑95 (MQ=255)                                                                 cAAACTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATc                           <  1:625853/100‑1 (MQ=255)                                                                  aaaCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCa                          <  1:995635/100‑1 (MQ=255)                                                                    aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                        <  1:719665/100‑1 (MQ=255)                                                                    aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                        <  2:111480/100‑1 (MQ=255)                                                                    aCTGACGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAAt                        <  1:435913/100‑1 (MQ=255)                                                                         cGGATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATAc                   <  2:547704/100‑1 (MQ=255)                                                                          ggATGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACg                  >  1:1342009/1‑100 (MQ=255)                                                                             tGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa               >  1:694512/1‑100 (MQ=255)                                                                             tGAAAATTTCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGaa               >  1:694498/1‑100 (MQ=255)                                                                                   tttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGtt         <  1:679956/100‑1 (MQ=255)                                                                                    ttCACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtc        >  2:658201/1‑100 (MQ=255)                                                                                      cACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc      <  2:863995/100‑1 (MQ=255)                                                                                      cACCCACATTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTtctc      >  1:90815/1‑100 (MQ=255)                                                                                          cacaTTGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTa  >  2:301315/1‑100 (MQ=255)                                                                                               |                                                                                               TTGCAGGACTCCTGGAAAATACAACTGGTATGCCCATTCCCCAAAAAACGGCTGACAAAGCTTTCAAACTGACGGATGAAAATTTCACCCACATCGTTTTTAACCCATTGCATAAAAATGGTTGACATAAACTTGCCATAAGCCGTGGGAGGCACAGAAAAATCAATGATACGGAATGTGTTCTCACTA  >  NC_000913/1580000‑1580188

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