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%	158.3 / NA	44	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 (25/19);  total (25/19)

GGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCT  >  NC_000913/1759541‑1759807                                                                                                                                        |                                                                                                                                     gggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                                  >  2:971364/1‑139 (MQ=255)        agagagTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTAt                                                                                                                            <  1:106014/139‑1 (MQ=255)                  aTCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTAcc                                                                                                                  >  1:36684/1‑139 (MQ=255)                        tACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTc                                                                                                            >  1:1192381/1‑139 (MQ=255)                        tACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTc                                                                                                            >  1:294312/1‑139 (MQ=255)                           aGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCaaa                                                                                                         >  2:377733/1‑139 (MQ=255)                                 cAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCa                                                                                                   >  2:1069220/1‑139 (MQ=255)                                       tgATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGc                                                                                             >  1:167667/1‑139 (MQ=255)                                        gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            <  1:1069220/139‑1 (MQ=255)                                        gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            <  1:694740/139‑1 (MQ=255)                                        gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            <  1:377733/139‑1 (MQ=255)                                         aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccccggta                                                                                                                            >  1:491832/1‑100 (MQ=255)                                         aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccccggta                                                                                                                            <  2:491832/105‑6 (MQ=255)                                                   aTGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                                  <  2:1098751/89‑1 (MQ=255)                                                   aTGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccc                                                                                                                                  >  1:1098751/1‑89 (MQ=255)                                                    tGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            <  2:685469/127‑1 (MQ=255)                                                    tGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            >  1:685469/1‑127 (MQ=255)                                                    tGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTa                                                                                >  2:141544/1‑139 (MQ=255)                                                     gCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                               <  1:141544/139‑1 (MQ=255)                                                            aTTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTc                                                                        >  1:678822/1‑139 (MQ=255)                                                                           ggCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                     <  1:1023163/127‑1 (MQ=255)                                                                           ggCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACt                                                                     >  2:1023163/1‑127 (MQ=255)                                                                            gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTg                                                        >  1:500271/1‑139 (MQ=255)                                                                            gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccccggta                                                                                                                            >  2:20086/1‑65 (MQ=255)                                                                            gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccccccggta                                                                                                                            <  1:20086/70‑6 (MQ=255)                                                                              tATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                      <  1:971364/139‑1 (MQ=255)                                                                                tGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGa                                                    >  2:1212610/1‑139 (MQ=255)                                                                                     ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATa                                               <  1:598370/139‑1 (MQ=255)                                                                                     ttaAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATa                                               >  1:1214818/1‑139 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                                       >  1:1277627/1‑46 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                                       <  2:1277627/46‑1 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            >  2:921707/1‑89 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            <  1:921707/89‑1 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGt                                                         <  1:597981/124‑1 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGt                                                         >  2:597981/1‑124 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                      <  1:1066357/127‑1 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                      >  2:1066357/1‑127 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                          >  1:485481/1‑139 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                          <  1:867738/139‑1 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccccc                                                                                                                                >  2:787266/1‑51 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑cccccc                                                                                                                                <  1:787266/52‑2 (MQ=255)                                                                                              aGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGATCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGccc                                      >  2:525533/1‑139 (MQ=255)                                                                                               gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCttt                                                                                          >  1:911623/1‑86 (MQ=255)                                                                                               gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCttt                                                                                          <  2:911623/86‑1 (MQ=255)                                                                                               gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            >  2:658703/1‑84 (MQ=255)                                                                                               gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            <  1:658703/84‑1 (MQ=255)                                                                                                                  cgccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATg                                                                                                 >  1:137258/1‑60 (MQ=255)                                                                                                                  cgccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATg                                                                                                 <  2:137258/60‑1 (MQ=255)                                                                                                                    ccgcgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCt                >  1:921464/1‑139 (MQ=255)                                                                                                                       cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            >  1:260724/1‑60 (MQ=255)                                                                                                                       cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                            <  2:260724/60‑1 (MQ=255)                                                                                                                       cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGca             >  2:629322/1‑139 (MQ=255)                                                                                                                       cgccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGca             >  2:1050182/1‑139 (MQ=255)                                                                                                                         ccAAGACCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGcaac           <  2:1214818/139‑1 (MQ=255)                                                                                                                              aCCAATGATGCCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGc      <  2:500271/139‑1 (MQ=255)                                                                                                                                  atgatgcCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGcttct  <  2:294312/139‑1 (MQ=255)                                                                                                                                  atgatgcCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGcttct  <  1:629322/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                     GGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCT  >  NC_000913/1759541‑1759807

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: ‑