Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,474,817	(G)5→4	coding (149/288 nt)	tusB ←	sulfurtransferase complex subunit TusB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,474,813	0	G	.	100.0%	172.2 / NA	41	coding (153/288 nt)	tusB	sulfurtransferase complex subunit TusB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base . (16/25);  total (16/25)

CACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGA  >  NC_000913/3474679‑3474944                                                                                                                                       |                                                                                                                                    cACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATggggca                                                                                                                                <  1:370300/140‑3 (MQ=255)        ccATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCa                                                                                                                        <  2:536780/140‑1 (MQ=255)         cATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAg                                                                                                                       <  1:816676/140‑1 (MQ=255)           tCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGAc                                                                                                                     >  1:487499/1‑140 (MQ=255)            ctggctggGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACt                                                                                                                    >  2:225625/1‑140 (MQ=255)              ggctggGGTGCTTAACCGTAAGTCTGACGAAATTAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACttt                                                                                                                  <  2:389514/140‑1 (MQ=255)              ggctggGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACttt                                                                                                                  <  2:2935/140‑1 (MQ=255)                 tggGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCaa                                                                                                               >  1:820262/1‑140 (MQ=255)                     gTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCa                                                                                                                        >  2:563263/1‑127 (MQ=255)                     gTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCa                                                                                                                        <  1:563263/127‑1 (MQ=255)                     gTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTa                                                                                                           >  2:369361/1‑140 (MQ=255)                         ttAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCgg                                                                                                       >  2:585136/1‑140 (MQ=255)                         ttAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCgg                                                                                                       <  2:69860/140‑1 (MQ=255)                                  aGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTGCGCAGACTTTCAAGGTAGCGGTTACCGTCa                                                                                              >  1:888895/1‑140 (MQ=255)                                           aaaTCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCa                                                                                     <  1:225625/140‑1 (MQ=255)                                                 gTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACg                                                                               <  1:320365/140‑1 (MQ=255)                                                    taGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCa                                                                            <  1:524620/140‑1 (MQ=255)                                                     aGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCAt                                                                           <  1:228338/140‑1 (MQ=255)                                                     aGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCAt                                                                           <  2:536208/140‑1 (MQ=255)                                                       tCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCt                                                                         >  2:677876/1‑140 (MQ=255)                                                       tCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCt                                                                         >  2:171655/1‑140 (MQ=255)                                                                 tgatgTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAataata                                                               >  2:407061/1‑140 (MQ=255)                                                                    tgTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGACTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCa                                                            >  2:53942/1‑140 (MQ=255)                                                                      tCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGt                                                          <  2:63397/140‑1 (MQ=255)                                                                        gTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTGCGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTc                                                        <  2:888895/140‑1 (MQ=255)                                                                         ttCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCg                                                       >  2:158652/1‑140 (MQ=255)                                                                         ttCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCg                                                       >  1:558586/1‑140 (MQ=255)                                                                          tCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGt                                                      <  2:544792/140‑1 (MQ=255)                                                                                 ttGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTc                                               <  1:585136/140‑1 (MQ=255)                                                                                  tGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCa                                              <  2:547118/140‑1 (MQ=255)                                                                                  tGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCa                                              <  1:238480/140‑1 (MQ=255)                                                                                   gACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCAc                                             <  1:115739/140‑1 (MQ=255)                                                                                    aCCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACt                                            <  2:820262/140‑1 (MQ=255)                                                                                        gTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGc                                        >  2:331956/1‑140 (MQ=255)                                                                                                gcgGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACgcag                                <  1:31529/140‑1 (MQ=255)                                                                                                                tCGTTCAGGGCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATcc                <  1:407061/140‑1 (MQ=255)                                                                                                                       gggCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGt              >  2:483374/1‑135 (MQ=255)                                                                                                                       gggCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGt              <  1:483374/135‑1 (MQ=255)                                                                                                                       gggCATAGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGcc         <  1:331956/140‑1 (MQ=255)                                                                                                                             aGACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTg   >  1:700268/1‑140 (MQ=255)                                                                                                                              gACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGa  <  2:487499/140‑1 (MQ=255)                                                                                                                              gACCTTAAT‑GGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGa  <  1:677876/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                    CACCAGGCCATCTGGCTGGGGTGCTTAACCGTAAGTCTGACGAAATCAGTATAGTCAATGAGAATGATGTCGTTCGAAATTTGACCAGTCAAACCGCGGGCAATAAGGTCTTCGTTCAGGGCATAGACCTTAATGGGGGCATTACGCAGACTTTCAAGGTAGCGGTTACCGTCAACTGCGGCAGTTACGCCATCTTGCAATAATAGCAGTTCGTCTCCTTCACTGAGCAGACGCAGCAGCGCAGCAAAATCCGTCAGCCAGGGTGA  >  NC_000913/3474679‑3474944

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