Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	994,795	(G)5→6	coding (38/792 nt)	ssuC ←	aliphatic sulfonate ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	994,790	1	.	G	100.0%	125.3 / NA	34	W15R (TGG→CGG)	ssuC	aliphatic sulfonate ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (16/18);  total (16/18)

TGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCA‑GGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAA  >  NC_000913/994654‑994922                                                                                                                                          |                                                                                                                                     tgcCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCA‑ggg                                                                                                                                   >  1:943485/1‑140 (MQ=255)  gcCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCA‑gggg                                                                                                                                  >  1:598714/1‑140 (MQ=255)  gcCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCA‑gggg                                                                                                                                  >  2:205526/1‑140 (MQ=255)         ttCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAAcg                                                                                                                            <  1:887054/140‑1 (MQ=255)         tccgccgcTGGCGGAGAGCGTCCAGAACGCCGTCCCCCCTCCTCCCGGGGAAGCCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAAcg                                                                                                                            <  2:598714/138‑1 (MQ=255)                  ggCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAAc                                                                                                                   <  1:959627/140‑1 (MQ=255)                   gCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAAcc                                                                                                                  >  2:287008/1‑140 (MQ=255)                      gagaGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCAct                                                                                                               >  1:822711/1‑140 (MQ=255)                           cGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCActctgtc                                                                                                          >  1:381045/1‑135 (MQ=255)                            gTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCa                                                                                                         <  2:400068/140‑1 (MQ=255)                             tCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATa                                                                                                                     <  2:826668/127‑1 (MQ=255)                             tCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATa                                                                                                                     >  1:826668/1‑127 (MQ=255)                             tCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCAc                                                                                                        >  2:655351/1‑140 (MQ=255)                                  aaCGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCg                                                                                                   <  1:248073/140‑1 (MQ=255)                                    cGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGtt                                                                                                 <  2:964955/140‑1 (MQ=255)                                                    tCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGt                                                                                                  >  1:344985/1‑123 (MQ=255)                                                    tCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGt                                                                                                  <  2:344985/123‑1 (MQ=255)                                                    tCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCt                                                                                 >  2:188042/1‑140 (MQ=255)                                                     ccGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCtt                                                                                >  1:73450/1‑140 (MQ=255)                                                        gTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAg                                                                              <  1:886036/139‑1 (MQ=255)                                                        gTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAg                                                                              >  2:886036/1‑139 (MQ=255)                                                            aGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTg                                                                         >  1:617314/1‑140 (MQ=255)                                                            aGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTg                                                                         <  1:188042/140‑1 (MQ=255)                                                             ggCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTgc                                                                        >  2:499356/1‑140 (MQ=255)                                                             ggCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTgc                                                                        >  1:773854/1‑140 (MQ=255)                                                             ggCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTgc                                                                        >  1:85494/1‑140 (MQ=255)                                                              gCAAAATCCGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTgcg                                                                       <  2:751124/140‑1 (MQ=255)                                                                        cgTCGATAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACg                                                             >  1:536749/1‑140 (MQ=255)                                                                        cgTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACg                                                             <  1:287008/140‑1 (MQ=255)                                                                        cgTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACg                                                             <  1:1153801/140‑1 (MQ=255)                                                                                           gaggagGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGcc                                          <  2:943485/140‑1 (MQ=255)                                                                                            aggaggCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCa                                         <  2:73450/140‑1 (MQ=255)                                                                                                      ttGCCACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGcc                               >  2:681824/1‑140 (MQ=255)                                                                                                         ccACACCGCCACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTg                            <  2:85494/140‑1 (MQ=255)                                                                                                                  cACGATGCCCACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAg                   <  1:361560/140‑1 (MQ=255)                                                                                                                           cACCGGTAAAAACCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAgcggctgcgg          <  2:822711/140‑1 (MQ=255)                                                                                                                                   aaaaCCAGGGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGaa  >  1:747021/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                     TGCCACAGTTCGCCGCTGGCGGAGAGCGTCCAGAACGCCGTCACCACTCCTTCCGGTGAAGGCAAAATACGCGTCGACAGCCAGCCAACCGAGGAGGCCAGTTGCCACACCGCCACGATGCCCACCGGTAAAAACCA‑GGGGGCAACGCGCAATAACCACTTCTTCACTGGCGTTGCCATTATTTCTTCCTTAGCTTTGCGCGACTTTACGGGGGATAAAATCATTCGCCACCGCTTCGCCTTGCGGATTCAGCGGCTGCGGCTGGGGAA  >  NC_000913/994654‑994922

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