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%	116.6 / NA	37	A281T (GCG→ACG)	mdtL	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/18);  total (19/18)

GCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCG  >  NC_000913/3892322‑3892584                                                                                                                                      |                                                                                                                                  gcgcGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACgggg                                                                                                                              <  2:1054835/139‑1 (MQ=255)   gcgGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGAt                                                                                                                            <  1:996902/139‑1 (MQ=255)      gTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGCTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCAc                                                                                                                         >  1:705329/1‑139 (MQ=255)                   caTTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCa                                                                                                            >  1:386733/1‑139 (MQ=255)                     ttATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGgcgacg                                                                                                                                 >  1:411102/1‑116 (MQ=255)                            gcTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCAt                                                                                                   >  2:428846/1‑139 (MQ=255)                             cTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGCTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATg                                                                                                  <  2:705329/139‑1 (MQ=255)                                   gCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGttt                                                                                            <  1:206434/139‑1 (MQ=255)                                    cTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtc                                                                                           >  1:785847/1‑139 (MQ=255)                                              cATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                          >  2:932210/1‑130 (MQ=255)                                              cATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                          <  1:932210/130‑1 (MQ=255)                                                    cGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACg                                                                           <  1:428846/139‑1 (MQ=255)                                                        tCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGtt                                                                                                               <  1:1057363/99‑1 (MQ=255)                                                        tCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGtt                                                                                                               >  2:1057363/1‑99 (MQ=255)                                                                       tttGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTtct                                                        <  1:393232/139‑1 (MQ=255)                                                                        ttGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACg                                                                           >  2:24333/1‑119 (MQ=255)                                                                        ttGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACg                                                                           <  1:24333/119‑1 (MQ=255)                                                                                            ccACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGa                                   <  2:386733/139‑1 (MQ=255)                                                                                                        atgatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGt                       >  1:138050/1‑139 (MQ=255)                                                                                                         tgatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGttt                                              <  1:1040136/115‑1 (MQ=255)                                                                                                         tgatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGttt                                              >  2:1040136/1‑115 (MQ=255)                                                                                                          gatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTa                     >  1:1747/1‑139 (MQ=255)                                                                                                              aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCActct                                                                                                                      >  1:798209/1‑38 (MQ=255)                                                                                                              aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCActct                                                                                                                      <  2:798209/38‑1 (MQ=255)                                                                                                              aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGtt                                                                                     >  1:196217/1‑71 (MQ=255)                                                                                                              aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGtt                                                                                     <  2:196217/71‑1 (MQ=255)                                                                                                              aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCg                                                      >  1:612179/1‑102 (MQ=255)                                                                                                              aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCg                                                      <  2:612179/102‑1 (MQ=255)                                                                                                              aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGc                 >  1:91026/1‑139 (MQ=255)                                                                                                              aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGc                 >  2:246185/1‑139 (MQ=255)                                                                                                              aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGc                 >  1:13837/1‑139 (MQ=255)                                                                                                                cTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCtt                                                                                                                     <  2:212785/37‑1 (MQ=255)                                                                                                                cTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCtt                                                                                                                     >  1:212785/1‑37 (MQ=255)                                                                                                                cTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCg                                                      <  2:434700/100‑1 (MQ=255)                                                                                                                cTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCg                                                      >  1:434700/1‑100 (MQ=255)                                                                                                                          ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATg                                                                                                  >  1:946514/1‑46 (MQ=255)                                                                                                                          ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATg                                                                                                  <  2:946514/46‑1 (MQ=255)                                                                                                                             ttCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAcgcg  <  2:91026/139‑1 (MQ=255)                                                                                                                                  ggcgacGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                          <  2:542996/46‑1 (MQ=255)                                                                                                                                  ggcgacGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                          >  1:542996/1‑46 (MQ=255)                                                                                                                                      |                                                                                                                                  GCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCG  >  NC_000913/3892322‑3892584

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