Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,759,680	(C)5→6	coding (844/1221 nt)	sufS ←	cysteine desulfurase, stimulated by SufE, selenocysteine lyase, PLP‑dependent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,759,675	1	.	C	100.0%	167.8 / NA	49	G283G (GGC→GGG)	sufS	cysteine desulfurase, stimulated by SufE, selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base C (27/22);  total (27/22)

TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGC  >  NC_000913/1759537‑1759803                                                                                                                                            |                                                                                                                                 ttGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                                   <  2:772700/139‑1 (MQ=255)                            tACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTc                                                                                                        >  2:405909/1‑139 (MQ=255)                                 aTTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAAcc                                                                                                   >  2:1782322/1‑139 (MQ=255)                                 aTTCCAGCTGTGATAGCGCATAACGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGAGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAAcc                                                                                                   >  1:1875006/1‑139 (MQ=255)                                   tCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                              >  1:346542/1‑109 (MQ=255)                                   tCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                              <  2:346542/109‑1 (MQ=255)                                      aGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCACCAAGACCAatgatg                                                                                                                                   <  2:124274/102‑1 (MQ=255)                                      aGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCACCAAGACCAatgatg                                                                                                                                   >  1:124274/1‑102 (MQ=255)                                            gATAGCGCATAATGCATCAGATTCTGTTCATACTCGTCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                        <  2:1181323/139‑1 (MQ=255)                                            gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                        <  1:1782322/139‑1 (MQ=255)                                            gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                        <  2:1875006/139‑1 (MQ=255)                                            gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                        <  1:1283450/139‑1 (MQ=255)                                                     tAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAg                                                                               >  2:1496160/1‑139 (MQ=255)                                                        tGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTa                                                                            >  1:1657748/1‑139 (MQ=255)                                                         gCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGTCGAAACATACTCCAGCGTCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACGGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                           <  2:520578/139‑1 (MQ=255)                                                             cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGcc                                                                       <  1:1496160/139‑1 (MQ=255)                                                                aTTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTc                                                                    >  1:1001102/1‑139 (MQ=255)                                                                  tCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCAc                                                                  >  2:539932/1‑139 (MQ=255)                                                                      ttCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                 >  1:1250344/1‑136 (MQ=255)                                                                      ttCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                 <  2:1250344/136‑1 (MQ=255)                                                                      ttCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAg                                                              >  2:1518417/1‑139 (MQ=255)                                                                           aCTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                        >  1:1714588/1‑108 (MQ=255)                                                                           aCTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                        <  2:1714588/108‑1 (MQ=255)                                                                             tCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACg                                                       <  2:1810357/139‑1 (MQ=255)                                                                               ggCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGt                                                     <  1:1518417/139‑1 (MQ=255)                                                                                  tATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                  <  1:137645/139‑1 (MQ=255)                                                                                         ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                              <  1:555360/55‑1 (MQ=255)                                                                                         ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                              >  2:555360/1‑55 (MQ=255)                                                                                           aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGa                                         >  1:961317/1‑139 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                           >  1:1316112/1‑102 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                           <  2:1316112/102‑1 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                 >  1:1513269/1‑112 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                 <  2:1513269/112‑1 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                      >  2:1327801/1‑139 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                      >  2:681962/1‑139 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                      >  2:1104802/1‑139 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                      >  2:1099746/1‑139 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                      >  1:1367779/1‑139 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccccc                                                                                                                            >  2:179951/1‑51 (MQ=255)                                                                                              ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccccc                                                                                                                            <  1:179951/52‑2 (MQ=255)                                                                                                ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGc                                    >  2:1808492/1‑139 (MQ=255)                                                                                                ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cc                                                                                                                                >  2:1873983/1‑46 (MQ=255)                                                                                                ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cc                                                                                                                                <  1:1873983/46‑1 (MQ=255)                                                                                                 cAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                 <  1:1863782/109‑1 (MQ=255)                                                                                                 cAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                 >  2:1863782/1‑109 (MQ=255)                                                                                                         aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTAtt                                                                                                                       <  2:1708281/47‑1 (MQ=255)                                                                                                         aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTAtt                                                                                                                       >  1:1708281/1‑47 (MQ=255)                                                                                                         aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCttt                                                                                      <  1:1179945/80‑1 (MQ=255)                                                                                                         aaaCATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCttt                                                                                      >  2:1179945/1‑80 (MQ=255)                                                                                                              tACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCAcggc                      >  1:1088639/1‑139 (MQ=255)                                                                                                                  ccAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCa                  <  1:1099746/139‑1 (MQ=255)                                                                                                                     gcgccgcgccAAGACCAATGATGCCCCCCGGTATTgg                                                                                                                     >  1:580464/1‑37 (MQ=255)                                                                                                                     gcgccgcgccAAGACCAATGATGCCCCCCGGTATTgg                                                                                                                     <  2:580464/37‑1 (MQ=255)                                                                                                                        ccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGTCTGACGGTGGc                                                  <  1:1090921/101‑1 (MQ=255)                                                                                                                        ccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                  >  2:1090921/1‑101 (MQ=255)                                                                                                                         cgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                                 <  2:333379/117‑1 (MQ=255)                                                                                                                         cgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcccc                                 >  1:333379/1‑117 (MQ=255)                                                                                                                           cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                           <  2:42193/73‑1 (MQ=255)                                                                                                                           cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                           >  1:42193/1‑73 (MQ=255)                                                                                                                                  aCCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGc  <  1:681962/139‑1 (MQ=255)                                                                                                                                      atgatgcCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                        <  1:917417/49‑1 (MQ=255)                                                                                                                                      atgatgcCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                        >  2:917417/1‑49 (MQ=255)                                                                                                                                       tgatgcCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGcaaca      <  1:550755/130‑1 (MQ=255)                                                                                                                                       tgatgcCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGcaaca      >  2:550755/1‑130 (MQ=255)                                                                                                                                            |                                                                                                                                 TTGTGGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGC  >  NC_000913/1759537‑1759803

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