Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,892,455	G→A	A281T (GCG→ACG)	mdtL →	multidrug efflux system protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,892,455	0	G	A	100.0%	90.4 / NA	29	A281T (GCG→ACG)	mdtL	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (16/13);  total (16/13)

CGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGC  >  NC_000913/3892348‑3892589                                                                                                            |                                                                                                                                       cgcTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCa                                                                                                         <  2:508510/139‑1 (MQ=255)   cTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGACGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATg                                                                                                       >  2:756024/1‑139 (MQ=255)         gCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGttt                                                                                                 <  2:876401/139‑1 (MQ=255)            ggCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtctc                                                                                              <  1:756024/139‑1 (MQ=255)                tCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGtt                                                                                          <  1:50660/139‑1 (MQ=255)                  aGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTg                                                                                        >  2:519350/1‑139 (MQ=255)                       gACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATc                                                                                   >  2:748894/1‑139 (MQ=255)                          cGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACg                                                                                >  1:107880/1‑139 (MQ=255)                                      cACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGcc                                                                    >  1:241488/1‑139 (MQ=255)                                        cGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCgg                                                                  >  1:969215/1‑139 (MQ=255)                                                 cgcTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGt                                                         <  2:930730/139‑1 (MQ=255)                                                   cTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAg                                                       >  2:179764/1‑139 (MQ=255)                                                                    acgtacgtTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATg                                      >  2:209308/1‑139 (MQ=255)                                                                       tacgtTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGt                                   >  2:568096/1‑139 (MQ=255)                                                                                  tCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAgg                        <  1:687167/139‑1 (MQ=255)                                                                                    aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                               >  1:280606/1‑66 (MQ=255)                                                                                    aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                               <  2:280606/66‑1 (MQ=255)                                                                                    aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGc                      >  1:858241/1‑139 (MQ=255)                                                                                    aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGc                      >  1:320822/1‑139 (MQ=255)                                                                                    aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGc                      <  1:46339/139‑1 (MQ=255)                                                                                      cTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCg                    >  2:304924/1‑139 (MQ=255)                                                                                      cTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCg                    <  1:304924/139‑1 (MQ=255)                                                                                        cGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACg                                                                                >  2:271938/1‑77 (MQ=255)                                                                                        cGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACg                                                                                <  1:271938/77‑1 (MQ=255)                                                                                        cGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTACCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGtt                  >  2:222890/1‑139 (MQ=255)                                                                                         gCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGTGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGttt                 >  1:519002/1‑139 (MQ=255)                                                                                               gtTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTa           <  1:396174/139‑1 (MQ=255)                                                                                                   ttCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGTGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAcgcg       <  2:519002/139‑1 (MQ=255)                                                                                                   ttCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAcgcg       <  2:73175/139‑1 (MQ=255)                                                                                                        ggcgacGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGc  >  2:2498/1‑139 (MQ=255)                                                                                                            |                                                                                                                                       CGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGC  >  NC_000913/3892348‑3892589

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