Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	470,643	G→A	*591* (TAG→TAA)  S3N (AGT→AAT)	mdlA → mdlB →	putative multidrug ABC transporter ATPase putative multidrug ABC transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	470,643	0	G	A	100.0%	183.6 / NA	57	*591* (TAG→TAA)  S3N (AGT→AAT)	mdlA mdlB	putative multidrug ABC transporter ATPase putative multidrug ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/28);  total (29/28)

TCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCT  >  NC_000913/470513‑470777                                                                                                                                   |                                                                                                                                       tCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGcc                                                                                                                                >  1:185579/1‑139 (MQ=255)      cAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACtg                                                                                                                           >  1:354314/1‑139 (MQ=255)                        gtgCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAgcgc                                                                                                         <  2:1130450/139‑1 (MQ=255)                           cTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTg                                                                                                      <  2:185579/139‑1 (MQ=255)                                        gCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCGAGCGCCTGTTAGCGTACGGtt                                                                                         >  1:600070/1‑139 (MQ=255)                                        gCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGtt                                                                                         >  1:572268/1‑139 (MQ=255)                                            cTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGcc                                                                                     <  2:600070/139‑1 (MQ=255)                                                    gcgATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTa                                                                             <  1:77585/139‑1 (MQ=255)                                                         atGTATCACTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCg                                                                        >  2:557628/1‑139 (MQ=255)                                                          tGTATCGCTATCAACAACTGGAGGTGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGcc                                                                                     <  1:887197/125‑1 (MQ=255)                                                          tGTATCGCTATCAACAACTGGAGGTGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGcc                                                                                     >  2:887197/1‑125 (MQ=255)                                                                gCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTa                                                                                              <  2:70548/110‑1 (MQ=255)                                                                gCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTa                                                                                              >  1:70548/1‑110 (MQ=255)                                                                    tCACCAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAAtt                                                                                                                                      <  2:740406/66‑1 (MQ=255)                                                                    tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGTGAGGAGGCCGTCGATGCGTAAtt                                                                                                                                      >  1:740406/1‑66 (MQ=255)                                                                    tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                             >  2:325829/1‑139 (MQ=255)                                                                    tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                             >  1:1105820/1‑139 (MQ=255)                                                                    tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                             >  1:956264/1‑139 (MQ=255)                                                                    tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                             >  1:170606/1‑139 (MQ=255)                                                                      aacaacTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           >  1:764379/1‑139 (MQ=255)                                                                      aacaacTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           >  1:376841/1‑139 (MQ=255)                                                                           cTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                                                      <  2:1105820/139‑1 (MQ=255)                                                                            tGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGAGTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtga                                                     >  2:813809/1‑139 (MQ=255)                                                                            tGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAgc                                                                                                           <  1:185346/85‑1 (MQ=255)                                                                            tGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAgc                                                                                                           >  2:185346/1‑85 (MQ=255)                                                                              gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                                                      >  1:610015/1‑136 (MQ=255)                                                                              gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                                                      >  1:106307/1‑136 (MQ=255)                                                                              gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                                                      <  2:106307/136‑1 (MQ=255)                                                                              gAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCtg                                                      <  2:610015/136‑1 (MQ=255)                                                                                   ggcgCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           <  2:460662/126‑1 (MQ=255)                                                                                   ggcgCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           >  1:460662/1‑126 (MQ=255)                                                                                   ggcgCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGatgatg                                                >  1:868670/1‑137 (MQ=255)                                                                                   ggcgCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGatgatg                                                <  2:868670/137‑1 (MQ=255)                                                                                     cgcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGc                                                                                                                       <  2:31572/64‑1 (MQ=255)                                                                                     cgcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGc                                                                                                                       >  1:31572/1‑64 (MQ=255)                                                                                      gcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCaa                                                                                                             >  1:1005161/1‑73 (MQ=255)                                                                                      gcTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCaa                                                                                                             <  2:1005161/73‑1 (MQ=255)                                                                                              acgCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTgcggcgg                                   <  2:376841/139‑1 (MQ=255)                                                                                               cgCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGcc                                                                                                        <  1:308070/69‑1 (MQ=255)                                                                                               cgCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGcc                                                                                                        >  2:308070/1‑69 (MQ=255)                                                                                                 cTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACtctc                                                                                                               >  2:407158/1‑60 (MQ=255)                                                                                                 cTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACtctc                                                                                                               <  1:407158/60‑1 (MQ=255)                                                                                                 cTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGcggcggcgg                                <  2:170606/139‑1 (MQ=255)                                                                                                      gAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGaa                           <  2:572268/139‑1 (MQ=255)                                                                                                            cgcgAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           >  1:999464/1‑101 (MQ=255)                                                                                                            cgcgAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                           <  2:999464/101‑1 (MQ=255)                                                                                                                   aggCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCgc              <  2:354314/139‑1 (MQ=255)                                                                                                                      ccGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAg                      <  2:1032911/128‑1 (MQ=255)                                                                                                                      ccGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAg                      >  1:1032911/1‑128 (MQ=255)                                                                                                                        gTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCtt         <  2:764379/139‑1 (MQ=255)                                                                                                                            aTGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCa     <  2:365473/139‑1 (MQ=255)                                                                                                                             tGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAg    >  1:656442/1‑139 (MQ=255)                                                                                                                               cGTAATTTTAGTCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTa                                                                             <  2:110243/64‑1 (MQ=255)                                                                                                                               cGTAATTTTAGTCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTa                                                                             >  1:110243/1‑64 (MQ=255)                                                                                                                               cGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                                              >  2:284483/1‑79 (MQ=255)                                                                                                                               cGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                                              <  1:284483/79‑1 (MQ=255)                                                                                                                               cGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCt  <  2:656442/139‑1 (MQ=255)                                                                                                                                   |                                                                                                                                       TCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATCAGCT  >  NC_000913/470513‑470777

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