Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,759,757	A→T	I256N (ATC→AAC)	sufS ←	L‑cysteine desulfurase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,759,757	0	A	T	100.0%	125.7 / NA	37	I256N (ATC→AAC)	sufS	L‑cysteine desulfurase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (18/19);  total (18/19)

AGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGC  >  NC_000913/1759628‑1759882                                                                                                                                  |                                                                                                                              aGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAAcc                                                                                                                     <  1:162730/140‑1 (MQ=255)   ccccAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGc                                                                                                                   >  2:81303/1‑140 (MQ=255)      cAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGcccc                                                                                                                >  2:70082/1‑140 (MQ=255)        gcgcCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCt                                                                                                              >  1:69942/1‑140 (MQ=255)            cGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTccc                                                                                                          >  2:183775/1‑140 (MQ=255)                    ccTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACggcggc                                                                                                  <  2:7717/139‑1 (MQ=255)                    aCTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACggcggc                                                                                                  <  2:216220/140‑1 (MQ=255)                      tCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCAt                                                                                                >  2:315769/1‑140 (MQ=255)                        cAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCAtct                                                                                              >  2:127767/1‑140 (MQ=255)                           cgccgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCAt                                                                                                >  1:147546/1‑135 (MQ=255)                           cgccgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCAt                                                                                                <  2:147546/135‑1 (MQ=255)                           cgccgcgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCt                                                                                           <  1:127767/140‑1 (MQ=255)                                cgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCAcggcg                                                                                                    <  1:182096/126‑1 (MQ=255)                                cgccAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCAcggcg                                                                                                    >  2:182096/1‑126 (MQ=255)                                         cAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGcttct                                                                             <  1:81303/140‑1 (MQ=255)                                         cAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGcttct                                                                             <  2:69942/140‑1 (MQ=255)                                                                 tACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCg                                                     >  2:21933/1‑140 (MQ=255)                                                                  aCCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCgg                                                    <  2:78629/140‑1 (MQ=255)                                                                  aCCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCgg                                                    >  1:78629/1‑140 (MQ=255)                                                                       cTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTgggg                                               <  1:21933/140‑1 (MQ=255)                                                                       cTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTgggg                                               <  1:82942/140‑1 (MQ=255)                                                                       cTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTgggg                                               <  1:315769/140‑1 (MQ=255)                                                                       cTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTgggg                                               <  1:315552/140‑1 (MQ=255)                                                                        ttCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTg                                                  <  1:299526/136‑1 (MQ=255)                                                                        ttCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTg                                                  >  2:299526/1‑136 (MQ=255)                                                                         tCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGTTTCTTTCACATAAAGAATGCCAATTCCGGTg                                                  >  2:351726/1‑135 (MQ=255)                                                                         tCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGTTTCTTTCACATAAAGAATGCCAATTCCGGTg                                                  <  1:351726/135‑1 (MQ=255)                                                                                 ccATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGtt                                     <  1:111360/140‑1 (MQ=255)                                                                                      gTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGc                                <  2:322008/140‑1 (MQ=255)                                                                                         cTTTGGTCCAGGTAGTGCCCTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCg                             <  1:70082/140‑1 (MQ=255)                                                                                           ttGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGga                           <  2:315801/140‑1 (MQ=255)                                                                                                             tCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCa         >  1:108243/1‑140 (MQ=255)                                                                                                             tCACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCa         >  2:96427/1‑140 (MQ=255)                                                                                                              cACTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCaa        >  1:67231/1‑140 (MQ=255)                                                                                                               aCTCAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAAt       >  1:199529/1‑140 (MQ=255)                                                                                                                  cAGGCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCa    >  1:289673/1‑140 (MQ=255)                                                                                                                    ggCTGACGGTGGCGTTCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCTCGGAGAACACGTAAAAGTCGCAATCCAgc  >  1:148794/1‑140 (MQ=255)                                                                                                                                  |                                                                                                                              AGCCCCAGCGCCGAAACATACTCCAGCGCCGCGCCAAGACCAATGATGCCCCCGGTATTGGGTGTACCGGCTTCAAACCGCCATGGTGCTTTGGTCCAGGTAGTGCCTTCACTCAGGCTGACGGTGGCGATCATAGAACCGCCCCCTTCCCACGGCGGCATCTCCTGCAACAAGGCTTCTTTCACATAAAGAATGCCAATTCCGGTGGGGCCATACAGTTTATGCCCGGAGAACACGTAAAAGTCGCAATCCAGC  >  NC_000913/1759628‑1759882

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