Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	994,244	T→G	I197L (ATT→CTT)	ssuC ←	aliphatic sulfonate ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	994,244	0	T	G	100.0%	165.3 / NA	44	I197L (ATT→CTT)	ssuC	aliphatic sulfonate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (30/14);  total (30/14)

ACAGCTGCGCGCTGACGTCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAATCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTCCACCAACCCGCGATCGAT  >  NC_000913/994100‑994385                                                                                                                                                 |                                                                                                                                              aCAGCTGCGCGCTGACGTCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGt                                                                                                                                         >  1:154909/1‑151 (MQ=255)       gcgcgcTGACGTCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGc                                                                                                                                     <  1:733604/149‑1 (MQ=255)         gcgcTGACGTCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCcac                                                                                                                                  >  2:401644/1‑150 (MQ=255)         gcgcTGACGTCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCcac                                                                                                                                  >  2:272435/1‑150 (MQ=255)          cgcTGACGTCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACAt                                                                                                                                >  2:513480/1‑151 (MQ=255)           gcTGACGTCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATc                                                                                                                               >  2:15634/1‑151 (MQ=255)                 gTCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGccc                                                                                                                         >  2:276806/1‑151 (MQ=255)                  tCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGccc                                                                                                                         >  1:311053/1‑150 (MQ=255)                   cAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCa                                                                                                                        >  2:905429/1‑150 (MQ=255)                       cAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGc                                                                                                                    >  1:908877/1‑150 (MQ=255)                        aGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCa                                                                                                                   >  2:624177/1‑150 (MQ=255)                           ttGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAac                                                                                                               >  1:83472/1‑151 (MQ=255)                                     gggCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCCTCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGAc                                                                                                     <  1:502812/151‑1 (MQ=255)                                      ggCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGAcc                                                                                                    >  2:824530/1‑151 (MQ=255)                                          tAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATa                                                                                                 >  1:820115/1‑150 (MQ=255)                                                            accaccACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACca                                                                               <  2:222814/150‑1 (MQ=255)                                                             ccaccaCGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACcagg                                                                             >  2:369755/1‑151 (MQ=255)                                                             ccaccaCGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACcagg                                                                             >  1:37727/1‑151 (MQ=255)                                                                 cacGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACcaggcagg                                                                         <  2:607663/151‑1 (MQ=255)                                                                    gTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATc                                                                      >  1:878309/1‑151 (MQ=255)                                                                           tGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATgg                                                                <  2:703872/150‑1 (MQ=255)                                                                             cAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGAt                                                              >  1:185630/1‑150 (MQ=255)                                                                                 aaCTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACa                                                         <  2:769292/151‑1 (MQ=255)                                                                                             ttCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCgg                                             >  1:806848/1‑151 (MQ=255)                                                                                                  cGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCATCAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAAt                                        <  2:362633/151‑1 (MQ=255)                                                                                                       gATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCAt                                    >  1:452387/1‑150 (MQ=255)                                                                                                       gATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGGATGCCGGATAATCCATa                                   <  2:590018/151‑1 (MQ=255)                                                                                                        aTAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCAt                                    >  1:403484/1‑149 (MQ=255)                                                                                                         tAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGc                                 >  2:886067/1‑151 (MQ=255)                                                                                                           aCCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTa                               <  1:229741/151‑1 (MQ=255)                                                                                                           aCCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTa                               <  2:779825/151‑1 (MQ=255)                                                                                                            ccAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTAc                              >  2:737405/1‑151 (MQ=255)                                                                                                            ccAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTAc                              >  2:90455/1‑151 (MQ=255)                                                                                                               aTGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTAcgcg                           <  2:683875/151‑1 (MQ=255)                                                                                                                tGCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACgcgc                          >  1:270516/1‑151 (MQ=255)                                                                                                                 gCCTGAATTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGcc                         >  1:397904/1‑151 (MQ=255)                                                                                                                       aTTGGCAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTcc                   >  2:716555/1‑151 (MQ=255)                                                                                                                            cAGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTCCACCaa              <  2:803415/151‑1 (MQ=255)                                                                                                                             aGAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTCCACCAAc             >  2:1771/1‑151 (MQ=255)                                                                                                                              gAAATGGTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTCCACCAAcc            >  1:30686/1‑151 (MQ=255)                                                                                                                                    gTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTCCACCAACCCGcgat      <  1:94903/151‑1 (MQ=255)                                                                                                                                    gTTTCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTCCACCAACCCGcgat      >  1:18792/1‑151 (MQ=255)                                                                                                                                       tCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTCCACCAACCCGcgatcg    <  2:154908/150‑1 (MQ=255)                                                                                                                                       tCGGCAACAAGCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGA‑CATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTCCACCAACCCGcgatcgat  >  2:324335/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                              ACAGCTGCGCGCTGACGTCAGCCAGTTTGCCGAGCAGGGCGTAAAGAATAATAGCGACTACCACCACGTCCGTTTGCAAAAACTCCCGCGCATTCATCGCCAGATAACCAATGCCTGAATTGGCAGAAATGGTTTCGGCAACAATCAGCGTCAGCCACATCAGGCCCAACGCAAAACGCACGCCGACCATAATTGAGGGCAGGGCACCAGGCAGGATCACATGGATAAACAGCGGTATGCCGGATAATCCATAGCTACGCGCCATCTCCACCAACCCGCGATCGAT  >  NC_000913/994100‑994385

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