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%	105.9 / NA	34	*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 (18/16);  total (18/16)

TATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATC  >  NC_000913/470511‑470773                                                                                                                                     |                                                                                                                                   tatCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAg                                                                                                                              <  1:15387/139‑1 (MQ=255)   tCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGcc                                                                                                                            >  1:786742/1‑139 (MQ=255)   tCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCCTAATTTTAGcc                                                                                                                            >  2:77448/1‑139 (MQ=255)        cAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACtg                                                                                                                       >  2:608898/1‑139 (MQ=255)           cGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTgg                                                                                                                    >  2:725886/1‑139 (MQ=255)              ggCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCg                                                                                                                 <  1:725886/139‑1 (MQ=255)                           tgCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGcc                                                                                                    <  1:563241/139‑1 (MQ=255)                             cTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTg                                                                                                  >  2:79227/1‑139 (MQ=255)                             cTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTg                                                                                                  >  2:636364/1‑139 (MQ=255)                              tGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGt                                                                                                 >  2:443707/1‑139 (MQ=255)                                   caacaaAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCg                                                                                            >  2:231404/1‑139 (MQ=255)                                          gCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGTGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTg                                                                                                  <  2:427652/126‑1 (MQ=255)                                          gCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGTGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTg                                                                                                  >  1:427652/1‑126 (MQ=255)                                                ggTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGt                                                                               >  2:233911/1‑139 (MQ=255)                                                             gTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCt                                                                  <  1:636364/139‑1 (MQ=255)                                                              tATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTg                                                                 >  1:597047/1‑139 (MQ=255)                                                                    tATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGAtt                                                           <  2:203667/139‑1 (MQ=255)                                                                     aTCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTg                                                          <  1:231404/139‑1 (MQ=255)                                                                      tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                         >  2:186753/1‑139 (MQ=255)                                                                      tCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGc                                                         <  1:103042/139‑1 (MQ=255)                                                                        aacaacTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCgg                                                       >  2:671677/1‑139 (MQ=255)                                                                          caacTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTc                                                     <  1:671677/139‑1 (MQ=255)                                                                                  ggcggcGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTgatgat                                             >  2:170062/1‑139 (MQ=255)                                                                                   gcggcgATAGACGGCGCTCCGGAAAATCGCGCGTAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACCGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGatgatg                                            <  1:568512/139‑1 (MQ=255)                                                                                                   cTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGATTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGcggcggcgg                            >  1:163341/1‑139 (MQ=255)                                                                                                        gAAAATCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGaa                       <  1:186753/139‑1 (MQ=255)                                                                                                             tCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTa                                                                                          >  1:612338/1‑67 (MQ=255)                                                                                                             tCGCGAGGAGGCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTa                                                                                          <  2:612338/67‑1 (MQ=255)                                                                                                                  aggaggCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCATACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGc             <  1:608898/139‑1 (MQ=255)                                                                                                                   ggaggCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTggg                                                               <  2:293423/88‑1 (MQ=255)                                                                                                                   ggaggCCGTCGATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTggg                                                               >  1:293423/1‑88 (MQ=255)                                                                                                                             gATGCGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATc  <  1:79227/139‑1 (MQ=255)                                                                                                                                 cGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCg                                                                                >  2:58896/1‑57 (MQ=255)                                                                                                                                 cGTAATTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCg                                                                                <  1:58896/57‑1 (MQ=255)                                                                                                                                     |                                                                                                                                   TATCGCCCAGCGTGGCAATCATGATGTGCTGGCACAACAAAGCGGCTGGTATCGCGATATGTATCGCTATCAACAACTGGAGGCGGCGCTCGACGACGCTCCGGAAAATCGCGAGGAGGCCGTCGATGCGTAGTTTTAGCCAACTGTGGCCGACTCTCAAGCGCCTGTTAGCGTACGGTTCGCCGTGGCGTAAACCGCTGGGGATTGCGGTCCTGATGATGTGGGTTGCGGCGGCGGCAGAAGTCAGTGGGCCGCTGCTTATC  >  NC_000913/470511‑470773

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