Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,158,566	T→C	S61P (TCT→CCT)	mdtC →	multidrug efflux system, subunit C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,158,566	0	T	C	100.0%	138.9 / NA	42	S61P (TCT→CCT)	mdtC	multidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/21);  total (21/21)

CGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGTCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGG  >  NC_000913/2158423‑2158713                                                                                                                                                |                                                                                                                                                    cggcgaTTTTACCGTCGGTTGCCATTACCCTGTGCGGAATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTAAGCGCGTCGCTGCCCGGTGCGTCACCAGAAGCAATGGCGCCTTCCg                                                                                                                                               <  1:833858/147‑1 (MQ=255)    cgaTTTTACTGTCGGTTGCCATTACCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTg                                                                                                                                            <  1:83063/150‑1 (MQ=255)    cgaTTTTACTGTCGGTTGCCATTACCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGc                                                                                                                                           <  2:274587/151‑1 (MQ=255)     gaTTTTACTGTCGGTTGCCATTACCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGcc                                                                                                                                          >  2:421849/1‑151 (MQ=255)        tttACTGTCGGTTGCCATTACCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCAc                                                                                                                                        <  2:456008/150‑1 (MQ=255)           aCTGTCGGTTGCCATTACCCTGTGCGGCATCCTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCAcgcc                                                                                                                                     <  2:150240/150‑1 (MQ=255)            cTGTCGGTTGCCATTACCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCAcgccgc                                                                                                                                   >  2:610451/1‑151 (MQ=255)              gTCGGTTGCCATTACCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTg                                                                                                                                 >  1:622542/1‑151 (MQ=255)                    tGCCATTACCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAgcgc                                                                                                                           >  1:604913/1‑151 (MQ=255)                           aCCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACtt                                                                                                                     <  2:340081/150‑1 (MQ=255)                            cccTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTgg                                                                                                                   <  2:145994/151‑1 (MQ=255)                             ccTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTggg                                                                                                                  >  2:227529/1‑151 (MQ=255)                               tgtgCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGgcg                                                                                                                >  2:853037/1‑151 (MQ=255)                                gtgCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGgc                                                                                                                 <  2:699492/149‑1 (MQ=255)                                   cGGCATACTGGGCTGCCCTATGGTGCCGCTCGCCCCGCTGCCGCGGTTCGAGTTTCCGGTTATTTTCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAGCAATGGCGCCTTCCGTTTCCACGCCGCTGGAGCGCTCACTTGGgcg                                                                                                                <  1:794276/147‑1 (MQ=255)                                     gCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCGGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGc                                                                                                          <  2:683166/151‑1 (MQ=255)                                                ttCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTg                                                                                               >  2:722365/1‑151 (MQ=255)                                                  ccGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGa                                                                                              <  2:519995/150‑1 (MQ=255)                                                        gctgcCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGAcc                                                                                       <  2:835959/151‑1 (MQ=255)                                                                 cGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGt                                                                              >  2:794265/1‑151 (MQ=255)                                                                           ccGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGgca                                                                    <  1:582349/151‑1 (MQ=255)                                                                                    gATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGt                                                            >  2:146281/1‑150 (MQ=255)                                                                                     aTTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGtat                                                          >  2:334856/1‑151 (MQ=255)                                                                                       tttCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGtatt                                                         >  1:748089/1‑150 (MQ=255)                                                                                       tttCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGtatt                                                         <  2:345397/150‑1 (MQ=255)                                                                                         tCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGtattat                                                       <  1:199108/150‑1 (MQ=255)                                                                                         tCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTattatt                                                      >  2:668663/1‑151 (MQ=255)                                                                                          ccGGTTATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTAttt                                                     >  1:170399/1‑151 (MQ=255)                                                                                               gATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCa                                                 <  2:573088/150‑1 (MQ=255)                                                                                                  tATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGtt                                              <  1:840034/150‑1 (MQ=255)                                                                                                     cGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTg                                            <  2:540102/149‑1 (MQ=255)                                                                                                                cGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCggg                                >  1:845292/1‑150 (MQ=255)                                                                                                                cGCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCgg                                 <  1:35322/149‑1 (MQ=255)                                                                                                                 gCTGCCCGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGat                              <  2:678973/151‑1 (MQ=255)                                                                                                                     cccGGTGCGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCa                          <  1:112958/151‑1 (MQ=255)                                                                                                                            cGTCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGgcgc                   >  2:848337/1‑151 (MQ=255)                                                                                                                              tCACCAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGgcgc                   >  1:420071/1‑149 (MQ=255)                                                                                                                                 ccAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAgcgc              >  1:629964/1‑151 (MQ=255)                                                                                                                                 ccAGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGGATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAgcgc              >  1:548971/1‑151 (MQ=255)                                                                                                                                   aGAAACAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGt            >  1:45542/1‑151 (MQ=255)                                                                                                                                       aCAATGGCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGAtg        >  2:73845/1‑151 (MQ=255)                                                                                                                                            ggCGCCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGc     >  2:666831/1‑149 (MQ=255)                                                                                                                                             gcgcCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAgg  <  2:690766/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                    CGACGATTTTACTGTCGGTTGCCATTACCCTGTGCGGCATACTGGGCTTCCGTATGCTGCCGGTCGCCCCGCTGCCGCAGGTCGATTTTCCGGTGATTATCGTCAGCGCCTCGCTGCCCGGTGCGTCACCAGAAACAATGGCGTCTTCCGTTGCCACGCCGCTGGAGCGCTCACTTGGGCGCATTGCCGGAGTCAGTGAAATGACCTCCAGCAGTTCGCTCGGCAGCACGCGTATTATTTTGCAGTTTGATTTTGACCGGGATATCAACGGCGCAGCGCGTGATGTGCAGG  >  NC_000913/2158423‑2158713

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