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

GGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGC  >  NC_000913/3892326‑3892589                                                                                                                                  |                                                                                                                                       ggTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCa                                                                                                                               <  2:605106/139‑1 (MQ=255)  gTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCAc                                                                                                                              >  1:239481/1‑139 (MQ=255)               caTTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCa                                                                                                                 >  1:672777/1‑139 (MQ=255)               caTTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCa                                                                                                                 >  2:432759/1‑139 (MQ=255)               caTTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCa                                                                                                                 >  1:15361/1‑139 (MQ=255)               caTTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCa                                                                                                                 <  2:1449604/139‑1 (MQ=255)                 ttATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGGCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCAcc                                                                                                               <  1:421186/139‑1 (MQ=255)                 ttATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCAcc                                                                                                               >  2:1390530/1‑139 (MQ=255)                         cTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATg                                                                                                       <  1:477818/139‑1 (MQ=255)                            aCCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCgg                                                                                                    <  2:261149/139‑1 (MQ=255)                            aCCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCgg                                                                                                    >  1:1295446/1‑139 (MQ=255)                            aCCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCgg                                                                                                    >  2:1319322/1‑139 (MQ=255)                                      tCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTcc                                                                                                           >  2:283240/1‑122 (MQ=255)                                      tCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTcc                                                                                                           <  1:283240/122‑1 (MQ=255)                                             gACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATc                                                                                   >  2:1320764/1‑139 (MQ=255)                                                  tttCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCa                                                                                                                 >  2:1064430/1‑104 (MQ=255)                                                  tttCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCa                                                                                                                 <  1:1064430/104‑1 (MQ=255)                                                      aTTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGAtt                                                                          <  2:1295446/139‑1 (MQ=255)                                                           ccACGCCATATGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTgcgc                                                                     >  2:986270/1‑139 (MQ=255)                                                            cACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCActct                                                                                                                           <  1:498062/84‑1 (MQ=255)                                                            cACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCActct                                                                                                                           >  2:498062/1‑84 (MQ=255)                                                            cACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGcc                                                                    >  2:927944/1‑139 (MQ=255)                                                                   tttGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTtct                                                             <  1:1390530/139‑1 (MQ=255)                                                                      gcgcTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCgg                                                          <  1:927944/139‑1 (MQ=255)                                                                           gggAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGtt                                                                                          <  2:291216/102‑1 (MQ=255)                                                                           gggAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGtt                                                                                          >  1:291216/1‑102 (MQ=255)                                                                           gggAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTCGCGCCGGTTTCTCGGTAGGt                                                     >  1:625170/1‑139 (MQ=255)                                                                                        ccACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCCCCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGa                                        >  1:1232073/1‑139 (MQ=255)                                                                                                gtTGATGATCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAg                                >  2:612581/1‑139 (MQ=255)                                                                                                    atgatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGt                            >  1:1407999/1‑139 (MQ=255)                                                                                                     tgatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCa                                 <  1:327326/133‑1 (MQ=255)                                                                                                     tgatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCa                                 >  2:327326/1‑133 (MQ=255)                                                                                                      gatCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTa                          >  1:270981/1‑139 (MQ=255)                                                                                                       atCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGc                                                                                                      >  2:251292/1‑62 (MQ=255)                                                                                                       atCACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGc                                                                                                      <  1:251292/62‑1 (MQ=255)                                                                                                         cACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCgg                                                                                                    <  1:788513/62‑1 (MQ=255)                                                                                                         cACCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCgg                                                                                                    >  2:788513/1‑62 (MQ=255)                                                                                                          aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATg                                                                                                       >  1:1455137/1‑58 (MQ=255)                                                                                                          aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATg                                                                                                       <  2:1455137/58‑1 (MQ=255)                                                                                                          aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGttt                                                                                                 <  2:1280463/64‑1 (MQ=255)                                                                                                          aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGttt                                                                                                 >  1:1280463/1‑64 (MQ=255)                                                                                                          aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                               >  2:587856/1‑66 (MQ=255)                                                                                                          aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                               <  1:587856/66‑1 (MQ=255)                                                                                                          aCCTCGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGc                      >  1:56047/1‑139 (MQ=255)                                                                                                              cGCAGGTGTTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGtt                  >  2:402441/1‑139 (MQ=255)                                                                                                                    tgtTATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTCGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCAtt            <  2:625170/139‑1 (MQ=255)                                                                                                                      ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGa                                        >  1:420056/1‑109 (MQ=255)                                                                                                                      ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGa                                        <  2:420056/109‑1 (MQ=255)                                                                                                                      ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAc          >  2:695678/1‑139 (MQ=255)                                                                                                                      ttATTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAc          >  1:1457085/1‑139 (MQ=255)                                                                                                                        aTTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCt                                                                                                              <  2:543192/37‑1 (MQ=255)                                                                                                                        aTTCCTGGCGACGGGGATCACTCTTGCCGTTTCACCt                                                                                                              >  1:543192/1‑37 (MQ=255)                                                                                                                         ttCCTGGCGACGGGTATCACTCTTGCCGTTTCAACTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAcgcg       <  2:1457085/139‑1 (MQ=255)                                                                                                                         ttCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAcgcg       <  2:56047/139‑1 (MQ=255)                                                                                                                         ttCCTGGCGACGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAcgcg       <  1:695678/139‑1 (MQ=255)                                                                                                                              ggcgacGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGttt                                                                                                 <  2:1090802/44‑1 (MQ=255)                                                                                                                              ggcgacGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGttt                                                                                                 >  1:1090802/1‑44 (MQ=255)                                                                                                                              ggcgacGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                               <  2:645111/46‑1 (MQ=255)                                                                                                                              ggcgacGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTtct                                                                                               >  1:645111/1‑46 (MQ=255)                                                                                                                              ggcgacGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGCCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGATCGGGc  >  1:411772/1‑139 (MQ=255)                                                                                                                                  |                                                                                                                                       GGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGC  >  NC_000913/3892326‑3892589

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