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	98.2%	183.1 / ‑3.4	54	coding (849/1221 nt)	sufS	cysteine desulfurase, stimulated by SufE, selenocysteine lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (29/24);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.35e-01

GGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAG  >  NC_000913/1759541‑1759801                                                                                                                                        |                                                                                                                                gggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                                  <  2:326961/135‑1 (MQ=255) gggCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAatgatg                                                                                                                                  >  1:326961/1‑135 (MQ=255)             gTGAGATCCGGTACAGATTCCAGCTGTGATAACGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGgt                                                                                                                 >  2:134895/1‑139 (MQ=255)              tGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGgtg                                                                                                                >  1:1357041/1‑139 (MQ=255)                      ggTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCt                                                                                                        >  1:441407/1‑139 (MQ=255)                        tACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTc                                                                                                      >  2:653056/1‑139 (MQ=255)                        tACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTc                                                                                                      >  1:252123/1‑139 (MQ=255)                            gATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAAc                                                                                                  <  1:134539/139‑1 (MQ=255)                               tCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGc                                                                                               >  2:93689/1‑139 (MQ=255)                                 cAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCa                                                                                             <  2:934806/139‑1 (MQ=255)                                        gATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                      <  2:252123/139‑1 (MQ=255)                                         aTAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCtt                                                                                     <  1:93689/139‑1 (MQ=255)                                             cgcATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGt                                                                                 >  1:181874/1‑139 (MQ=255)                                             cgcATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccc‑c                                                                                                                            >  2:656894/1‑95 (MQ=255)                                             cgcATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccc‑c                                                                                                                            <  1:656894/95‑1 (MQ=255)                                                     gCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCCGCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAg                                                                         <  1:134895/139‑1 (MQ=255)                                                         cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAAc                                                                                                  <  1:1360026/110‑1 (MQ=255)                                                         cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAAc                                                                                                  >  2:1360026/1‑110 (MQ=255)                                                         cAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGcc                                                                     <  2:535382/139‑1 (MQ=255)                                                          cggcgcTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                      <  2:47450/116‑1 (MQ=255)                                                          cggcgcTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                      >  1:47450/6‑121 (MQ=255)                                                                  ttCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTa                                                                          >  2:532406/1‑125 (MQ=255)                                                                  ttCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTa                                                                          <  1:532406/125‑1 (MQ=255)                                                                  ttCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAg                                                            >  2:390551/1‑139 (MQ=255)                                                                     aTACTCGGCTATGTTATTAAGCCCCAGCACCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCt                                                         >  1:301434/1‑139 (MQ=255)                                                                            gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                      >  1:252295/1‑103 (MQ=255)                                                                            gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                      <  2:252295/103‑1 (MQ=255)                                                                            gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCt                                                         <  1:751446/132‑1 (MQ=255)                                                                            gCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCt                                                         >  2:751446/1‑132 (MQ=255)                                                                              tATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                <  1:771029/139‑1 (MQ=255)                                                                              tATGTTATTAAGCCCCAGCACCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGc                                                <  2:301434/139‑1 (MQ=255)                                                                                tGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGa                                              >  1:851084/1‑139 (MQ=255)                                                                                       aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTAtt                                                                                                                     >  2:669718/1‑61 (MQ=255)                                                                                       aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTAtt                                                                                                                     <  1:669718/61‑1 (MQ=255)                                                                                       aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTgg                                                                                                                   <  2:631919/63‑1 (MQ=255)                                                                                       aaGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTgg                                                                                                                   >  1:631919/1‑63 (MQ=255)                                                                                          ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAAcc                                    >  2:1017784/1‑139 (MQ=255)                                                                                            ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAAc                                                                                                  >  2:213179/1‑75 (MQ=255)                                                                                            ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAAc                                                                                                  <  1:213179/75‑1 (MQ=255)                                                                                            ccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGc                                  >  1:1360901/1‑139 (MQ=255)                                                                                             cAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGt                                                                           <  1:375825/97‑1 (MQ=255)                                                                                             cAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGt                                                                           >  2:375825/1‑97 (MQ=255)                                                                                             cAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGcc                                 >  1:159916/1‑139 (MQ=255)                                                                                              aGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGccc                                >  2:40202/1‑139 (MQ=255)                                                                                              aGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccc‑ccc                                                                                                                          <  1:1052969/48‑2 (MQ=255)                                                                                              aGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑ccc‑ccc                                                                                                                          >  2:1052969/1‑47 (MQ=255)                                                                                                          tACTCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCAcggc                    >  1:649905/1‑139 (MQ=255)                                                                                                             tCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTc                                                                                <  1:419695/76‑1 (MQ=255)                                                                                                             tCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTc                                                                                >  2:419695/1‑76 (MQ=255)                                                                                                             tCCAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACggcggc                 <  1:653056/139‑1 (MQ=255)                                                                                                              ccAGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGTCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCa                <  1:1017784/139‑1 (MQ=255)                                                                                                                aGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTAc                                                                                                             >  1:1001699/1‑44 (MQ=255)                                                                                                                aGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTAc                                                                                                             <  2:1001699/44‑1 (MQ=255)                                                                                                                aGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCa                                                                                             >  1:254576/1‑60 (MQ=255)                                                                                                                aGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCa                                                                                             <  2:254576/60‑1 (MQ=255)                                                                                                                aGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                      >  1:261139/1‑67 (MQ=255)                                                                                                                aGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCt                                                                                      <  2:261139/67‑1 (MQ=255)                                                                                                                aGCGCCGCGCCAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCAtc              <  2:181874/139‑1 (MQ=255)                                                                                                                       cgccAAGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGca       >  1:1331636/1‑139 (MQ=255)                                                                                                                            aGACCAATGATGCCCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCTCGGCGGCATCTCCTGCAACAAg  <  2:851084/139‑1 (MQ=255)                                                                                                                                        |                                                                                                                                GGGCCATAGAGAGTGAGATCCGGTACAGATTCCAGCTGTGATAGCGCATAATGCATCAGATTCTGTTCATACTCGGCTATGTTATTAAGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATG‑CCC‑CCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAG  >  NC_000913/1759541‑1759801

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