Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,760,000	T→C	N175S (AAC→AGC)	sufS ←	L‑cysteine desulfurase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,760,000	0	T	C	98.2%	200.9 / ‑4.3	56	N175S (AAC→AGC)	sufS	L‑cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (32/23);  total (33/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

TAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCC  >  NC_000913/1759865‑1760135                                                                                                                                        |                                                                                                                                        tAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGga                                                                                                                                     <  2:379763/140‑1 (MQ=255)  aaaaGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGgag                                                                                                                                    >  1:33212/1‑140 (MQ=255)      gTCGCAATCCAGCGCCTGAACATCCACCGGACGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACAt                                                                                                                                <  1:707043/140‑1 (MQ=255)        cGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGa                                                                                                                              >  1:1469907/1‑140 (MQ=255)         gCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAg                                                                                                                             <  1:218140/140‑1 (MQ=255)                 gcgcCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGcc                                                                                                                     >  1:1408272/1‑140 (MQ=255)                     cTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCcagca                                                                                                                 <  2:1042602/140‑1 (MQ=255)                       gAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAgg                                                                                                               >  1:432372/1‑140 (MQ=255)                                 cGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTtct                                                                                                     >  1:725571/1‑140 (MQ=255)                                  ggATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTtctc                                                                                                    >  2:796399/1‑140 (MQ=255)                                   gatgatGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCtca                                                                                                   <  2:693756/140‑1 (MQ=255)                                      gatgCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCtcatca                                                                                                <  1:245512/140‑1 (MQ=255)                                       atgCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCaa                                                                                               <  1:1070637/140‑1 (MQ=255)                                       atgCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCaa                                                                                               >  1:640811/1‑140 (MQ=255)                                        tgCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCACGCTGGTGCGCAAACGTGATCATTTCCGCCAGTGGCTTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCaaa                                                                                              >  1:672053/1‑140 (MQ=255)                                          cATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACa                                                                                            >  1:534674/1‑140 (MQ=255)                                               ccGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTa                                                                                       >  2:1181293/1‑140 (MQ=255)                                                  ccTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCTAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGc                                                                                    >  2:1302280/1‑140 (MQ=255)                                                  ccTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGc                                                                                    >  1:1292539/1‑140 (MQ=255)                                                  ccTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGCGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGc                                                                                    >  1:182506/1‑140 (MQ=255)                                                   cTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCa                                                                                   >  2:737031/1‑140 (MQ=255)                                                      aGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCg                                                                                <  2:432372/140‑1 (MQ=255)                                                      aGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCg                                                                                >  1:501318/1‑140 (MQ=255)                                                       gcgcCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGt                                                                               >  1:1571527/1‑140 (MQ=255)                                                        cgcCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGtc                                                                              >  1:318986/1‑140 (MQ=255)                                                             tCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAg                                                                         <  2:521103/140‑1 (MQ=255)                                                                  caGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCa                                                                    >  2:1250520/1‑140 (MQ=255)                                                                          ttGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCa                                                            >  1:1229901/1‑140 (MQ=255)                                                                          ttGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGCAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCa                                                            >  2:1553959/1‑140 (MQ=255)                                                                           tGCGCCATGCTGGTGCGCAAGTGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCAt                                                           <  2:157955/140‑1 (MQ=255)                                                                                cATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGa                                                      <  2:725571/140‑1 (MQ=255)                                                                                cATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGa                                                      <  2:1408272/140‑1 (MQ=255)                                                                                  tGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGAtt                                                    >  1:639248/1‑140 (MQ=255)                                                                                          cgcAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTAGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGa                                            >  1:861839/1‑140 (MQ=255)                                                                                          cgcAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGa                                            >  2:1528302/1‑140 (MQ=255)                                                                                            cAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATc                                          >  2:1180197/1‑140 (MQ=255)                                                                                            cAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATc                                          <  2:672053/140‑1 (MQ=255)                                                                                            cAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATc                                          >  2:28876/1‑140 (MQ=255)                                                                                            cAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATc                                          >  2:158395/1‑140 (MQ=255)                                                                                            cAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATc                                          >  1:1314535/1‑140 (MQ=255)                                                                                            cAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATc                                          >  2:1411692/1‑140 (MQ=255)                                                                                              aGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATcac                                        >  1:1352893/1‑140 (MQ=255)                                                                                               gCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATcaca                                       <  2:33212/140‑1 (MQ=255)                                                                                               gCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATcaca                                       >  1:234430/1‑140 (MQ=255)                                                                                                 gTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACg                                     <  2:1314535/140‑1 (MQ=255)                                                                                                    aTCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAg                                  >  1:489915/1‑140 (MQ=255)                                                                                                      cATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGct                                >  1:636371/1‑140 (MQ=255)                                                                                                      cATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGct                                <  2:404378/140‑1 (MQ=255)                                                                                                      cATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGct                                <  1:104325/140‑1 (MQ=255)                                                                                                               cAGTGGATTTTCTGTGCCAAGCACGCTGGAGACATGGGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCaa                       <  1:1466206/140‑1 (MQ=255)                                                                                                                  tGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAAcgc                    <  2:1319266/140‑1 (MQ=255)                                                                                                                  tGGATTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAAcgc                    <  2:796221/140‑1 (MQ=255)                                                                                                                     aTTTTCTGTGCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTg                 <  1:925371/140‑1 (MQ=255)                                                                                                                              gCCAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCa        >  2:1098388/1‑140 (MQ=255)                                                                                                                                cAAGCACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATc      <  2:489915/140‑1 (MQ=255)                                                                                                                                    cACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGcc  <  1:546026/140‑1 (MQ=255)                                                                                                                                    cACGCTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGcc  <  1:1302280/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                        TAAAAGTCGCAATCCAGCGCCTGAACATCCACCGGATGATGCATCACCGCCTGAGCGCCATCCACCAGCACTTTTGCGCCATGCTGGTGCGCAAGCGTGATCATTTCCGCCAGTGGATTTTCTGTGCCAAGCACGTTGGAGACATGAGTAATTGCCAGCAGGCGAGTTTTCTCATCAAACAGCGTAGGCAGCGTCTCCAGTTGCAACGTACCATCGGGATTGAGCGGGATCACACGCAGCTCTGCGCCAACGCGTGCGCAAAGCATCTGCC  >  NC_000913/1759865‑1760135

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