Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,393,847	T→C	S804S (TCA→TCG)	yhdP ←	AsmA2 domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,393,847	0	T	C	100.0%	203.9 / NA	56	S804S (TCA→TCG)	yhdP	AsmA2 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (25/31);  total (25/31)

TCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCACTTTACCATCCCACGCCACGCTGCCACTCAATG  >  NC_000913/3393711‑3393968                                                                                                                                         |                                                                                                                          tCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGaa                                                                                                                        <  1:335085/140‑1 (MQ=255)   aGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAgg                                                                                                                      >  2:583424/1‑140 (MQ=255)   aGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAgg                                                                                                                      >  1:510659/1‑140 (MQ=255)   aGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAgg                                                                                                                      >  2:788701/1‑140 (MQ=255)    gCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGt                                                                                                                     >  1:465370/1‑140 (MQ=255)    gCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGt                                                                                                                     >  1:421321/1‑140 (MQ=255)           aCCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGAc                                                                                                              <  1:336849/140‑1 (MQ=255)               aGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCt                                                                                                          >  2:169915/1‑140 (MQ=255)               aGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCt                                                                                                          >  1:587423/1‑140 (MQ=255)               aGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCt                                                                                                          >  2:343770/1‑140 (MQ=255)                gCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTACCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTc                                                                                                         <  1:747386/140‑1 (MQ=255)                  aaCCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTcac                                                                                                       <  2:465370/140‑1 (MQ=255)                  aaCCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTcac                                                                                                       <  2:78250/140‑1 (MQ=255)                   aCCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTcaca                                                                                                      <  1:583424/140‑1 (MQ=255)                         ggCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTcaca                                                                                                      >  2:352687/1‑134 (MQ=255)                         ggCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTcaca                                                                                                      <  1:352687/134‑1 (MQ=255)                           cTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAg                                                                                              >  1:375184/1‑140 (MQ=255)                              ttGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATc                                                                                           >  1:277651/1‑140 (MQ=255)                                    tGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGtt                                                                                     <  2:587423/140‑1 (MQ=255)                                    tGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGtt                                                                                     <  2:669973/140‑1 (MQ=255)                                    tGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTCGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGtt                                                                                     <  2:591101/140‑1 (MQ=255)                                         aTCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGct                                                                                <  1:169915/140‑1 (MQ=255)                                          tCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctc                                                                               <  2:277651/140‑1 (MQ=255)                                          tCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctc                                                                               <  2:433396/140‑1 (MQ=255)                                           ccGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctct                                                                              <  1:21683/140‑1 (MQ=255)                                            cGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTa                                                                             >  2:40533/1‑140 (MQ=255)                                               aCCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctct                                                                              >  1:568518/1‑136 (MQ=255)                                               aCCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctct                                                                              <  2:568518/136‑1 (MQ=255)                                               aCCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGt                                                                          >  1:127235/1‑140 (MQ=255)                                                ccAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGtt                                                                         >  1:780102/1‑140 (MQ=255)                                                 cAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTta                                                                        <  2:223589/140‑1 (MQ=255)                                                    ccTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctct                                                                              <  2:129769/131‑1 (MQ=255)                                                    ccTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGctct                                                                              >  1:129769/1‑131 (MQ=255)                                                         ccGGTTAATTCAAAGCTGTTGAGATTGTCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGc                                                                >  2:601629/1‑140 (MQ=255)                                                             ttAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCa                                                            >  1:278327/1‑140 (MQ=255)                                                              tAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAg                                                           <  1:788701/140‑1 (MQ=255)                                                              tAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAg                                                           <  2:421321/140‑1 (MQ=255)                                                                  tCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATg                                                       >  2:438099/1‑140 (MQ=255)                                                                   cAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGa                                                      <  2:375184/140‑1 (MQ=255)                                                                   cAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGa                                                      <  2:780102/140‑1 (MQ=255)                                                                   cAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGa                                                      <  2:510659/140‑1 (MQ=255)                                                                       gCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGt                                                                   >  1:803363/1‑123 (MQ=255)                                                                       gCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGt                                                                   <  2:803363/123‑1 (MQ=255)                                                                                aTTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGa                                             >  1:809575/1‑136 (MQ=255)                                                                                aTTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGa                                             <  2:809575/136‑1 (MQ=255)                                                                                  tGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATg                                       <  2:127235/140‑1 (MQ=255)                                                                                   gCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGc                                      <  2:278327/140‑1 (MQ=255)                                                                                   gCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGc                                      >  2:859867/1‑140 (MQ=255)                                                                                      aTCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCa                                   <  1:859867/140‑1 (MQ=255)                                                                                      aTCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCa                                   <  1:571116/140‑1 (MQ=255)                                                                                      aTCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCa                                   <  1:438099/140‑1 (MQ=255)                                                                                                     tACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCACTTTACCATCccacg                    >  2:1896/1‑140 (MQ=255)                                                                                                     tACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCACTTTACCATCccacg                    >  1:310639/1‑140 (MQ=255)                                                                                                      aCCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCACTTTACCATCccacgc                   <  1:343770/140‑1 (MQ=255)                                                                                                            aGTGGTTCACCCGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCACTTTACCATCCCACGCCACGCt             <  1:1896/140‑1 (MQ=255)                                                                                                                       cGCAGGTTTGGCTAACGGCGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCACTTTACCATCCCACGCCACGCTGCCACTCAATg  >  1:54913/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                          TCAGCTTTTGACCGAGCAACCAGCGGCTATTGAAATGATTATCCGCACCAGCCTGTCCGGTTAATTCAAAGCTGTTGAGATTGCCATCAACCTTAACGTTTACCGCTAGTGGTTCACCCGCAGGTTTGGCTAACGGTGAAGGTAAGTGACTGCTCACATTCTTTAGATCGCCGTTCAGCTCTATGTTATAGGTCGCACCAGCATGATAAGGCAGATCAATGCCCACTTTACCATCCCACGCCACGCTGCCACTCAATG  >  NC_000913/3393711‑3393968

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