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

GCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGG  >  NC_000913/3892324‑3892588                                                                                                                                    |                                                                                                                                      gcgGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGAt                                                                                                                                <  2:157054/139‑1 (MQ=255)         tACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTg                                                                                                                        >  2:8687/1‑139 (MQ=255)                 caTTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCa                                                                                                                >  2:421071/1‑139 (MQ=255)                 caTTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCa                                                                                                                >  1:1294094/1‑139 (MQ=255)                           cTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATg                                                                                                      >  2:983001/1‑139 (MQ=255)                             gACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTGTTAAGCCACGTACGTTGATGATCACCCCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCg                                                                                                    <  2:1391164/139‑1 (MQ=255)                              aCCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACGTTCCCATGCgg                                                                                                   >  1:935328/1‑139 (MQ=255)                              aCCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCgg                                                                                                   >  2:282589/1‑139 (MQ=255)                                 gCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGttt                                                                                                <  2:279910/139‑1 (MQ=255)                                   tGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                              <  2:935328/139‑1 (MQ=255)                                   tGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                              <  1:983001/139‑1 (MQ=255)                                        tCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGtt                                                                                         >  1:707160/1‑139 (MQ=255)                                            cATGACCGTTTCATTCTCCACGCCATTTGTGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGt                                                                                     >  1:408175/1‑139 (MQ=255)                                             aTGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTa                                                                                    <  2:758020/139‑1 (MQ=255)                                                                cGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGcc                                                                   >  1:1074183/1‑137 (MQ=255)                                                                cGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGcc                                                                   <  2:1074183/137‑1 (MQ=255)                                                                     tttGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTtct                                                            <  1:421071/139‑1 (MQ=255)                                                                        gcgcTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCgg                                                         <  1:315733/139‑1 (MQ=255)                                                                         cgcTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGt                                                        <  1:1121762/139‑1 (MQ=255)                                                                         cgcTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGt                                                        <  1:8687/139‑1 (MQ=255)                                                                           cTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAg                                                      >  1:1057494/1‑139 (MQ=255)                                                                                         gCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCg                                        <  2:985207/139‑1 (MQ=255)                                                                                            acgtacgtTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATg                                     <  2:1294094/139‑1 (MQ=255)                                                                                             cgtacgtTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTctct                                                                                            >  1:460748/1‑83 (MQ=255)                                                                                             cgtacgtTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTctct                                                                                            <  2:460748/83‑1 (MQ=255)                                                                                             cgtacgtTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGA‑CACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAg                                   >  1:1364901/1‑139 (MQ=255)                                                                                               tacgtTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTgcgc                                                                    >  1:119016/1‑105 (MQ=255)                                                                                               tacgtTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTgcgc                                                                    <  2:119016/105‑1 (MQ=255)                                                                                                     gatgatCACCTCGCAGGTGTTATTCCTGGCGACTGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGt                                                        <  2:907926/111‑1 (MQ=255)                                                                                                     gatgatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGt                                                        >  1:907926/1‑111 (MQ=255)                                                                                                      atgatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGt                           >  1:1193463/1‑139 (MQ=255)                                                                                                       tgatCACCTCGCAGGTGTTATTCCTGGCGACGGGGA‑CACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTa                         <  2:1364901/139‑1 (MQ=255)                                                                                                           cACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGcc                                                                   >  2:1159460/1‑94 (MQ=255)                                                                                                           cACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGcc                                                                   <  1:1159460/94‑1 (MQ=255)                                                                                                            aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGgtg                                             <  2:1325009/115‑1 (MQ=255)                                                                                                            aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGgtg                                             >  1:1325009/1‑115 (MQ=255)                                                                                                            aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGc                     <  1:464912/139‑1 (MQ=255)                                                                                                              cTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCg                   >  1:1362557/1‑139 (MQ=255)                                                                                                              cTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCg                   >  1:280836/1‑139 (MQ=255)                                                                                                              cTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCg                   >  1:816936/1‑139 (MQ=255)                                                                                                                cGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCt                                                                                                             >  2:453945/1‑47 (MQ=255)                                                                                                                cGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCt                                                                                                             <  1:453945/47‑1 (MQ=255)                                                                                                                cGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGt                                                                >  1:1352860/1‑92 (MQ=255)                                                                                                                cGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGt                                                                <  2:1352860/92‑1 (MQ=255)                                                                                                                cGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGtt                 >  1:647482/1‑139 (MQ=255)                                                                                                                        ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtctc                                                                                             <  1:678100/55‑1 (MQ=255)                                                                                                                        ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtctc                                                                                             >  2:678100/1‑55 (MQ=255)                                                                                                                        ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                              >  1:1202588/1‑54 (MQ=255)                                                                                                                        ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                              <  2:1202588/54‑1 (MQ=255)                                                                                                                        ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGgt                                              <  1:1133911/102‑1 (MQ=255)                                                                                                                        ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGgt                                              >  2:1133911/1‑102 (MQ=255)                                                                                                                        ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAc         >  2:323119/1‑139 (MQ=255)                                                                                                                          aTTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCg                                                                                                    <  2:884047/46‑1 (MQ=255)                                                                                                                          aTTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCg                                                                                                    >  1:884047/1‑46 (MQ=255)                                                                                                                           ttCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAcgcg      <  2:280836/139‑1 (MQ=255)                                                                                                                           ttCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAcgcg      <  1:282589/139‑1 (MQ=255)                                                                                                                               tGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCggg  >  2:91586/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                      GCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGG  >  NC_000913/3892324‑3892588

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: ‑